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FOREWORD 


In previous years, the National Space Science Data Center (NSSDC) 
has generated a Catalog of Correlative Data. This catalog listed and, 
to some extent, described the correlative data holdings of NSSDC. How- 
ever, it was felt desirable to have a document of a more general nature, 
which would point to correlative data available at any data center or 
in any readily accessible publication. This Handbook of Correlative 
Data was generated to fulfill this need. 

This Handbook can be used to acquaint the reader with many solar 
and geophysical phenomena, to point the reader to more detailed discus- 
sions of the phenomena, and to direct him to the data. It is antici- 
pated that many readers will require this Handbook for only the last 
one or two of these purposes . 

A second edition of this Handbook will be generated when warranted. 
Accordingly, I would appreciate receiving comments from the users of 
this Handbook . Any recommendations for improving its usefulness will 
be most helpful in the generation of a second edition. 

I want to thank the NSSDC personnel who served as authors of the 
individual sections; their names appear on the pages introducing each 
of the nine major parts. I would like to acknowledge the assistance 
of other NSSDC personnel, particularly Miss Mary Anne Beeler, who were 
involved in the various stages of the generation of this Handbook . I 
would also like to thank scientists outside the NSSDC, S. J. Bauer, 

J. E. Jackson, G. D. Mead, and A. J. Zmuda, who reviewed certain sec- 
tions of this Handbook. 


January 15, 1971 
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INTRODUCTION 


The National Space Science Data Center (NSSDC) is issuing this 
Handbook of Correlative Data to inform scientists of the availability 
of data potentially useful as correlative data in space science studies. 
Although the major activities of NSSDC are concerned with the acquisi- 
tion, storage, retrieval, and dissemination of space science data, the 
significance of correlative data in many studies is recognized, and 
assistance in this area is therefore given to investigators . 

For the purposes of this document, correlative data are those data 
that are routinely gathered and that can be used to correlate with space 
science data. 

This Handbook of Correlative Data contains six major discipline 
oriented parts covering galactic cosmic rays, solar electromagnetic 
radiation, energetic solar protons, geomagnetism, the ionosphere, and 
the neutral atmosphere. In addition, a miscellaneous data part covers 
magnetospherically trapped particles, solar wind, airglow, aurora, 
calendar records, activity charts, and Jovian radio emission. Each 
section includes: 

• a brief description of the phenomenon 

• reference to more extensive discussions of the phenomenon 

• reference to discussion of measurement techniques 

• a brief discussion of available data 

• the time periods for which data exist 

• the medium on which data are stored 

• sources of more extensive data availability listings, if 
appropriate 

• sources from which the actual data can be obtained (both 
places and publications) 

In the cases of geomagnetism, neutral atmosphere, and magnetospherically 
trapped particles, models based on syntheses of large amounts of data 
are discussed. 
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Following the miscellaneous data part is a miscellaneous information 
part in which various coordinate systems useful in space science, the 
International Ursigram and World Days Service, and rocket data availa- 
bility are discussed. These fall outside the definition of correlative 
data, but were judged to be of sufficient usefulness to warrant their 
inclusion. 

All parts except for Part IX are intended to stand alone, and each 
has a separate list of references . 

The final part of the Handbook lists and discusses the major sources 
of correlative data. Also listed are all other documents referred to in 
the text as containing some correlative data. Sufficient detail is given 
on all sources to enable the reader to obtain the document and the data. 

If the sources cannot supply the needed data, NSSDC will endeavor to 
assist users in locating the data. It should be pointed out here that 
no attempt has been made to survey or reference the data catalogs of indi- 
vidual observatories. Station lists are included in the WDC-A Catalogue 
of Data on Solar-Terrestrial Physics . The World Data Center A for Upper 
Atmosphere Geophysics can provide the mailing addresses of the individual 
observatories as required. 

It is the policy of NSSDC to obtain and provide correlative data for 
the use of other NASA personnel or contractors. NSSDC visitors are welcome 
to utilize any part of the NSSDC correlative data holdings. All inquiries 
concerning correlative data or requests for copies of this Handbook or 
for the NSSDC Data Catalog of Satellite Experiments , which describes the 
available spacecraft space science data, should be sent to: 

National Space Science Data Center 
Code 601.4 

Goddard Space Flight Center 
Greenbelt, Maryland 20771 

Phone (301) 982-6695 
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INTRODUCTION 


The measurement of cosmic-ray intensity by ground-based monitors 
lends itself to the observation of unpredictable and rare events and to 
the study of the systematic behavior of the cosmic-ray flux. Because 
of the long duration of these cosmic-ray experiments, data of high 
statistical quality are secured. These experiments are confined to 
locations relatively deep in the atmosphere and therefore measure 
principally secondary cosmic rays. Primary effects can be extrapo- 
lated from measurements of secondary effects and from the atmospheric 
cascade theory. A discussion of cascade theory and cosmic-ray phe- 
nomena is given in Rossi (1952) , Hess and Mead (1965) , and Nishimura 
(1967) . 

Discussed in Part I are the various types of detectors used to 
measure secondary particles and the pertinent information concerning 
the data coverage and procurement . 

The map in figure 1 shows the location of cosmic-ray monitoring 
stations throughout the world. The WDC-A Catalogue of Data on Sdlar- 
Terrestrial Physics lists these stations in alphabetical order by 
station name. This list includes a statement of the projects going 
on at each station, the form of the collected data, and the time 
interval of coverage. The vertical cutoff rigidities, measured in 
Bv and calculated using the quiescent magnetic field model of Finch 
and Leaton for epoch 1955.0, are also given in this list for each 
observing station. 

Data taken with some of the detection devices described below are 
reported as bihour ly counting rates. Pressure data are also included. 

A given station may submit an additional data sheet for cosmic-ray data 
corrected for pressure effects. In the case of neutron monitor stations, 
this is the general practice. Pressure, temperature, and pressure 
distribution throughout the atmosphere above the station are all of 
importance for ionization chamber and cubical meson telescope detectors. 

A niomber of methods exist to correct for these effects. Hence, the 
experimenter normally supplies meteorological data, allowing the user to 
make his own corrections. In many cases, the World Data Center A (WDC-A) 
for Upper Atmosphere Geophysics or the National Space Science Data Center 
(NSSDC) will aid the requester in obtaining more detailed information 
about instrument calibration and specified analysis procedures that can 
be used to extract the desired information from the data. 
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Figure 1 . G>smic-Ray Monitoring Stations 




NEUTRON MONITORING DEVICES 


The purpose of the neutron pile is to detect, deep within the 
atmosphere, variations in the low-energy portion of the primary 
cosmic-ray spectrum. The rigidities of particles, mainly secondaries, 
reaching these detectors vary from less than 0.01 Bv at stations near 
the geomagnetic poles to about 15 Bv at near-equatorial stations. 
Low-energy cosmic rays cause nuclear disintegrations, or "stars,” in 
the monitor. These stars are composed of charged fragments and neu- 
trons typically in the energy range 100 Mev £ E £ 1 Bev. By monitoring 
the rate of star formation, one can directly measure the low-energy 
portion of the primary spectrum. The rate of star production can be 
obtained more conveniently and reliably by detecting the resulting 
neutrons than by detecting charged fragments. 

Tabulations of intensity of the neutron flux are presented as hourly 
or bihourly values of scaled count rate versus time and can be obtained 
from WDC-A. Data corrected for barometric pressure effects are noted 
in the WDC-A Catalogue of Data on Solar -Terrestrial Physics. The two 
basic types of neutron monitors will be referred to as the "neutron 
monitor" and the "super neutron monitor." 


Neutron Monitor 


The neutron monitor uses BF3 proportional counters to detect inci- 
dent neutrons. A description of this type of neutron monitor is given 
by Simpson (1957). Data from neutron monitors are available through 
WDC-A and span the time interval from January 1957 to the present. The 
data are available in various forms: magnetic tape, punched cards, hard 
copy, or microfilm. 


Super Neutron Monitor 

The super neutron monitor evolved during the IQSY when instruments 
with a higher counting rate were required. The total calculated count- 
ing rate capacity of the super neutron monitor is about 20 times 
greater than that of the neutron monitor. Data from super neutron 
monitors are available on magnetic tape and punched cards through WDC-A. 



They span the time interval from January 1962 to the present. A sample 
of typical super neutron monitor data distributed by WDC-A in hard-copy 
form would contain daily or hourly average counts per hour for a given 
month and year at a given station. A description of this type of moni- 
tor is given by Carmichael (1968) . 

Super neutron monitor cosmic-ray data from four stations (Churchill, 
Deep River, Climax, and Dallas) are available (from November 1964 to 
the present) in the Solar -Geophysical Data bulletin for the month 2 
months prior to publication. The data are in tabular form and include 
daily (UT) average counting rates (counts/hour) . The data from these 
four stations can be used to estimate the rigidity dependence of fluc- 
tuations occurring in the primary cosmic-ray spectrum. More detailed 
descriptions of the super neutron monitor and of this program are con- 
tained in the Solar-Geophysical Data "Descriptive Text" (February 1970) . 


MESON MONITORING DEVICES 


Mesons of various types are generated by collisions between inci- 
dent primary cosmic rays and atmospheric particles . Meson data are 
presented as hourly or bihourly values of scalar count rate versus time. 
These data must be corrected for barometric and temperature effects. 
There are several instruments used to measure this component of the 
cosmic-ray beam. 


Ionization Chamber 


The ionization chamber measures the total ionization produced by 
secondary cosmic rays capable of passing through approximately 10 cm 
of lead. It measures primarily the y-meson component of the secondary 
radiation. Hence, the response of the ionization chamber detector is 
similar to that of the counter telescopes, which are discussed later, 
but it is more sensitive to heavily ionizing events, e.g., bursts and 
showers . The ionization chamber is also more sensitive to nucleons 
and air showers than the counter telescope. Ionization chamber data 
are available on punched cards (1957 to 1959) or in hard-copy form 
(1960 to the present) through WDC-A at Boulder. The typical ionization 
chamber detector data distributed by WDC-A in hard-copy form contain 
bihourly mean values, probably uncorrected for barometric pressure, as 
a function of day, month, and year for a given station. This type of 
detector is described in reference 3. 
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Cubical Meson Telescope 


The early models of the cubical meson telescope consisted of three 
trays of Geiger tubes in coincidence. A detailed description of this 
telescope is given by Elliot (1957) . Carmichael (1968) describes a 
newer version (MT-64) which uses scintillation plastic. The detector 
responds to all incident particles that can penetrate the lead absorber 
and arrive within the solid angle of the detector; i.e., the detector 
responds to the hard component of secondary radiation. (This component 
is composed mainly of y-mesons near sea level.) 

Cubical meson telescope data are available on punched cards (1957 
to 1959) or hard copy (1960 to present) through WDC-A at Boulder. A 
sample of the typical cubical meson telescope data distributed by WDC-A 
in hard-copy form would contain bihourly count rates as a function of 
day, month, and year for a given station. 


Narrow-Angle Telescope 


The narrow-angle telescope is similar to the cubical meson tele- 
scope, except that the narrow-angle telescope has a much narrower 
response pattern about the vertical. A description of this detector 
is given in reference 3. 

Narrow-angle telescope data are available through WDC-A for the 
time interval 1957 to 1964 (1957 to 1959 on punched cards and 1960 to 
1964 on hard copy) . The typical narrow-angle telescope data distributed 
by WDC-A in hard-copy form contain bihourly rates as a function of day, 
month, and year for a given station. 


Underground Telescope 

The underground telescope is similar to the cubical meson telescope 
but is placed underground. This detector responds to even higher energies 
than the cubical meson telescope at sea level due to the added absorber 
between the production layer and the detector. A description of this 
apparatus is given by Vasilov et al. (1966), and further information 
concerning underground telescopes can be found in papers by Peacock 
et al . (1968) and Barton and Stockel (1968). 

Underground telescope data are available on punched cards (1957 to 
1959) or hard copy (1960 to 1968) through WDC-A. The typical under- 
ground telescope data distributed by WDC-A in hard-copy form contain 
bihourly rates as a function of day, month, and year for a given station. 
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SHOWER APPARATUS 


A shower apparatus consists of an array of detectors, covering 
several square meters, used to measure the component of cosmic rays 
having energies about 10^^ ev. The data obtained are highly complex 
and difficult to analyze without a good knowledge of the experimental 
setup and the corrections to be made. An extensive bibliography and 
a discussion of air showers and shower apparatus appear in Fujimoto 
and Hayakawa (1967). A limited amount of data is available through 
WDC-A at Boulder for the period 1957 to 1967. 
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INTRODUCTION 


The solar atmosphere is divided into three layers. The lowest of 
these is the cool, dense, weakly ionized photosphere. The highest is 
the hot, tenuous, almost fully ionized corona, of which the interplane- 
tary plasma is the continuous extension. The 10, 000-km- thick chromo- 
sphere serves as the transitional region. See Zirin (1966) . 

The interaction of the differentially rotating solar atmosphere 
with the solar magnetic field gives rise to solar active regions, which 
show statistically an 11-year periodicity in number and location. An 
individual active region may persist for several solar rotation periods 
(27 days); within an active region, enhancements of emissions called 
flares may last from seconds to hours. See Tandberg-H^nssen (1967). 

Electromagnetic radiation emitted from the solar atmosphere con- 
sists of three components corresponding to the three solar activity 
time scales. The quiet sun component is that component observed when 
all contributions due to individual active regions have been deleted; 
coronal variations with the solar cycle are manifest in this component. 
The slowly varying component, which corresponds to a solar rotational 
period, is due to individual active regions. The burst component, 
associated with solar flares, has a time scale of seconds to hours. 

There are four major periodicals that list solar data. Three of 
these are discussed in Part IX of this Handbook : Solar-Geophysical 

Data , Geophysics and Space Data Bulletin , and the Quarterly Bulletin 
on Solar Activity . The fourth source is the weekly Preliminary Report 
cind Forecast of Solar-Geophysical Activity prepared by the Space Dis- 
turbance Forecast Center of the National Oceanic and Atmospheric Admin- 
istration (NOAA) in Boulder, Colorado. This publication summarizes so- 
lar and geophysical events of the preceding week and forecasts solar 
related activity for the coming few weeks. Due to its preliminary na- 
ture, this publication is inappropriate for referencing in more defin- 
itive publications. 

The remainder of this part is divided into sections based on the 
wavelength of the observed solar electromagnetic radiation. The wave- 
length range covered varies from radio wavelengths (kilometers to 
millimeters) to X-ray wavelengths (angstroms and less) . Solar particle 
radiation is discussed in Part III in this Handbook. 
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RADIO OBSERVATIONS 


Solar radio emissions originate from both thermal and nonthermal 
processes in the solar atmosphere (chromosphere and corona) . These 
emissions are divided into the three previously described components, 
the quiet sun and slowly varying components, both of which have pri- 
marily theiinal origins, and the burst component, which has both non- 
thermal and thermal phases . A comprehensive review of solar radio 
astronomy is given by Kundu (1965). The instrumentation used to measure 
the solar radio emissions is described in Kraus (1966). 

Observatories engaged in measuring solar radio emissions are listed 
in table 1. Because of their length, the tables referenced in this 
section on Radio Observations are included at the end of the section, 
beginning on page 23. The first line of each station listing gives 
the following information in the order and format listed: station name; 

alternate names ; the International Astronomical Union abbreviation that 
is used in the Quarterly Bulle tin on Solar Activity; other abbreviations; 
and the latitude, longitude , and observing hours in UT. The observing 
hours are given, where possible, for two time periods: the first and 
fourth quarters of the year and the second and third quarters . 

The subsequent lines list: the frequency of observation in MHz; 
the type of observation (B = bursts are reported, F = flux measurements 
are reported, V = variability index is reported, P = polarization in- 
dices are reported, SS = one-dimensional strip scans are reported, 

MP = two-dimensional maps are given, and SP = spectrograms are given); 
the years that the apparatus was in operation (a year followed by a 
blank indicates that the apparatus is believed to be currently in 
operation, a '?' indicates that operation ceased at an unknown date); 
and the sources from which the data can be obtained. Data from Eastern 
Hemisphere stations did not appear in Solar-Geophysical Data before 
1969. In the case of the Canary Islands, Carnarvon, and Houston 
stations, the data can be obtained only through the WDC-A for Upper 
Atmosphere Geophysics. The data are usually available in hard-copy 
form. 


Flux, Variability, and Polarization Measurements 

Flux densities, in units of lO"^^ watts meter"^ Hz~^, of the sum 
total of the radio emissions from the entire visible solar disk are 
measured at numerous frequencies by many observatories. These flux 
densities are measured instantaneously at solar meridian passage (local 
noon). If a burst occurs during meridian passage, the appropriate flux 
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is interpolated. The fluxes are useful as indices of the level of solar 
activity and reflect the flux level due to the quiet sun and slowly 
varying components. The fluxes can also be the mean solar flux densi- 
ties observed during a specified interval of time, such as an observing 
day or a 3-hour observing period. 

Additionally, some stations indicate the relative variability of 
the flux during the observing period. A scale of 0 to 3 is used, 0 
indicating quiet conditions and 3 indicating violent variability. 

Some stations also indicate the polarization of the observed solar 
radio flux. This polarization is described as being right- or left- 
handed (in the radioelectric sense, a clockwise rotating approaching 
wave is defined as left-handed) and having a magnitude of greater than 
50%, less than 50%, or zero. The percent of polarization is defined 
as 100 times the difference of the right-handed and left-handed flux 
divided by the sum of the two components . 

More detailed descriptions of the flux measurements are given by 
Fokker (1968) and in the descriptive text accompanying Solar- Geophysical 
Data . Calibration techniques used to determine the absolute flux values 
are described by Medd and Covington (1958) and Priese (1969) . 

Stations engaged in absolute flux measurements are listed in order 
of observing frequency on table 2. If parentheses enclose the station 
name, it is believed that the station is no longer in operation (no 
data published for 1968 or later) . The letter P or V after a station 
name indicates that polarization (P) or variability (V) indices are 
given. Figure 2 illustrates the approximate coverage obtained by 
stations making solar flux observations since the IGY. Individual 
stations, although not identified on figure 2, may be identified by 
table 2. Since flux measurements are used to provide a background 
level of solar activity, the observing hours are not critical and 
hence are not illustrated. 

Ottawa daily noon 10.7-cm flux values (1947 to the present) are 
available on the magnetic tape generated by the European Space Research 
Organization (ESRO) (Lenhart, 1968). This tape also contains geomagnetic 
and sunspot indices and is available from WDC-A, Boulder, WDC-A, Rockville, 
or NSSDC. Ottawa and Sagamore Hill data (January 1967 to the present) 
on punched cards are also available from WDC-A, Boulder. 
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EAST -WEST SOLAR SCANS 

December 1969 



10,7 cm . 21 ,cm 

Fan - Beam wif4i 1 .5 minufes of arc Fan - Beam with 2 minutes of arc 

E-W Resolution E-W Resolution 

Figures. Example of a Strip Scan 
Daily Maps of the Sun 

The distribution of sources on the solar disk is measured daily by 
several observatories. This distribution may be measured in only one 
dimension (typically requiring 4 minutes of observing time), as in a 
strip scan in which all sources lying in a line perpendicular to a point 
on the strip are summed. It may also be measured in two dimensions 
(typically requiring 1 hour of observing’ time) , in which case the data 
are usually portrayed in the form of a map or a matrix. Examples of 
strip scans and solar maps are given in figures 3 and 4, respectively. 
Table 3 lists stations making one- or two-dimensional scans of the sun, 
ordered by frequency. 


Burst Observations 


Increases in the intensity of solar radio emissions occur in con- 
junction with solar flares (Smith and Smith, 1963, and Wild et al., 
1963). These solar emissions, called bursts, events, or radio flares, 
are classified by their spectral and temporal profiles. The char- 
acter of bursts observed at meter wavelengths is markedly different 
from the character of the bursts observed at microwave wavelengths. 
Spectral information is most important in classifying solar bursts 
in the 5- to 500-MHz range, (see figure 5), while the temporal profile 
alone is sufficient to classify the higher frequency microwave bursts. 
Further, the classification scheme used for microwave bursts is not 
universal. (See Fokker, 1968, and the "Descriptive Text" of Solar- 
Geophysical Data for high frequency microwave temporal profiles.) 
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TIME (minutes) AFTER START OF FLARE 
jure 5. Spectra of Low-Frequency Bursts 




Single Frequency Burst Observations 


Stations listing solar radio bursts are given on table 4, and the 
coverage obtained since the IGY is illustrated on figure 6. The 
observing hours currently covered are shown on figures 7a and 7b. The 
radio flare listings usually include the start time, the time of peak 
intensity, the end time or duration, the peak flux density, the average 
flux density, and an indication of the overall character of the event. 
Occasionally, temporal profiles are given for microwave bursts. In 
addition to the hard-copy data indicated in table 1, WDC-A, Boulder, 
holds single frequency burst data on punched cards from the following 
observatories: Ottawa (April 1966 to the present); Sagamore Hill 

(April 1966 to the present); Pennsylvania State University (April 1966 
to the present); Boulder (April 1966 to July 1966, March 1967 to Jan- 
uary 1968, July 1968 to the present) ; Washington State University 
(July 1966 to December 1967) ; Sao Paulo (October 1967 to the present) ; 
and San Miguel (October 1967 to the present) . 


Burst Spectrograms 

Spectral observations in the 5- to 500-MHz range are made by 
observatories equipped with receivers that sweep the frequency range 
of interest. These observatories are listed on table 5, and the 
coverage obtained is illustrated on figure 8. The data are in the 
form of tabular listings, but spectrograms may be available. 

Recently, satellites have made possible radio burst observations 
at frequencies less than the lower limit imposed on ground-based 
observers by the earth's ionosphere (i.e., about 10 MHz). These 
observations are made with multichannel receivers. The only satel- 
lite data currently available are the data from the University of 
Michigan OGO 3 radio astronomy experiment, and these are available 
at the National Space Science Data Center on microfilm. They cover 
the frequency range of 2 to 4 MHz (3 to 4 MHz after October 1966) 
in the time period from June 1966 to September 1967. Data from 
the Goddard Space Flight Center experiment on the ATS 2 satellite 
covering .45 to 3.0 MHz in seven steps will be available soon at 
NSSDC. These data cover the period from April to July 1967. 
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TABLE 1 


OBSERVATORY LIST 


Station lAU 

Name Alternate (s) Abbreviation (s) Coordinates 


Observing 

Hours 


ABASTUMANI, Abastumany / ABA, ABST / 42N, 43E / 8-12, 6-12 

221 (or 209) / B, F / 57-62, 64- / WDC-A, QBSA, SGD 


Frequency (MHz) Type of Years 

Observation Observing 


Sources 


ABASTUMANI, Abastumany / ABA, ABST / 42N, 43E / 8-12, 6-12 

221 (or 209) / B, F / 57-62, 64- / WDC-A, QBSA, SGD 

ABASTUMANY, see ABASTUMANI 

ALGONQUIN, see OTTAWA 

ANN ARBOR, see McMATH-HULBERT 

ARCETRI / ARC*, ARCE / 44N, HE / 7-15** 

186 / F / 57-60 / SGD 9285 / B / 69- / SGD 

225 / F / 63-64 / SGD 1420 / B / 69- / SGD 

ATHENS, Penteli, Pentili / ATH* / 38N, 24E / 7-15 

2980 / B, F / 65-66 / G§SD, WDC-A 

BERLIN-ADLERSHOF / BER, see HEINRICH HERTZ INSTITUTE 

BIG PINE, Owens Valley / BPN* / 37N, 118W / 16-24** 

450-1000 / SP / 60-61 / SGD, WDC-A 

BJURAKAN (via Moscow IZMIRAN) / BJU / 40N, 44E / 6-9, 6-12 

191, 203, or 209 / B, FV / 57-59 / QBSA, WDC-A 

BLINDERN, see OSLO 

BOEING, see SEATTLE 

BOLOGNA / BOL / 44N, HE / 8-15, 8-15 

327 / B, F / 64- / QBSA, WDC-A 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate time 
of solar meridian passage. 
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TABLE 1 (continued) 


BONN, Stockert Rad. Obs., Eschweiler / BON* / 51N, 07E / 8-16** 

36 / F / 67 / WDC-A 

35,000 / B / 68 / WDC-A 

BORDEAUX, Floirac / BOR, BORD / 44N, OIW / 8-16** 

930 / B, F / 68- / WDC-A, SGD 

BOULDER, see HIGH ALTITUDE OBSERVATORY and NOAA 

BUENOS AIRES, see OSS. DI FISICA COSMICA 

BUKAVU, see I.R.S.A.C. 

CAMBRIDGE, see CAVENDISH 

CANARY ISLANDS (via HOUSTON) / CAN* / 28N, 15W / 8-18 

4995 / B / 67- / via WDC-A 

2695 / B / 67- / via WDC-A 

1420 / B / 67- / via WDC-A 

CARNARVON (via HOUSTON) / CAR* / 25S, 114E / 23-9 

4995 / B / 66- / via WDC-A 

2695 / B / 66- / via WDC-A 

1420 / B / 66- / via WDC-A 

CAVENDISH, Cambridge / CAV / 52N, OOE / 10-15 

178 / B, F / 57-58 / WDC-A 
81 / B, F, V / 57-58 / WDC-A 

CHALMERS INSTITUTE OF TECHNOLOGY, Onsala / CIT / 57N, HE / 6-18, 6-18 
150 / B, F, V / 57-58 / QBSA 

COLORADO, UNIVERSITY OF, see HIGH ALTITUDE OBSERVATORY 

CORNELL UNIVERSITY, Ithaca / COR / 42N, 76W / 14-21, 12-21 

430 / B / 65 / WDC-A 

201 / B, F, V / 57-58 / QBSA, SGD, WDC-A 

CRACOW, Krakow / CRA / SON, 19E / 9-12 

810 / B, F, P / 57-60 / QBSA, WDC-A 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


CRIMEAN ASTROPHYS. OBS., Simferopol, Siraeis / CRL, SIM, CRIM / 

44N, 34E / 9-12 

3100 / B, F / 59-62, 64- / QBSA, WDC-A, SGD 

1000 / B, F / 66 / QBSA, WDC-A 

208 / B, F, V / 57-62, 65-? / QBSA, WDC-A 

CULGOORA / CUL* / 30S, 147E / 22-8 
10-210 / SP / 67- / WDC-A 

8-2000 / SP / 69- / SGD, WDC-A 

DELHI, New Delhi / DEL* / 28N, 77E / 3-11** 

2000 / F / 66- / WDC-A 

DOMINION RADIO AST. OBS., see PENTICTON 

DOURBES, see UCCLE 

ESCHWEILER, see BONN 

FLEURS, see SYDNEY 

FLOIRAC, see BORDEAUX 

FORT DAVIS, see HARVARD 

GIZA / GIZ* / 30N, 31E / 6-14** 

GORKY, Zimenki / GOR, GORK / 56N, 44E / 6-12, 5-13 

19,000 / F, V / 59 / QBSA 

9100 or 9375 / B, F / 57-60, 64-65, 67- / WDC-A, QBSA, SGD 

3850 / B, F / 67- / QBSA, WDC-A, SGD 

2950 or 3000 / B, F, V / 57-59, 67- / QBSA, WDC-A, SGD 

650 / B, F, V / 57-59, 67- / QBSA, WDC-A, SGD 

200 / B, F / 57-59, 67- / QBSA, WDC-A, SGD 

100 / B, F, V / 67- / QBSA, WDC-A, SGD 

GRAFTON, Troy, Rensselaer Polytechnic Inst. / GRA / 42N, 77W / 13-21** 

18 / B / 57-62 / QBSA 

GRAHAMSTOWN, see RHODES UNIV. 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


HALEAKALA, Hawaii / HAW / 20N, 156W / 17-28, 16-29 

200 / B / 57-59 / QBSA, SGD, WDC-A 

107 / B / 64-67 / QBSA, SGD, WDC-A 

18 / B / 58-62, 64-66 / SGD, WDC-A (via ESSA) 

HARESTUA, see OSLO 

HARVARD, Fort Davis / HAR / SON, 103,W / 14-23, 13-24 

10-580 / SP / 57- / SGD, QBSAV WDC-A 

2100-3900 / SP / 60, 61 / WDC-A 

HAWAII, see HALEAKALA 

HEINRICH HERTZ INSTITUTE, Berlin-Adlershof / HHI, BER, BERL / 

52N, 13E / 9-14, 8-15 

9490 / B, F / 57- / WDC-A, QBSA, SGD 

3000 / B / 57-62, 65- / QBSA, SGD 

2000 / B / 57-64- ? / QBSA 

1470 / B, F / 57- / WDC-A, QBSA, SGD 

900 / B / 57- ? / WDC-A, QBSA 

HIGH ALTITUDE OBSERVATORY, Univ. of Colorado / 40N, 105W / 14-23, 13-24 
Boulder / HAO, BOUL 

7.6-41 / SP / 59-68 / SGD, QBSA, WDC-A 

7-80 / SP / 68- / SGD, QBSA, WDC-A 

27 / B / 57- ? / WDC-A 
18 / B / 57- / WDC-A, SGD 

HIRAISO - HIR, HIRA / 36N, 400E / 22-7, 21-8 

500 / B, F / 63- / QBSA, WDC-A, SGD 

200 / B, F, V / 57- / QBSA, WDC-A, SGD 

HOLLANDIA / HOL / 38, 141E / 21-8 

500 / B, F / 57-62 / QBSA, WDC-A 

200 / B, F / 57-62 / QBSA, WDC-A 

HOUSTON, Manned Spacecraft Center / HOU* / 30N, 95W / 14-24 

(operates NASA-SPAN, includes Canary Islands and Carnarvon) 

4995 / B / 66- / WDC-A 

2695 / B / 66- / via WDC-A 

1420 / B / 66- / via WDC-A 


♦Abbreviation assigned by NSSDC. 
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TABLE 1 (continued) 


HUANCAYO / HUN*, HUAN / 12S, 75W / 13-21** 

9400 / B / 69- / SGD 

HUMAIN, see UCCLE 

HYDERABAD / HYD* / 17N, 78E / 3-11** 

30 / B / 58 / WDC-A 

IRKUTSK, Siberian IZMIRAN / IRK, IRKU / 52N, 104E / 2-9, 3-9 

9500 or 9300 / B, F / 64- / QBSA, WDC-A, SGD 
208 / B, F, V / 57-? / QBSA, WDC-A 

I.R.S.A.C., Bukavu / IRS / / 6-15 

169 / B, F, V / 57-59 / QBSA 

ITHACA, see CORNELL UNIVERSITY 

IZMIRAN, Moscow, Krasnaya Fakhra / IZM*, IZMI / 56N, 37E / 6-14** 

(Institute of Terrestrial Magnetism, 

Ionosphere, and Radio Wave Propagation - 
coordinating agency for Abastumani, 

Bjurakan, Cracow, Gorky, Kislovodsk, 

Moscow, and Crimea) 

202 / B, F / ? / SGD 

JODRELL, Nuffield Rad. Ast. Lab, Manchester / JOD / 53N, 03W / 8-16** 
3000 / B, F, P / 57-59 / QBSA, WDC-A 

2000 / B, F / 57-59 / QBSA, WDC-A 

200 / B, F / 57-59 / QBSA, WDC-A 

80 / B, F, V / 57-59 / QBSA, WDC-A 

KIEL / KIE*, KIEL / 54N, lOE / 8-15, 6-17 

1420 / B, F / 67- / QBSA, WDC-A, SGD 

420 / B, F / 67- / QBSA, WDC-A, SGD 

240 / B, F / 67- / QBSA, WDC-A, SGD 

KIEV / KEV*, KIEV / SON, 30E / 6-14** 

204 or 210 / F / 64- / WDC-A, SGD 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


KISLOVODSK, Kislovadak (via IZMIRAN) / KIS / 44N, 44E / 9-12, 6-12 
15,000 / B, F, V / 64- / QBSA, WDC-A, SGD 

6100 / B, F, P / 64- / QBSA, WDC-A, SGD 

178 / B, F, V / 57-65 / QBSA, WDC-A 

KODAIKANAL / KOD*, KODA / ION, 77E / 3-11** 

100 / B / 66- / WDC-A, SGD 

KRAKOW, see CRACOW 

KRASNAYA PAKHRA, see IZMIRAN 

LWIRO / LWI* / 02S, 28E / 6-14** 

169 / B, F / 57-59 / WDC-A 

MACAO, Macau / MAC* / 22N, 113E / 0-8** 

MACKENZIE, UN I VERS I DADE, see SAO PAULO 

MANILA / MNL, MANI / 14N, 121E / 0-8** 

8800 / B, F / 67- / G§SD, WDC-A, SGD 

4995 / B, F / 67- / G8SD, WDC-A, SGD 

2695 / B, F / 67- / G8SD, WDC-A, SGD 

1415 / B, F / 67- / G8SD, WDC-A, SGD 

27 / B / 62, 63, 67, 68 / WDC-A 

18 / B / 62 / WDC-A 

MANNED SPACECRAFT CENTER, see HOUSTON 

MCMATH-HULBERT, Pontiac, Ann Arbor / MCM*, MCMA / 42N, 83W / 13-21** 
100-580 / SP / 58-61 / WDC-A 

18 / B / 57-62, 64- / WDC-A, SGD 

MITAKA, see TOKYO 

MOSCOW (via IZMIRAN) / MOS / 55N, 37E / 7-12, 6-12 
600 / B, F, V / 57-58 / QBSA, WDC-A 

545 / B, F, V / 59-62 / QBSA, WDC-A 

208 or 202 / B, F, V / 57- / QBSA, WDC-A 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


NAGOYA, Toyokawa / NAG, TYKW / 35N, 137E / 23-6, 22-7 

9400 / B, F / 57- / QBSA, WDC-A, SGD 
4000 / SS / 57-68 / QBSA, WDC-A 

3750 / B, F, P / 57-68 / QBSA, WDC-A, SGD 

2000 / B, F, P / 57- / QBSA, WDC-A, SGD 
1000 / B, F, P / 57-69 / QBSA, WDC-A, SGD 

NANCAY / NAN / 47N, 02E / 11-13 

408 / SS / 65- / QBSA, WDC-A 

169 / SS / 57- / QBSA, WDC-A 

NATIONAL BUREAU OF STANDARDS, see NOAA, / NBS 

NAVAL RESEARCH LABORATORY / NRL* / 38N, 77W / 13-21** 

9530 / B, F / 58-59 / SGD, WDC-A 

3200 / B, F / 58-59 / SGD, WDC-A 

NEDERHORST DEN BERG, see NETHERLANDS 


NERA-UTRECHT, see NETHERLANDS 

NETHERLANDS, Nederhorst Den Berg, NERA -Utrecht / NED, NERA / 

52N, 5E / 8-15, 6-17 

9500 / B / 59-63, 65- / QBSA, WDC-A, SGD 

2980 / B / 57- / QBSA, WDC-A, SGD 

610 / B / 65- / QBSA, WDC-A, SGD 

545 / B / 57-65 / QBSA, WDC-A 

250 / B / 57- ? / QBSA, WDC-A 

252 / SS / 65-68 / WDC-A 

200 / B, F, V, P / 57- / QBSA, WDC-A, SGD 

136 / SS / 65- / WDC-A, QBSA 

NEUSTRELITZ, Operated by Heinrich Hertz Institute / NEU* / 

53N, 13E / 7-15** 


9140 

/ 

B, F 

/ 61- / WDC-A, SGD 

2920 

/ 

B / 

57-61, 65- ? / WDC-A 

2000 

/ 

B, F 

/ 64- ? / WDC-A 

1490 

/ 

B, F 

/ 60- / WDC-A, SGD 


NEW COPERNICUS UNIVERSITY, see TORUN 
NEW DELHI, see DELHI 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


NOAA, previously ESSA, Boulder / NBS, BOUL / 40N, 105W / 14-23, 13-24 
460 or 470 / B, F / 57-58 / SGD, QBSA, WDC-A 
184 / B / 66- / SGD, WDC-A, QBSA 

167 / B, F, V / 57-60 / SGD, QBSA, WDC-A 

180 / B / 60-67 / SGD, QBSA 

NUFFIELD RADIO AST. LAB., see JODRELL 

ONDREJOV, see PRAGUE 

ONSALA, see CHALMERS INSTITUTE OF TECHNOLOGY 

OSLO, Harestua, Blindem / OSL, HARS / 60N, lOE / 9-15, 5-19 

225 or 215 / B, F, V / 57- / G§SD, QBSA, WDC-A, SGD 

200 / B, F, V / 57-61 / QBSA 

OSSERVATORIO DI FISICA COSMICA, San Miguel / OFC, SANM / 34S, 58W / 12-20*’ 
408 / B, F / 66- / SGD, QBSA, WDC-A 

OTTAWA, Algonquin / OTT, ARO, OTTA / 54N, 75W / 13-21, 11-23 

2800 / SS / 68- / WDC-A 

2800 / B, F / 57- / SGD, QBSA, WDC-A 

500 / B / 57-59 / WDC-A 

50 / B / 57-59 / WDC-A 

OWENS VALLEY, see BIG PINE 

PALO ALTO, see STANFORD 

PARAMARIBO / PAR / 05N, 55W / 11-21 

545 / B / 57-68 / QBSA, WDC-A 

200 / B, F / 57-67 / QBSA, WDC-A 

PENNSYLVANIA STATE UNIVERSITY, University Park / PSU, PENN / 

41N, 78W / 13-21, 11-23 

10,700 / B, F / 64- / SGD, QBSA, WDC-A 

2690 / B, F / 64- / SGD, QBSA, WDC-A 

960 / B, F / 64- / SGD, QBSA, WDC-A 

328 / B, F / 65- / SGD, QBSA, WDC-A 

PENTICTON, Dominion Radio Ast. Obs. / PENN, PENT / 49N, 120W / 16-24** 

2700 / B, F / 64- / QBSA, WDC-A 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


PENTILI, see ATHENS 

PIRKULI / PRK* / 40N, 48E / 5-13** 

234 / B / 59-62 / WDC-A 

PONTIAC, see MCMATH-HULBERT 

POTSDAM, Tremsdorf / AOP, POTS / 52N, 13E / 9-14, 8-16 

234 / B, F, V / 57- / QBSA, WDC-A, SGD 

111 / B, F / -58, 61— / QBSA, WDC-A, SGD 

23 / B, F / 57- / QBSA, WDC-A, SGD 

PRAGUE, Ondrejov / PRA, ONDR / 49N, 14E / 7-14, 6-15 

9400 / B, F / 61-62, 66- / WDC-A, SGD 

808 / B, F, P / 59-62, 66- / WDC-A, QBSA, SGD 

536 / B, F, V / 57- / QBSA, WDC-A, SGD 

261 / B, F / 58- / QBSA, WDC-A, SGD 

231 / B, F, V / 58-61 / WDC-A, QBSA 

PROSPECT HILL, Waltham / PRO* / 42N, 72W / 13-21** 

35,000 / MP / 67- / G8SD, WDC-A 

PULKOVO / PUL / 60N, 33E / 9-12 

9375 / B, F / 57, 58 / QBSA, WDC-A 

PULLMAN, Washington State University / PMN*, WASH / 47N, 117W / 16-24** 
486 / B / 65-68 / SGD, WDC-A 

RENSSELAER POLYTECHNIC INST., see GRAFTON 

RHODES UNIVERSITY, Grahams town / RHO / 33S, 26E / 6-14** 

125 / B / 58 / QBSA, WDC-A 

RIGA / RIG / 57N, 24E / 11-14, 12-15 

220 / B, F / 57-60, 63-66 / QBSA, WDC-A 


ROMA, see ROME 

ROME, Roma / ROM* / 41N, 12E / 7-15** 

27 / B / 63-66 / WDC-A 

18 / B / 63-65 / WDC-A 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


SACRAMENTO PEAK, Sunspot / SAC* / 32N, 105W / 15-23** 

18 / B / 57-61 / WDC-A 

SAGAMORE HILL RADIO OBS., AFCRL / SAG, SGMR / 42N, 72W / 13-21, 11-23 


Hand Iton 
15,400 / B, 

F 

/ 68- 

/ G8SD, SGD, QBSA, WDC-A 

8800 / B, F 

/ 

65- / 

G^SD, SGD, QBSA, WDC-A 

4995 / B, F 

/ 

65- / 

G§SD, SGD, QBSA, WDC-A 

2695 / B, F 

/ 

65- / 

G§SD, SGD, QBSA, WDC-A 

1415 / B, F 

/ 

65- / 

G§SD, SGD, QBSA, WDC-A 

606 / B, F 

/ 

65- / 

G8SD, SGD, QBSA, WDC-A 

245 / B, F 

/ 

69- / 

G§SD, SGD, WDC-A 

19-41 / SP 

/ 

64- / 

G8SD, SGD, WDC-A 


SAN MIGUEL, see OSS. DI FISICA COSMICA 

SAO PAULO, Umuarama Obs., Univ. Mackenzie, / SAO, SAOP / 22S, 46W / 10-19 
7000 / B, F / 67- / G8SD, SGD, QBSA, WDC-A 

SEATTLE, Boeing / SEA / 47N, 122W / 16-2 

221 / B / 62-63, 65 / QBSA, WDC-A 

222 / SS / 62-65 / SGD 

SIBERIAN IZMIRAN, see IRKUTSK 
SIMEIS, see CRIMEAN ASTROPHYS. OBS. 

SIMFEROPOL, see CRIMEAN ASTROPHYS. OBS. 

SLOUGH / SLG, SLOU / SIN, OOE / 9-15, 8-17 

71.000 / B / 68- / QBSA, WDC-A, SGD 

19.000 / B / 68- / QBSA, WDC-A, SGD 

9400 / B / 68- ? / QBSA, WDC-A 

2800 / B / 68- / QBSA, WDC-A, SGD 

STANFORD, Palo Alto / STN* / 37N, 122W / 20-21 

2700 / MP / 60- / SGD, WDC-A 

STOCKERT RADIO OBS., see BONN 

SUNSPOT, see SACRAMENTO PEAK 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


SYDNEY, Fleurs / SYD / 36S, 151E / 23-5, 22-6 

1420 / B, F / 57-61, 67-68 / QBSA 

1420 / SS / 58-62 / QBSA 

1420 / MP / 57-62, 64- / WDC-A, SCO 

690 / SS / 65- / SCO, WDC-A, QBSA 

600 / B, F / 57-58 / QBSA 

10-250 / SP / 57-68 / QBSA, WDC-A 

TOKYO, Mi taka / TOK, TOKO / 35N, 139E / -6, 23-7 

35.000 / B, F / 68 / WDC-A 

17.000 / B, F, P / 64- / QBSA, WDC-A, SGD 

9500 / B, F, P / 57-64 / QBSA, WDC-A 

3000 / B, F / 57-63 / QBSA, WDC-A 

612 / B, F / 67, 69- / WDC-A, SGD 

408 / B / 66-67 / WDC-A 

326 / B / 66-67 / WDC-A 

227 / B, F, V, P / 57-68 / QBSA, WDC-A 

200 / B, F, V, P / 57-63 / QBSA 

100 / B, F / 57-63 / QBSA, WDC-A 

67 / B, F / 57-63 / QBSA, WDC-A 

TORONTO / TTO* / 43N, 75W / 13-21** 

520 / B / 57-59 / WDC-A 

TORUN, New Copernicus Univ. / TOR / 53N, 18E / 6-14, 9-15 

127 / B, F, V / 58- / QBSA, WDC-A 

TOYOKAWA, see NAGOYA 

TREMSDORF, see POTSDAM 

TRIESTE / TRI, TRST / 45N, 13E / 10-14 

239 / B, F / 66- / G§SD, QBSA, WDC-A, SGD 

TROY, see GRAFTON 

TUBINGEN, see WEISSENAU 

UCCLE, Humain, Dourbes / UCC, UCCL / SON, 4E / 8-15, 6-17 

600 / B, F, V / 57- / QBSA, WDC-A, SGD 

169 / B, F / 57-61 / QBSA, WDC-A 

UMUARAMA, see SAO PAULO 

*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 1 (continued) 


UNIVERSITY PARK, see PENNSYLVANIA STATE UNIVERSITY 

UPPER VAN NORMAN / UVN* / 34N, II8W / 16-24** 

100,000 / MP / 67- / WDC-A 

USSURIJSK, Ussurisk, Voroshilov / USS, VORO / 43N, 132E / 23-3 

208 / B, F / 59- / QBSA, WDC-A, SGD 

VILLA ELISA / VEL* / 34S, 58W / 12-20** 

2695 / F / 67-68 / WDC-A 

10-210 / SP / 66-67 / WDC-A 

VOROSHILOV, see USSURIJSK 

WALTHAM, see PROSPECT HILL 

WASHINGTON STATE UNIVERSITY, see PULLMAN 

WEISSENAU, Tubingen / WEI / 47N, 9E / 7-16, 8-15 

1000 / B, F / 68- / WDC-A 

611 / B, F / 68- / WDC-A 

30-1000 / SP / 68- / SGD, WDC-A 

46-540 / SP / 66-69 / QBSA, WDC-A 

ZIMENKI, see GORKY 


Abbreviations for Sources : 

QBSA - Quarterly Bulletin on Solar Activity 

SGD - Solar-Geophysical Data 

WDC-A - World Data Center A Catalogue of Data on Solar -Terrestrial 
Physics 

G§SD -Geophysics and Space Data Bulletin 


*Abbreviation assigned by NSSDC. 

**Estimated by adding and subtracting 4 hours from the approximate 
time of solar meridian passage. 
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TABLE 2 


OBSERVATORIES MAKING FLUX MEASUREMENTS 


Frequency (MHz) Observatory* 


35,000 

Tokyo 

19,000 

(Gorky-V) 

17,000 

Tokyo 

15,400 

Sagamore Hill 

15,000 

Kislovodsk-P 

10,700 

Pennsylvania State University 

9,530 

(Naval Research Laboratory) 

9,500 

Tokyo, Irkutsk-P 

9,490 

Heinrich Hertz Institute 

9,400 

Prague, Nagoya-P 

9,375 

(Pulkovo) 

9,140 

Neustrelitz 

9,100 

Gorky 

8,800 

Manila, Sagamore Hill 

7,000 

Sao Paulo 

6,100 

Kislovodsk-P 

4,995 

Manila, Sagamore Hill 

3,800 

Gorky 

3,750 

Nagoya-P 

3,200 

(Naval Research Laboratory) 

3,100 

Crimea 

3,000 

Tokyo, (Jodrell-V) 

2,980 

(Athens) 

2,950 

Gorky-V 

2,800 

Ottawa 

2,700 

Penticton 

2,695 

Manila, Sagamore Hill, Villa Elisa 

2,690 

Pennsylvania State University 

2,000 

Delhi, Nagoya-P, Neustrelitz, (Jodrel 

1,490 

Heinrich Hertz Institute, Neustrelitz 

1,415 

Kiel, Manila, Sagamore Hill, (Sydney) 


*Stations in parentheses are not operational at present; P or V 
indicates polarization or variability indices given, respectively. 
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TABLE 2 (continued) 


Frequency (MHz) Observatory* 


1,000 

Nagoya-P, Weissenau, (Crimea) 

960 

Pennsylvania State University 

930 

Bordeaux 

810 

(Cracow-V) 

808 

Prague -V 

650 

Gorky-V 

612 

Tokyo 

611 

Weissenau 

606 

Sagamore Hill 

600 

Uccle-V, (Moscow, Sydney) 

545 

(Moscow-V, Paramaribo) 

536 

Prague -V 

520 

(Toronto) 

500 

Hiraiso, (Hollandia) 

467 

(NBS, now NOAA) 

430 

(Cornell) 

420 

Kiel 

408 

Osservatorio di Fisica Cosmica 

328 

Pennsylvania State University 

327 

Bologna 

261 

Prague 

240 

Kiel 

239 

Trieste 

234 

Potsdam-V, (Pirkuli) 

231 

( Prague -V) 

227 

Tokyo- V, P 

225 

Oslo-V, (Arcetri) 

223 

Seattle 

221 

Abastumani 

220 

Crimea-V, (Riga) 


*Stations in parentheses are not operational at present; P or V 
indicates polarization or variability indices given, respectively. 
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TABLE 2 (continued) 


Frequency (MHz) Observatory* 


210 Kiev 

209 Irkutsk-V, (Abastumani-V, Bjurakan) 

208 Ussurijsk 

202 Moscow-V 

201 (Comell-V) 

200 Gorky, Hiraiso-V, Netherlands -V, P, (Jodrell, Hollandia) , 

Paramaribo, Tokyo-V, P, Oslo-V 
191 (Bjurakan-V) 

186 (Arcetri) 

178 Kislovodsk-V, (Cavendish) 

169 (I .R.S.A.C.-V, Uccle, Lwiro) 

167 (NBS, now NOAA) 

150 (Chalmers Institute of Technology) 

127 Torun-V 

111 Potsdam 

100 Gorky- V, (Tokyo) 

81 (Cavendish-V) 

80 (Jodrell -V) 

67 (Tokyo) 

36 (Bonn) 


*Stations in parentheses are not operational at present; P or V 
indicates polarization or variability indices given, respectively. 
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TABLE 3 


OBSERVATORIES MAPPING THE SUN 




Frequency (MHz) 


Observatory 


Strip scan 

4,000 

Nagoya 

Strip scan 

2,800 

Ottawa 

Strip scan 

1,420 

Sydney 

Strip scan 

690 

Sydney 

Strip scan 

408 

Nancay 

Strip scan 

252 

Netherlands 

Strip scan 

222 

Seattle 

Strip scan 

169 

Nancay 

Strip scan 

136 

Netherlands 

Map 

100,000 

Upper Van Norman 

Map 

35,000 

Prospect Hill 

Map 

2,700 

Stanford 

Map 

1,420 

Sydney 
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TABLE 4 


OBSERVATORIES LISTING SINGLE FREQUENCY BURSTS 


Frequency (MHz) 

Observatory* 

71,000 

Slough 

35,000 

Bonn, Tokyo 

19,000 

Slough 

17,000 

Tokyo 

15,400 

Sagamore Hill 

15,000 

Kislovodsk 

10,700 

Pennsylvania State University 

9,530 

(Naval Research Laboratory) 

9,500 

Irkutsk, Netherlands, (Tokyo) 

9,490 

Heinrich Hertz Institute 

9,400 

Nagoya, Prague, Slough 

9,375 

(Pulkovo) 

9., 140 

Neustrelitz 

9,100 

Gorky 

8,800 

Manila, Sagamore Hill 

7,000 

Sao Paulo 

6,100 

Kislovodsk 

4,995 

Canary Islands, Carnarvon, Houston, Manila, 
Sagamore Hill 

3,850 

Gorky 

3,750 

Nagoya 

3,200 

(Naval Research Laboratory) 

3,100 

Crimea 

3,000 

Heinrich Hertz Institute, Netherlands, (Tokyo 
Jodrell) 

2,980 

(Athens) 

2,950 

Gorky 

2,920 

Neustrelitz 

2,800 

Ottawa, Slough 

2,700 

Penticton 

2,695 

Canary Islands, Carnarvon, Houston, Manila, 
Sagamore Hill 

2,690 

Pennsylvania State University 


*Stations in parentheses are not in operation. 
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TABLE 4 (continued) 


Frequency (MHz) Observatory* 


2,000 

Heinrich Hertz Institute, Nagoya, Neustrelitz ,. 
(Jodrell) 

1,490 

Heinrich Hertz Institute 

1,420 

Canary Islands, Carnarvon, Houston, Kiel, Sydney 

1,415 

Manila, Sagamore Hill 

1,000 

Nagoya, Weissenau, (Crimea) 

960 

Pennsylvania State University 

930 

Bordeaux 

900 

(Heinrich Hertz Institute) 

810 

(Cracow) 

808 

Prague 

650 

Gorky 

612 

Tokyo 

611 

Weissenau 

610 

Netherlands 

606 

Sagamore Hill 

600 

Sydney, Uccle, (Moscow) 

545 

(Netherlands, Paramaribo, Moscow) 

536 

Prague 

500 

Hiraiso, (Hollandia, Ottawa, Toronto) 

486 

Washington State University 

467 

(NBS, now NOAA) 

430 

(Cornell) 

420 

Kiel 

408 

Osservatorio di Fisica Cosmica 

328 

Pennsylvania State University 

327 

Bologna 

326 

Tokyo 

320 

(Toronto) 

260 

Prague 

250 

Netherlands 

240 

Kiel 

239 

Trieste 

234 

Potsdam, (Pirkuli) 

231 

Prague 

227 

Tokyo 


*Stations in parentheses are not in operation. 
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TABLE 4 (continued) 


Frequency (MHz) 


Observatory^ 


Oslo, (Arcetri) 

Abastumani, Seattle 
(Crimea) , (Riga) 

Crimea, Ussuri jsk, (Irkutsk) 

Mos cow 

(Cornell) 

Hiraiso, Netherlands, Tokyo, (Hollandia, Honolulu, 
Jodrell, Paramaribo, Oslo) 

(Bjurakan) 

(Arcetri) 

NOAA 

(Cavendish, Kislovodsk) 

(I.R.S.A.C., Lwiro, Uccle) 

(NBS, now NOAA) 

(Chalmers Institute of Technology 
(Torun) 

Rhodes University 

Potsdam 

(NOAA) 

(Haleakala) 

Gorky , Kodaikanal , (Tokyo) 

(Cavendish, Giza) 

(Jodrell) 

(Tokyo) 

(Macao) 

(Ottawa) 

(Hyderabad) 

HAO , (Mani 1 a , Rome ) 

Potsdam 

HAO, Manila, McMath-Hulbert , (Grafton, Haleakala, 
Rome, Sacramento Peak) 


"^Stations in parentheses are not in operation. 
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TABLE 5 


OBSERVATORIES MAKING BURST SPECTROGRAMS 


Frequency Band (MHz) Observatory* 


3,900 - 

2,100 

(Harvard) 

2,000 - 

8 

Culgoora 

1,000 - 

450 

(Big Pine) 

580 - 

10 

Harvard 

580 - 

100 

(McMath-Hulbert) 

540 - 

46 

Weissenau 

250 - 

10 

(Sydney) 

210 - 

10 

Culgoora, Villa Elisa 

80 - 

7 

HAO 

41 - 

19 

Sagamore Hill 

41 - 

7.6 

HAO 


*Stations in parentheses are not in operation. 
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SOLAR OPTICAL DATA 


The photosphere emits 99% of the total solar radiant flux with a 
black-body spectrum that peaks at 5500 A, near the center of the optical 
wavelength band. Weaker chromospheric and coronal optical emissions 
are observed either over the solar disk using a characteristic spectral 
line or over the limb in white light with the solar disk occulted. 
Extensive discussions of solar optical emissions and observational 
techniques are found in Zirin (1966). Ward (1956) gives a comprehen- 
sive discussion of solar optical phenomena and the various indices of 
solar activity. 

Solar optical data are conveniently divided into six basic types 
associated with sunspots, plages, magnetic fields, flares, hydrogen a 
emission, and the corona. A worldwide network of observatories per- 
forms one or more of these observations. The WDC-A Catalogue of Data 
on Solar-Terrestrial Physics lists these observatories and their 
locations and normal temporal coverage since 1957. Figure 9 is a map 
that locates the observatories. 


Sunspots 

Sunspots, first observed in the fourth century B.C., are optically 
dark photospheric regions of low temperature and high magnetic field. 

The occurrence and characteristics of sunspots exhibit an 11-year 
periodicity. Individual sunspots are frequent, but most spots appear 
in groups (solar active regions). Comprehensive discussions of sun- 
spots and observational techniques are given by Bray and Loughhead 
(1965) and by Zirin (1966). 

The most common index of solar activity throughout the history of 
solar observation is the sunspot number. The Wolf sunspot number was 
introduced in 1848 to give a quantitative measure of solar activity. 

The Wolf number was defined as R = lOg+f, where g is the number of spot 
groups observed and f is the number of individual spots observed. 

Despite the increase in the observing capability of modem telescopes, 
the current relative sunspot number definition assures continuity with 
the Wolf number. Relative sunspot numbers are now defined as R = K(10g+f), 
where K is an observatory- dependent factor that effects the change of 
scale required to assure continuity. 
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The Zurich relative sunspot number, R^, is the most widely used 
form of the normalized set of sunspot numbers. It is based principally 
on the observations of three sites of the Swiss Federal Observatory, 
and, when necessary, on data supplied by a number of international 
sites (Waldmeier, 1961) . 

Daily values are published in the lAU Quarterly Bulletin on Solar 
Activity and in Solar-Geophysical Data . These numbers are also available 
for the period 1932 to the present on the magnetic tape generated by the 
European Space Research Organization. This tape also contains many geo- 
magnetic indices (Lenhart, 1968). Copies of this tape are held at WDC-A 
for Upper Atmosphere Geophysics, WDC-A for Geomagnetism, and at NSSDC. 

The relative sunspot numbers, as well as the sunspot area, published in 
the Quarterly Bulletin on Solar Activity from January 1957 to December 
1967 are kept on 35-mm microfilm at NSSDC and are updated as new data 
become available. Daily values of R^ and (the American relative 
sunspot number) are held on punched cards at WDC-A, Boulder, for the 
period October 1965 to the present. 

Yearly and monthly means of Zurich relative sunspot numbers for the 
periods 1700 to 1960 and 1749 to 1960, respectively, are tabulated and 
plotted in Waldmeier (1961) . This comprehensive work also contains 
tabulations and plots of daily sunspot numbers for the period 1818 to 
1960 and epochs of minima and maxima of sunspot activity from 1610 to 
1960. The remaining data (1960 to present) can be found, with a few 
exceptions, in Solar-Geophysical Data and in the Journal of Geophys- 
ical Research (until February 1969) . 

Daily sunspot numbers, obtained by 48 stations for various years 
between 1957 and the present, are kept on hard copy at WDC-A, Boulder. 

The complete station list and more detailed information on the avail- 
ability and coverage of the data are given in the WDC-A Catalogue of 
Data on Solar-Terrestrial Physics. 


In addition, 22 observatories provide tabulated data on the posi- 
tion, area, and classifications of sunspots. Two of these sites, Boulder 
and Mt. Wilson, also prepare daily solar hemisphere maps depicting the 
position and size of sunspot groups (figure 10) . Both the tabulated 
data and the hemisphere maps, covering various years between 1957 and 
the present, are available on hard copy at WDC-A, Boulder. The WDC-A 
Catalogue of Data on Solar-Terrestrial Physics should be consulted to 
determine the coverage for each station. Some of the data from the 
22 stations are selected to be published monthly in Solar-Geophysical 
Data and in the USSR Academy of Sciences Solar Data. Starting with 
January 1969, data on the solar active regions from Solar-Geophysical 
Data have been placed on magnetic tape and are available through WDC-A, 
Boulder. For the period February 1927 to December 1947, sunspot count. 
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Figure 10. Daily Solar Hemisphere Map Depicting Position 
and Size of Sunspot Groups 
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position, and area data prepared by the U.S. Naval Observatory, Washington, 
D.C., were published in the Monthly Weather Review . From 1948 to the 
present, the "Positions, Areas, and Counts of Sunspots for (Month, Year)" 
reports were published by the U.S. Naval Observatory in the form of 
monthly circulars and can be obtained from the observatory or, through 
May 1969, from NSSDC. 

Sunspot groups that meet certain criteria (maximum area during disk 
passage > 500 millionths of solar hemisphere, a or 3 magnetic classifi- 
cation, and/or associated with a major solar flare) have been isolated, 
and data on these groups have been published for the years 1954 to 1963 
in the NASA/MSC Catalogue of Solar Activity . Some of the parameters , 
given in tabulated form, are the Mt. Wilson sunspot group identification 
number, the Greenwich. sunspot group identification number, the McMath 
plage number, the sunspot group mean latitude and Carrington longitude 
during disk passage (see Coordinate Systems section. Part VIII, in this 
Handbook) , the maximum area of the group in millionths of the solar 
hemisphere, the mean area of the group, the mean magnetic class and 
field strength, the Zurich and Mt. Wilson classifications for each day 
during disk passage, and the date on which the sunspot group was first 
observed. 

Photographs of the solar disk taken by the White-Light Photographic 
Patrol at Sacramento Peak Observatory, New Mexico, are used to obtain 
information about sunspots and sunspot groups . These data are published 
in tabulated form in the Geophysics and Space Data Bulletin and cover 
the period 1964 to the present. Among the parameters given are the 
time of observation, the Zurich classification, the total number of 
spots in each group, the coordinates of the center of mass of each 
group vmder study, the ratio of total sunspot area (umbra + penumbra) 
to total umbral area, and the quality of the frame. 

"Royal Observatory Bulletins," which are published annually in 
issues of Photoheliographic Results , give information about sunspots 
and sunspot groups as a function of solar rotation number. The bulletins 
give the daily positions and areas of sunspots, the compilations of sun- 
spot group parameters (duration, area, and position) , the daily total 
areas of sunspots and faculae by hemisphere, and the mean areas and 
latitudes of sunspots and faculae. WDC-A, Boulder, has copies of these 
publications covering the period 1951 to 1960 . 

The Solar Activity Summaries (1951-1959) and the Preliminary 
Reports of Solar Activity (1960 to the present) , prepared by the High 
Altitude Observatory, Boulder, Colorado, give summaries of solar 
activity as a function of solar rotation number. The data, generally 
depicted on a Mercator- type map, include the location and size of 
sunspots and sunspot groups along with other parameters of solar activity. 
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In some reports the spot area and count of sunspot groups are tabulated 
along with such parameters as the flare index, the solar radio noise, 
the plage area, and the coronal emission for specific solar rotation 
numbers. Copies of these reports are kept at WDC-A, Boulder, and at 
NSSDC in hard-copy form. 

Visual observations of the solar disk and drawings of the observed 
sunspot distribution are made by several different observatories ; two 
of these, Boulder and Mt. Wilson, have already been mentioned. A third 
such site is the Fraunhofer Institut in Germany, which publishes daily 
hemispheric maps of the solar disk depicting various phenomena, such 
as the location and size of sunspot groups associated with the active 
sun. Duplicates of these maps in hard-copy form for January 1957 to 
June 1961 and January 1962 to the present are available through WDC-A, 
Boulder. 

t 

The Trieste Astronomical Observatory in Italy prepares- white -light 
visual sunspot drawings; the data, however, are extracted from these 
pictures before publication and given in tabulated form in the Geophysics 
and Space Data Bulletin . The parameters listed include the time and 
quality of observation, the total number of observed sunspot groups, 
the total number of spots observed, an adjusted (K = .91) relative 
sunspot number, and the sky conditions at the time of observation. 


Plages 


The chromospheric fine structure around sunspots is marked by bright 
areas called flocculi or plages. When viewed using the K line of calcium, 
these areas are called calcium plages. Further details on plages, as well 
as details on observational techniques (spectroheliographs) , may be found 
in Zirin (1966) and in Kuiper (1953) . 

Data from daily calcium plage observations made by the McMath-Hulbert 
Observatory are available through WDC-A, Boulder, and include the date of 
observation, the viewing quality, the wavelength used for the observation, 
the plage area number, a subjective activity index value, the complete 
position coordinates, the total area of the associated faculous region, 
the total area of sunspots in the active region, and an estimate of the 
plage intensity. WDC-A, Boulder, has data on 16-mm microfilm for January 
1955 through 1967 and punched card data for 1962 to 1964 and for February 
1968 to the present. 

Monthly summaries of the McMath-Hulbert plage observations for 1962 
to 1968 and January 1969 to the present can be obtained from WDC-A, 
Boulder, on punched cards and magnetic tape, respectively. The summaries, 
containing generally the same parameters as the daily reports, are also 
available in Solar-Geophysical Data for January 1955 to the present. 
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Data related to plage regions which meet specific criteria of size, 
intensity, and flare activity from 1954 through 1963 are given in the 
NASA/MSC Catalogue of Solar Activity . This publication lists, among 
other parameters, the McMath plage number, the plage category, the mean 
latitude during disk passage, the date first observed, the number of 
days visible, the average maximum area, the intensity on a 1 to 5 
scale of the total number of flares during the disk passage of the 
plage, the Mt. Wilson sunspot group identification numbers of all spot 
groups covered by the plage, and the Mt. Wilson magnetic classification 
of the sunspot groups . 

Visual observations of calcium plage distribution are made on a 
daily basis at McMath-Hulbert Observatory and are published as hemi- 
spheric maps (see figure 11) in Solar-Geophysical Data . Alongside 
each map is a listing of the plage regions (ordered according to region 
number), a subjective scaling of the plage intensity (1 to 5), the area 
in millionths of the solar disk, and the quality of the observation. 

When sky conditions have prevented the McMath-Hulbert Observatory from 
observing the plage regions, the Catania Astrophysical Observatory at 
Catania, Italy, has supplied available drawings. Similar but more 
complete maps of solar activity (including plages) are made by the 
Fraunhofer Institut and are available through WDC-A, Boulder, for the 
periods from January 1957 to June 1961 and from January 1962 to the 
present. All of the above maps that are available from WDC-A are in 
hard-copy form. 

Photographs taken in calcium K light (3933.7 A) on a daily basis 
by the Rome Astronomical Observatory and published in the Photographic 
Journal of the Sun are kept on 35 -mm microfilm at WDC-A, Boulder. The 
pictures, cataloged according to solar rotation number, cover the period 
from January 1968 to the present. 

Plage data from several other observing sites, while more irregular 
in coverage and availability, generally contain some of the more basic 
parameters; however, the formats vary widely. For complete details 
concerning the location of these sites and the availability of such 
data at WDC-A, Boulder, the user should refer to the WDC-A Catalogue 
of Data on Solar-Terrestrial Physics. 


Solar Magnetic Fields 


The solar magnetic field is approximately dipolar at high latitudes 
and is highly irregular at middle and low latitudes. This field apparently 
gives rise to a coupling between solar differential rotation and the 
11-year cycle of solar activity; the field also governs the evolution 
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of individual active regions. Field magnitudes range from a few gauss 
in quiet regions to a few thousand gauss in some active regions. Com- 
plete discussions of the solar magnetic fields appear in Bumba (1967) 
and Godoli (1967) . 

Longitudinal (line-of-sight) components of the solar magnetic 
field are measured by raagnetographs which use the Zeeman effect. Mag- 
netographs are discussed in von Kliiber (1967) and in Zirin (1966). 

The Mt. Wilson magnetograph uses the Fe(I) line (X 5250.216) and 
scans perpendicular to the solar central meridian with a resolution of 
18 arc-seconds. A complete solar scan requires 1.5 hours. The results 
of these measurements are daily solar magnetograms, an example of which 
is shown in figure 12. These computer-plotted isogauss drawings use 
solid and dotted lines for positive and negative field polarities, 
respectively. (Magnitudes are read as shown.) 

Since September 1966, daily Mt. Wilson magnetograms have been 
published in Solar- Geophysical Data; they are now being published 2 
months after the time of observation. Large-scale (17-cm-diameter) 
copies of the daily magnetograms, as well as magnetic tapes containing 
the data from which the magnetograms are constructed, are available 
from WDC-A, Boulder. 


Flares 


A flare is a short-lived sudden increase in the intensity of radi- 
ation emitted in the neighborhood of sunspots. It is best seen in Ha, 
and usually occurs in the chromosphere. Flares are characterized by a 
rise time of the order of minutes and a decay time of the order of tens 
of minutes. The total energy expended in a typical flare is about 10^° 
ergs; the magnetic field is extraordinarily high, reaching values of 
10* to 10** gauss. Optical flares are usually accompanied by radio and 
X-ray bursts, or flares, and occasionally by high-energy particle 
emissions. The optical brightness and size of the flare are indicated 
by a two-character code called the importance. The first character, a 
number from 1 to 4, indicates the apparent area. For areas of less than 
1, an "s" is used to designate a subflare. The second character indicates 
relative brightness: b for bright, n for normal, and f for faint. The 
most recent general discussion of solar flares is found in Smith and 
Smith (1963). 
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Figure 12. Daily Solar Magnetogram 


52 


Hydrogen-a solar flare data are collected by numerous solar patrol 
stations over the globe. The February 1970 Solar-Geophysical Data 
Descriptive Text lists 53 stations. Data from these stations are found 
in the Solar-Geophysical Data bulletin and in the lAU Quarterly Bulletin 
on Solar Activity . 

The Solar-Geophysical Data bulletin contains two listings of solar 
flare reports — those published the month following the flare occurrence 
and those published 6 months after occurrence. The former list includes 
all reported flares of importance greater than one, giving the observa- 
tory, the date, the start time, the maximum and end times, the solar 
latitude and longitude, the McMath plage region, the duration in minutes, 
the importance, the observation condition and type, the area, and the 
maximum width of the Ha intensity. A 1-month listing of subflares is 
also provided indicating date, time, and location. A full discussion 
of these parameters is found in the Solar-Geophysical Data Descriptive 
Text which appears in February of each year. 

The listing published in Solar-Geophysical Data 6 months after 
occurrence consists of those flare reports filtered according to a 
program used at the Observatoire de Paris, 92 Meudon, France, to pro- 
duce the lAU Quarterly Bulletin on Solar Activity flare listing. The 
format of the 6-month listing is the same as that of the 1-month flare 
listing described above. The grouped flare data appearing in the 
6 -month listing are also available from WDC-A, Boulder, on punched 
cards and magnetic tape for the period 1955 to the present. 

Written in French, the "Eruptions Chromospheriques" section of the 
lAU Quarterly Bulletin on Solar Activity contains flare data for approxi- 
mately 14 months after flare occurrence. 


Hydrogen-Alpha Spectroheliograms 

Ha spectroheliograms are generally taken with telescopes equipped 
with a half-angstrom bandwidth Halle filter as described in chapter 2 
of Zirin (1966) . Daily photographs of solar disk images 15 cm in diam- 
eter and sequence photographs of the same type taken at approximately 
10-second intervals can be obtained from the World Data Center A for 
Upper Atmosphere Geophysics, which has taken over data formerly distrib- 
uted by NASA's Solar Particle Alert Network. 
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Solar Corona 


The corona, the outer layer of the sun's atmosphere, is a hot 
plasma (1 to 2 million degrees K) consisting mainly of fully ionized 
hydrogen and helium with smaller amounts of heavier, highly ionized 
atoms . The high temperatures are apparently caused by the dissipation 
of sound waves and hydromagnetic waves generated below the corona. 

There are three major spectroscopic components of importance in 
the corona. The first is the L or E (emission) component, which consists 
of 31 major emission lines lying between 3328 A and 10,728 A. Coexist- 
ing with this emission spectrum is a weak, continuous spectrum com- 
posed of two components: the K (electron) component and the F (dust) 

component. The E component, however, is the most widely used in 
determining coronal properties . 

The various lines that comprise the emission component are generally 
divided into three classes based on the degree of excitation. Class I 
consists of the lowest excitation lines, which are most prominent near 
sunspot minimum and in regions of low activity. Class II lines are 
intense near plages and sunspots, particularly at sunspot maximum. The 
lines of Class III are less frequently seen, appearing only in regions 
of very high ionization such as around loop prominences and directly 
above large sunspots. 

In practice, only the Class I red Fe X line (X 6374 A) and the 
Class II green Fe XIV line (X 5303 A) are observed on a regular basis 
by coronagraph stations. The distributions of the above two lines, as 
well as of the Class III yellow Ca XV line (X 5694 A) when bright enough, 
are studied throughout the entire solar cycle to obtain information 
on temperature variations in the corona. 

More detailed accounts of the solar corona and its properties can 
be found in Shklovskii (1965), Rosch (1967), and Zirin (1966). 

The coronagraph, which is used to observe the coronal emission 
lines, consists of a special telescope with an occulting disk placed 
at the focus of the objective lens to block out the direct light of the 
sun. For a comprehensive treatment of the theory and operation of the 
coronagraph refer to Zirin (1966), McMath (1953), and Evans (1953). 

The World Data Center A, Boulder, has daily maps of the sun's 
corona, supplied by the Fraunhofer Ins ti tut, for January 1957 through 
the present. In addition, WDC-A has the following coronal data availa- 
ble; numerous HAO reports that list coronal emission line intensities 


54 



as a function of solar rotation number for a variety of periods from 
1961 to 1965, the lAU Quarterly Bulletin on Solar Activity from January 
1957 to June 1966 (giving daily intensities and isophotes of coronal 
activity) and from July 1957 through December 1959 (giving coronal 
indices) , and coronal line indices in Solar-Geophysical Data for July 
1957 through December 1958 and daily charts with three intensity levels 
from January 1969 to the present. 

Spectrograph ic coronal observations from 1964 to 1969 are given in 
the Geophysics and Space Data Bulletin. 


SOLAR X RAYS 


Solar X-ray emissions, like solar radio emissions, originate from 
both thermal and nonthermal processes that take place primarily in the 
solar corona. Solar X-ray emissions can be divided into three compo- 
nents: the quiet sun component, the slowly varying component, and the 

burst component. Nonthermal X-ray emissions occur only in the hard 
X-ray emissions (photon energies greater than about 10 kev or wave- 
lengths less than about 1 A) accompanying solar flares. Solar X-ray 
astronomy has been reviewed by Mandelshtam (1965), Underwood (1968), 
and Krimigis and Wende (1970). 

Solar X rays are directly observed with satellite-bome detectors 
and indirectly observed through their effects on the ionosphere. The 
occurrences of solar X-ray bursts are indicated by the sudden iono- 
spheric disturbances (SID) discussed in Part V of this Handbook . 

Patrol types of solar X-ray measurements must be carried out using 
instruments flown on satellites. The resulting data that are available 
are listed on table 6. (All figures and tables for this section are 
included here beginning on page 61.) The first line of each experiment 
listing gives the satellite common name and the experimenter's affilia- 
tion. The second line gives the wavelength passband, in angstroms; the 
type of observations (B = burst listings, F = flux measurement listing, 
PL = flux plotted, and MAP = two-dimensional mapping); the years the 
experiment was operational (a blank after a year indicates that the 
apparatus is believed to be currently in operation) ; and the sources 
of the data. Some of the quoted passbands have recently been ques- 
tioned (Wende, 1970) . Because solar X-ray instrumentation must be 
flown on satellites, the data quite often contain periodic and sig- 
nificant time intervals when- no data were obtained due to the satel- 
lite's being eclipsed. 
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Ion chambers, proportional counters, Geiger tubes, and scintil- 
lators are used to measure the integral solar X-ray flux within 
specific wavelength or energy passbands. A review of such instru- 
mentation is given by Boyd (1965) and Neupert (1969), and some cali- 
bration techniques are discussed by Wende (1969). 


Flux Measurements 


Solar X-ray_fluxes , integrated over the entire solar disk in units 
of ergs cm ^ sec ^ over a specified passband, have been obtained from 
several satellites. These measurements are listed in table 7 ordered 
by passband, and the coverage is illustrated in figure 13. These data 
are occasionally presented in the form of plots of the X-ray flux as 
a function of time (e.g., Drake et al . , 1969). These data are listed 
in table 8, and the coverage available is listed in figure 14. The 
X-ray fluxes are useful and sensitive indicators of the level of solar 
activity. 


Daily Maps of the X-Ray Sun 


As of August 1969, X-ray spectroheliograms taken in the 9.1- to 
10.50-A passband by an OSO 5 experiment have been printed with a 2-month 
lag time in Solar- Geophysical Data (Prompt Reports) by the University 
of Leicester/University College, London. A soft X-ray telescope of 
the type described by Giacconi and Rossi (1960) is being used, and a 
period of 5 minutes is required to scan the sun (32 arc-minutes in 
diameter) with an inherent spatial resolution of 1.6 arc -minutes . An 
example of such a spectroheliogram is given in figure 15. 


Burst Observations . 


Information on X-ray burst observations is given in table 9, and 
the coverage is illustrated in figure 16. X-ray burst data are pre- 
sented in a manner similar to that of radio bursts. The data lists 
usually include the start time, the time of peak intensity, the end 
time or duration, and the maximum flux or ratio of peak flux to back- 
ground flux, as well as descriptive remarks. Lists of observing times 
are given when there are significant time gaps in the data. 
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TABLE 6 


SOLAR X-RAY PATROLS 


Spacecraft Experimenter 

EXPLORER 33, EXPLORER 35, UNIVERSITY OF IOWA 

2-12 / B, F, PL / 7/66- / NSSDC, SGD (B only), WDC-A 


Wavelength (A) Type of Dates Sources 

Observation 


EXPLORER 33, EXPLORER 35, UNIVERSITY OF IOWA 

2-12 / B, F, PL / 7/66- / NSSDC, SGD (B only), WDC-A 

INJUN 1, INJUN 3, UNIVERSITY OF IOWA 

2-12 / F* / 6/61 - 10/63 / NSSDC, WDC-A 

OGO 1, OGO 3, UNIVEI^ITY OF MINNESOTA 

0.25-1.2 / B / 9/64 - 12/67 / NSSDC, WDC-A 

OSO 3, UNIVERSITY OF MICHIGAN 

8- 12 / B, F / 3/67 - 4/68 / NSSDC, WDC-A (B only) 

OSO 5, UNIVERSITY OF LEICESTER, UNIVERSITY COLLEGE, LONDON 

9- 10 / MAP / 9/69 / SGD, WDC-A 

SOLRAD Series, NAVAL RESEARCH LABORATORY 

44-60 / B, F, PL / 64- / SGD, WDC-A 

8-20 / B, F, PL / 64- / SGD, WDC-A 

0-8 / B, F, PL / 64- / SGD, WDC-A 

0- 3 / B / 65- / SGD 

VELA 1, LOS ALAMOS SCIENTIFIC LABORATORY 
0.5-4 / B / 8/63 / SGD, WDC-A 

0.5-10 / B / 8/63 / SGD, WDC-A 


*Data coverage very poor. 

Sources : 

WDC-A - World Data Center A for Upper Atmosphere Geophysics 
NSSDC - National Space Science Data Center 

SGD - Solar-Geophysical Data 
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TABLE 7 


SOLAR X-RAY FLUX MEASUREMENTS 
(Single Points or Time Averages) 


WAVELENGTH RANGE 
(ANGSTROMS) 


SPACECRAFT, EXPERIMENT SPONSOR 


0-8 
2-12 
2-12 
8-12 
8-20 
44 - 60 


SOLRADS, Naval Research Laboratory 
Injuns 183, University of Iowa 
Explorers 33 8 35, University of Iowa 
OSO 3, University of Michigan 
SOLRADS, Naval Research Laboratory 
SOLRADS, Naval Research Laboratory 


WAVELENGTH RANGE 
(ANGSTROMS) 


TABLE 8 

PLOTS OF SOLAR X-RAY FLUXES 

SPACECRAFT, EXPERIMENT SPONSOR 


0-8 
2-12 
8-20 
44 - 60 


SOLRADS, Naval Research Laboratory 
Explorers 33 8 35, University of Iowa 
SOLRADS,. Naval Research Laboratory 
SOLRADS, Naval Research Laboratory 


TABLE 9 


SOLAR X-RAY BURST OBSERVATIONS 


WAVELENGTH RANGE 
(ANGSTROMS) 


SPACECRAFT, EXPERIMENT SPONSOR 


0 - 3 

0 - 8 
0,25 - 1.2 

0.5 - 4 

0.5 - 10 

2 - 12 
8 - 12 
8 - 20 
44 - 60 


SOLRADS, Naval Research Laboratory 
SOLRADS, Naval Research Laboratory 
OGO 183, University of Minnesota 
VELA, Los Alamos Scientific Laboratory 
VELA, Los Alamos Scientific Laboratory 
Explorers 33 8 35, University of Iowa 
OSO 3, University of Michigan 
SOLRADS, Naval Research Laboratory 
SOLRADS, Naval Research Laboratory 
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INTRODUCTION 


Significant numbers of charged particles, mainly electrons and 
protons, are accelerated at the sun during many solar flares. Some 
of these particles are subsequently detected within the interplane- 
tary medium and at the earth. These particles serve as useful 
indicators of solar, terrestrial, and interplanetary conditions. 
Often, effects from these various regions are not separable, but 
in certain events they are. (See, for example, Cline and McDonald, 
1968.) A recent review paper with extensive references is that 
of Kavanagh et al. (1970). 

In addition to their relevance to scientific studies, particles 
with energies above about 10 Mev are also of interest because suf- 
ficiently large fluxes pose radiation hazards to space systems, 
particularly manned systems . 


SATELLITE OBSERVATIONS 


Solar Proton Monitoring Experiment 

A series of Solar Proton Monitoring Experiments (SPME) , initiated 
by scientists at Johns Hopkins University/ Applied Physics Laboratory 
and at NASA Goddard Space Flight Center, is intended to provide reli- 
able flux and spectral measurements of solar protons over at least 
half a solar cycle and over a wide flux range. Each SPME consists 
of an array of solid-state detectors that measure the combined count- 
ing rate due to fluxes of protons and alpha particles in each of the 
three energy ranges: ^10 Mev/nucleon, ^30 Mev/nucleon, and ^60 

Mev/nucleon. The alpha particle component is usually very small and 
may be thought of as introducing an ambiguity of less than 10% into 
the proton count rate. Further discussions of this point and of the 
SPME in general can be found in Kohl (1968) and in the "Descriptive 
Text" of Solar-Geophysical Data (February 1970) . 



The SPME packages have been successfully flown on the highly 
eccentric orbit Explorers 34 (IMP-F, May 24, 1967, to May 3, 1969) 
and 41 (IMP-G, June 21, 1969, to present, Oct. 1970) and on the 
1100-kra-altitude polar orbiting ITOS 1 (January 23, 1970, to present, 
Oct. 1970). These packages measure omnidirectional fluxes. A similar 
package is scheduled for launch on IMP- I (highly eccentric orbit) in 
early 1971. 

From IMP data, hourly averaged fluxes are generated using up to 
22 "good" data points for the ^ 60-Mev and ^ 30-Mev channels and up to 
44 "good" data points for the ^ 10-Mev channel. These hourly averages 
are published monthly in Solar-Geophysical Data , with a 6-month time 
lag. The data appear in the forms of plots and tabular listings. 


Pioneers 6 and 7 


The counting rates of the galactic and solar cosmic-ray protons 
observed by the University of Chicago experiments (Fan et al., 1968) 
on Pioneers 6 and 7 (launched December 12, 1965, and August 11, 1966, 
respectively) have been published since April 1969, with a 1-month 
lag, in Solar-Geophysical Data . The rates are typically given once 
per day for each of three energy channels: 0.6 to 13 Mev, 13 to 175 
Mev, and > 175 Mev. A letter flag indicates whether the flux was 
rising(R), steady(S) , or falling(F) at the time when the fluxes were 
monitored. Information on geometrical factors and quiescent counting 
rates and on the effects of electron and alpha particle contamination 
is given as footnotes to the tabular data listing. 


Pioneers 8 and 9 


The counting rates of the galactic and solar cosmic-ray protons 
observed by University of New Hampshire detectors on Pioneers 8 and 
9 (launched December 13, 1967, and November 8, 1968, respectively) 
have been published since January 1970, with a 1-month lag, in 
Solar-Geophysical Data . Solid-state telescopes are used to measure 
protons above 13.9 Mev and 64 Mev on Pioneer 8 and above 13.9 Mev 
and 40 Mev on Pioneer 9. Typically, count rates are given about four 
times per day in a tabular listing. 

In addition, hourly and daily averaged count rates are available 
in hard-copy form from Space Science Center, University of New Hampshire, 
Durham, New Hampshire. The time period of available data extends from 
launch until at least as recently as 1 year prior to the current date. 
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ATS 1 


Three Li -drifted cubical silicon detectors, designed by Aerospace 
Corporation personnel, separately measure omnidirectional electron and 
proton fluxes on ATS 1. The ATS 1 satellite was launched in December 
1966 and placed into a geostationary orbit at 6.6 earth radii (geocen- 
tric) at 150° W longitude. The experiment is further described in Blake 
et al. (1968) and in the descriptive text of Solar-Geophysical Data . 

Beginning in February 1970, Solar-Geophysical Data has published a 
tabular listing of the hourly averaged flux values, with a 1-month lag, 
for each of the two proton energy channels observed (5 to 21 Mev and 21 
to 70 Mev) . 


GROUND-BASED DATA 


Solar protons reaching the polar cap regions cause enhanced ioniza- 
tion which in turn leads to enhanced radio wave absorption. Such events 
are referred to as Polar Cap Absorption (PCA) events. From a knowledge 
of the absorption characteristics, both qualitative and quantitative 
information about the incident solar protons can be gained. Techniques 
used in studying relevant ionization enhancements are described below. 


Riometer 


The riometer detector is described in Part V of this Handbook . 

Polar cap riometer absorption is usually associated with incident solar 
proton flux (Adams and Masley, 1966). There are models relating riometer 
absorption readings for daytime polar cap absorption events to the 
corresponding particle flux above a given threshold energy. For example, 
Juday and Adams (1969) give the relationship between omnidirectional 
particle flux, J(cm~^sec" ^) , and riometer absorption reading. A, (in 
decibels) as 

J(> 11 Mev) - 295 A^ (30 MHz) 

for daytime measurements with essentially zero magnetic cutoff. The 
value of this proportionality constant has recently been called into 
question for early phases (through maximum) of events (Smart and Shea, 
1970) . 
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Several riometer stations are actively sending data to WDC-A for 
Upper Atmosphere Geophysics at Boulder, Colorado. The data are avail- 
able in the form of strip charts or film copies of strip charts for 
the time periods from 1957 to present. The stations cover geographic 
latitudes from 54° N to 80° N and 66° S to 90° S and are listed in the 
WDC-A Catalogue of Data on Solar -Terrestrial Physics . 


Ionospheric Vertical Incidence Soundings 

The ionospheric vertical incidence sounding experiments and re- 
sultant f-plots are described in Part V in the section on lonosondes. 
The f-plots (©specially f-min) can be used to indicate the presence 
of solar protons as Polar Cap Absorption events (Bailey, 1964). Data 
are available from WDC-A, Boulder, for the time period from 1957 to 
the present. The stations cover geographic latitudes from 58° N to 
87° N and 66° S to 90° S. A list of the stations is given in the 
WDC-A data catalog, which also includes an indication of the complete- 
ness of time coverage from each station. 


VHP Forward Scatter 


Information on the energy spectrum and temporal variations of 
incident protons and relativistic electrons can be derived from an 
analysis of the temporal and magnetic latitude dependences of the 
intensity of scattered radio waves (Bailey, 1964) . See Part V of 
this Handbook for a more detailed explanation of the forward scatter 
technique . The WDC-A data catalog contains information about the very 
limited amount of forward scatter data available. 
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INTRODUCTION 


Large-scale current systems, mainly internal to and secondarily 
exterior to the earth, give rise to the geomagnetic field. These 
current systems possess temporal variability, which is reflected in 
geomagnetic variations. Extensive discussions of geomagnetism can be 
found in Sugiura and Heppner (1968), in Matsushita and Campbell (1967), 
and in Chapman and Bartels (1940). 


MAGNETOGRAMS , DIGITIZATIONS 


Surface observations of the geomagnetic field and of its temporal 
variations have been made by a large number of observatories. A 269- 
member observatory list, containing geographic and geomagnetic coordi- 
nates of each station active for at least part of the 1957 to present 
time period, may be found in the WDC-A Catalogue of Data on Solar-Ter- 
restrial Physics . A map of stations (figure 17) is included here to 
show the geographic extent of data coverage. By using the techniques 
discussed in McComb (1952), components of the geomagnetic field are 
measured by magnetographs and are exhibited as functions of time on 
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magnetograms; a sample magnetogram is shown in figure 18. The magnetic 
components are defined using the local coordinate system shown in figure 
19. In figure 19, £ is the total magnetic vector, ^ is its vertical 
component, H is its total horizontal component, D is its declination 
angle, and I is its inclination angle. The less often used X and Y are 
the horizontal north and east components, respectively, of H. 

In the traces for each of the three components, fields arising from 
different sources can be distinguished. For example, the horizontal 
component H(t) can be written as 

H = Ho + Sq(H) + L(H) + D(H), 

where Ho is the- main field with its secular variation, Sq(H) is the quiet 
daily solar variation in H, L is the quiet lunar variation in H (usually 
insignificant) , and D(H) is the irregular variation in H (or the distur- 
bance field H component) . These variations and the subclassifications 
are well described in Chapman (1964) . 

Using station raagnetograms , hourly averaged values of the magnetic 
components have been generated at some stations. Further, 2.5-minute 
digitized values of magnetic components are being generated from some 
station raagnetograms at the National Oceanic and Atmospheric Administra- 
tion/U.S. Coast and Geodetic Survey (NOAA)/(USCGS) as a joint under- 
taking of that organization and NASA (Cain et al., 1966). Annual means 
are also available. 

WDC-A for Geomagnetism at Rockville, Maryland, has extensive holdings 
on microfilm of normal-run magnetograms (recording speed = 20 mm/hr) , 
rapid-run raagnetograms (4 mm/min) , hourly averages , and station K indices 
described later in this part of the Handbook . WDC-A for Geomagnetism 
also holds magnetic tapes containing the 2.5-minute digitizations, hourly 
averages, and annual means. NSSDC holds copies of a subset of the WDC-A 
microfilms and copies of all of the tapes. 

Lists of the available data since 1957 held by WDC-A for Geomagnetism 
are given by station in the WDC-A Catalogue of Data on Solar-Terrestrial 
Physics. Updated availability lists on the 2.5-minute digitizations are 
generated monthly and are available from Mr. Carl Absten, Geomagnetic 
Division, U.S. Coast and Geodetic Survey, NOAA, Washington Science 
Center, Rockville, Maryland 20852 (phone 301-496-8168) . A list of data 
available from NSSDC will be provided on request. The actual data on 
either microfilm or tape can be obtained by writing Mr. William Paulishak, 
World Data Center A for Geomagnetism, Environmental Data Service, NOAA, 
Rockville, Maryland 20852 (phone 301-496-8160). 
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NOTE: 


T 


Magnetogram generated at Point Barrow, Alaska, on February 28, 1968, 
The required scaling information is provided elsewhere on the microfilm 
from which this magnetogram was taken. 


Figure 18. Sample Magnetogram 
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INDICES 


A number of activity indices have been defined to facilitate the 
study of geomagnetic activity. Local indices are determined at indi- 
vidual geomagnetic observatories from station magnetograms. Planetary 
indices are obtained by performing appropriate worldwide averages of 
local indices. Detailed descriptions of many of the indices and their 
histories are given in Bartels (1957), Lincoln (1967), and Dwarkin and 
Mescher (1967) . A functional diagram taken from Dwarkin and Mescher 
(1967), illustrating the relationships among some of the various geo- 
magnetic activity indices, is given in figure 20. 


K Index 


The K index, adopted in 1939, is intended to measure the geomagnetic 
effect of solar particle flux. The index is determined at individual 
stations for each 3-hour interval (UT) by using the magnetic trace from 
station magnetograms after subtracting the effects of the quiet daily 
solar variation (Sq) , lunar variation (L) , solar flare effects (Sqa) , 
and the aftereffects of the disturbance field (DR; R for ring current) . 

For each field component in each 3-hour interval, a range, R, in 
gammas, is defined as the difference between the maximum and minimum 
values of the adjusted magnitude of the component observed in that 
time interval . McComb (1952) contains an equation by which the 
"component" D is converted from angular measure into a specified 
number of gammas . Then K is defined in terms of the largest of the 
three R values (two R values after 1964, the Z component being specifi- 
cally excluded in this determination). The relation between R and K is 
defined by one of about 10 loosely logarithmic scales, in each of which 
K runs between 0 and 9 . The appropriate scale is determined by the 
station latitude so that all stations will report a similar K value 
for a given event. The scale appropriate to each station is characterized 
by the lower limit of R (between 300y and 2800 y) for K = 9, and is given 
for most stations in each issue of the lAGA Bulletin No. 12 series. As 
an excunple, at a mid- latitude station, if R in gammas is between 

0-5 5-10 10-20 20-40 40-70 70-120 120-200 200-330 330-500 >500, 
then K=01 2 3 4 5 6 7 8 9. 
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Figure 20. Functional Diagram of Relationships Among Geomagnetic Activity Indices 























The World Data Center A for Geomagnetism maintains the K values 
for many observatories on microfilm. The catalog of data issued by 
WDC-A contains information on the availability of K values listed by 
station since 1957. The K values themselves for the year 3 years 
prior to publication are given in the lAGA Bulletin No. 12 series . 


Kpj^ - Fredericksburg K Index 


The Fredericksburg K index obtained in Fredericksburg, Virginia, 
has been the standard of geomagnetic activity for the United States. 
Between January 1957 and February 1969, the KpR values were given in 
the Journal of Geophysical Research for the month 5 months prior to 
publication. Prior to 1957, the K index derived from the Cheltenham, 
Maryland, magnetic records was also published in the Journal of Geo- 
physical Research . KpR values are currently listed only in the "Pre- 
liminary Report and Forecast of Solar Geophysical Activity" prepared 
by the Space Disturbance Forecast Center of NOAA, Boulder, Colorado. 


Ks - Standardized K Index 


The standardized K index, Ks, is obtained from K by using a scale 
dependent upon both station and local time . The Ks index was originated 
to remove local effects from the K indices . The Ks index runs from 
0 to 9. It is expressed in thirds of an integer; e.g., the interval 
1.5 to 2.5 is divided into thirds and denoted 2 , 2^, 2^. 

The Ks index is used almost exclusively for determining the Kp 
index . 


Kp - Planetary 3-Hour Range Index 

The planetary Kp index is determined by arithmetically averaging 
Ks for 12 stations located at geomagnetic latitudes between 48° and 63°. 
As with Ks, Kp can assume 28 values as follows: 
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Note that Oo and 9© denote one-sixth of a unit interval, while the 
other values denote one-third of a unit interval . 

Values of Kp are available from numerous sources, including: 

1. A BCD magnetic tape generated by the European Space 
Research Organization (ESRO) , copies of which are 
held at WDC-A for Geomagnetism, Rockville, Md., and 
at NSSDC, for the time period from 1932 until 1 month 
prior to the most repent update. See Lenhart (1968). 
Currently, one update is issued every 3 months. 

2. Solar-Geophysical Data for the month 2 months prior to 
publication. 

3. The lUGG/IAGA Three-monthly Bulletin on Geomagnetic 
Indices and Rapid Variations for the 3 -month period 
about 8 months prior to publication. These bulletins 
may be obtained from the International Service of 
Geomagnetic Indices, Royal Netherlands Meteorological 
Institute, DeBilt, The Netherlands. 

4. lAGA Bulletin No. 12 for the 3 years prior to publication 
and lAGA Bulletin No. 18 for the period 1932 to 1961. 

5. Journal of Geophysical Research for the month 5 months 
prior to publication, but only between 1951 and February 
1969. 

6. "Abbreviated Calendar Records" (see Miscellaneous Data 
section) . 

In addition, monthly and annual mean values of Kp, for the period 1932 
to 1962, and daily mean values of Kp, for the period 1957 to 1962, are 
published in the Handbook of Geophysics and Space Environments (Shea L. 
Valley, ed., AFCRL-USAF, 1965). 


Kn, Ks, Km Indices 


To reflect the differences in geomagnetic activity occurring in 
the earth's Northern and Southern Hemispheres, planetary indices Kn 
and Ks have been derived using data obtained by Northern and Southern 
Hemisphere observatories, respectively. Km is a mean value and is 
usually very close in magnitude to Kp. More detailed discussions of 
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these indices and a listing of their numerical values for the time 
period 1964 to 1967 are given in Mayaud (1968) . The little used 
Southern Hemisphere Ks index is not to be confused with the previ- 
ously discussed standardized station index Ks; the context in which 
Ks is used should make clear which meaning is intended. 


ak, Ak Indices 

The ak and Ak indices, introduced in 1951, are determined from 
the K index in a way that they are linearly related to the geomagnetic 
disturbance amplitude. As such, they are more appropriate in per- 
forming many statistical averages than are the K indices themselves. 

The ak values are determined from K according to 

K=0 1 2 34 5 6 7 89 

ak = 0 3 7 15 27 48 60 140 240 400 

The ak index may also be considered an equivalent range of the 
most disturbed component, given in units fy where the factor f is the 
ratio of the station lower limit for K = 9, to 250. For example, if 
the range is R = lOOy while the lower limit for K = 9 is 500, then 
K = 5, and, by the conversion scale, ak = 48. The factor f is now 2, 
so that ak = 96y gives an equivalent range. 

The Ak index is called the equivalent daily amplitude and is the 
arithmetic average of the eight ak values per Greenwich day for a 
given observatory. Neither ak nor Ak values are published on a 
regular basis . 


ap, Ap Indices 

The 3-hour planetary index, 
according to the following table: 

Kp = 0 0^ 1 1 12 2^ 2^ 

ap = 0 23456 79 

Kp = 7 7 7 8 8 

ap = 111 132 154 179 207 


ap, is determined from the Kp index 


333444555666 
-0 + -0 + -0 + -0 + 


11 15 18 22 27 32 39 48 56 67 80 94 


8 9 9 

+ - o 

236 300 400 
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The arithmetic average of the eight ap values per Greenwich day is 
Ap, the daily equivalent planetary amplitude. Since ap and Ap are lin- 
early related to disturbance amplitude while Kp is logarithmically re- 
lated, the use of ap and Ap in statistical averaging is recommended over 
the use of Kp or of daily sums of Kp. 

Sources of Ap values are the same as those for the Kp indices, except 
that Ap values do not appear in the lUGG/IAGA Three-monthly Bulletin on 
Geomagnetic Indices and Rapid Variations . Values of ap are found on 
the ESRO tape and in lAGA Bulle tin No. 12 and lAGA Bulletin No. 18, with 
the same respective time coverages as for Kp values . 


C - Daily Magnetic Character Index 

The C index is a daily measure of magnetic activity assigned for each 
magnetic observatory. For a given station, C = 0, 1, or 2 according to 
whether the Greenwich day is quiet, moderately disturbed, or highly 
disturbed. Some stations have fixed range limits for C = 1. These sta- 
tions and the limits are listed in the lAGA Bulletin No. 12 series. Each 
Bulletin in the series also contains daily C values listed by station 
for the year 3 years prior to publication. 


Ci - International Daily Character Figure 

The Ci index is the arithmetic mean of the C values reported by 
individual observatories for each Greenwich day. The Ci values are 
available from 1884 to the present. 

Monthly and annual mean values of Ci are available from 1884 to 1964 
in the AFCRL Handbook of Geophysics and Space Environments . Daily val- 
ues for Ci, also given in that document, are cited for the period 1957 
to 1963. Daily values are also available in the lUGG/IAGA Three-monthly 
Bulletin on Geomagnetic Indices and Rapid Variations for the 3-month 
period about 8 months prior to publication, in Solar-Geophysical Data 
for the month 2 months prior to publication, in the Journal of Geophysical 
Research (prior to February 1969) for the month 5 months prior to publi- 
cation, and in the lAGA Bulletin No. 12 series for the year 3 years 
prior to publication. The lAGA Bulletin No. 12 series also lists the 
monthly and annual mean values since 1900, with a cutoff date of 3 years 
prior to publication. 
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Cp - Daily Planetary Character Figure 

The Cp index is a measure of magnetic activity that is determined 
for each Greenwich day from the sum of the eight ap values for that day 
according to: 


8Ap 

= Zap 

< 22 

34 

44 

55 

66 

78 

90 

104 

120 

139 

164 

190 

Cp 

= 

0.0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

Zap 

< 228 

273 

320 

379 

453 

561 

729 

1101 

1399 

1699 

1999 

2399 

3200 

Cp 

= 1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 

2.0 

2.1 

2.2 

2.4 

2.5 


Although daily values of Ci and Cp rarely differ by more than 0.2, 
Cp is more reliable (Dwarkin and Mescher, 1967). Sources of Cp values 
are the same as those for the Kp indices, except that Cp values do not 
appear in lUGG/lAGA Three-monthly Bulletin on Geomagnetic Indices and 
Rapid Variations or in the "Abbreviated Calendar Records . " 


Cq Index 

The Cg index expresses the daily geomagnetic activity by a single 
digit. The values for Ci and Cp are contracted to form the Cg values 
as follows: 



0.0 

0.2 

0.4 

0.6 

0.8 

1.0 

1.2 

1.5 

2.0 

Ci or Cp = 

to 

to 

to 

to 

to 

to 

to 

to 1 .9 

to 

0.1 

0.3 

0.5 

0.7 

0.9 

1.1 

1.4 

1.8 

2.5 

Cg 

0 

1 

2 

3 

4 

5 

6 

7 8 

9 


Values of Cg appear occasionally in Solar-Geophysical Data (e.g., 
the period 1965 to 1968 is covered in Solar-Geophysical Data , 294, Feb. 
1969) . Cg values are also found on the ESRO tape (see discussion in 
Kp section) for the interval 1932 until 1 month prior to the most 
recent update. 


Dst - Axisymmetric Component of Geomagnetic Disturbance 


A Dst(X; t) field may be defined for any field component, X, as the 
average of the D(X; t) field values as observed by stations of similar 
geomagnetic latitude but varying geomagnetic longitude. In practice, only 
Dst for the horizontal field, Dst(H; t) has been considered. At a given 
station Ds(H; t) = D(H; t) - Dst(H; t) . 
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The Dst index, the magnitude of the Dst(H) field, is intended pri- 
marily to indicate the effect of magnetospheric ring currents, especially 
during magnetic storms. Thus, observatories used for the determination 
of Dst are distributed as uniformly as possible in longitude and lie at 
low, but not equatorial, latitude. The low latitudes minimize auroral 
zone effects. Equatorial latitudes are avoided to minimize equatorial 
electro jet effects. 

A more detailed discussion of the Dst(H) index and its derivation 
are presented in Sugiura (1964) . This reference also contains a list of 
hourly Dst(H) values for the IGY (July 1957 to December 1958) and a list 
of the eight stations from which these Dst(H) values were determined. 

Hourly Dst(H) values have been generated using data from Hermanus, Honolulu, 
and San Juan; these values are found in Sugiura and Hendricks (1967) for 
the years 1961 through 1963, in Sugiura and Cain (1969) for the years 1964 
through 1967, and in Sugiura and Cain (1970) for 1968.* No Dst values have 
yet been generated for the 1959/1960 time interval. Copies of these doc- 
uments can be obtained from the authors or through NSSDC. The Dst values 
are also available on a BCD magnetic tape, copies of which are kept at 
WDC-A, Rockville, and at NSSDC. 


U, u, Ui Measures 

The U, u, and ui measures are the older indices related to ring cur- 
rent effects. Their use has largely been superseded by the use of the 
Dst index. At a given station, the U measure is the absolute magnitude 

of the difference of the daily mean values of the H component for two 

successive days as assigned to the second day. These U values are geo- 
metrically reduced to equivalent equatorial u values. The ui measure is 
related to the u measure in a nonlinear fashion such that the larger val- 
ues of u are deemphasized; ui correlates well with sunspot numbers. 

Values for ui related to u are as follows: 

u = 0.3 0.5 0.7 0.9 1.2 1.8 2.1 2.7 3.6 

ui = 0 20 40 57 79 108 108 132 140 

(u in units of lOy) 


*Post deadline note: Dst through 7/70 now available through NSSDC. 
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Published values of the u measure are based on averages of values 
of one to six near-equator equatorial stations. Monthly and annual 
mean values for U are given in the AFCRL Handbook of Geophysics and 
Space Environments for the years 1872 to 1949. In Terrestrial Mag- 
netism and Atmospheric Electricity , 37, 1-52, 1932, the following 
values are listed: 

U - monthly and annual means for the years 1872 to 1930. 

u - annual means for the years 1835 to 1871. 

Uj - monthly and annual means for the years 1872 to 1930. 


AE - Auroral Electrojet Activity Index 


The AE index is a direct measure of the amplitude of the auroral 
electrojet and is obtained as follows. A AH(t) is determined at selected 
auroral zone stations as the difference between the observed H(t) and 
the quiet time H(t). The AH(t) values for all observatories are plotted 
versus time on a single graph. The AU and AL curves are then drawn 
versus time through the maximum and minimum AH values, respectively. 

At each point in time, AE is defined as AE = AU - AL. Another index, 

Ao, is defined by (AU + AL)/2, and is the displacement of the midpoint 
between AU and AL from the reference level. 

The significance of AE and a more detailed discussion of the method 
of derivation are given by Davis and Sugiura (1966). 

Hourly values of AE have been computed from magnetic data gathered 
at 11 stations (10 for 1957) fairly well distributed in geomagnetic longi- 
tude and located mainly in the Northern Hemisphere. These computations 
have been done for the period from July 1957 to December 1964 by T. N. 
Davis et al. at the University of Alaska. The results have been pub- 
lished in the University of Alaska Geophysical Institute Report (UAG-R) 

No. 192, 194, 195, 196, 197, 198, 199, and 200. The reports contain 
plots of the AE index and tabulations of AE, AU, AL, and Ao. 

WDC-A for Geomagnetism currently has available one IBM 7094 binary 
tape per calendar year containing 2. 5 -minute AE values for the period 
from September 1964 through December 1968 and one BCD tape containing 
hourly AE values for the period from July 1957 through December 1968. 

The 2.5-minute AE values are computed at NSSDC (in cooperation with D. 
Fairfield) and are based on available Northern Hemisphere 2.5-minute 
digitized magnetogram data. The hourly AE values are based on hourly 
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averages of the 2.5-minute digitized data. At present, 11 stations 
fairly well distributed in geomagnetic longitude are represented in 
the AE values for 1964 and 1965. For 1966 and 1967, data from six 
stations, four of which lie in the longitude range 256 to 347 degrees 
geomagnetic, have been used. For 1968, four stations in the preceding 
longitude range and one outside the range have been used. As 
additional 2.5-minute digitizations become available, these AE values 
will be updated and extended forward in time. 

In addition to these AE values on magnetic tape, WDC-A also has 
plots and tabulations of these AE values on microfilm and AU and AL 
values on both magnetic tape and microfilm. NSSDC also has copies of 
these data. 


Q Index 


The Q index is a measure of geomagnetic activity assigned by high 
latitude (geomagnetic latitude > 58“ ) magnetic observatories for each 
15-minute interval centered around 00, 15, 30, and 45 minutes past the 
hour. The index is designed to study auroral and ionospheric phenomena 
over a time scale smaller than that possible with Kp. The Q scale is 
loosely logarithmic, with possible values from 0 to 11. 

The index is defined as the sum of the maximum positive deviation 
from the normal nondisturbance curve and the absolute value of the 
maximum negative deviation. The more disturbed of the two horizontal 
elements is chosen for this determination. As an example, if the extreme 
deviations from normal are +200y and +50y, the K index would be based 
on a range of 150y, whereas Q is based on the 200y maximum deviation 
from normal. (The negative deviation is zero in this case.) The uniform 
conversion scale for all stations is: 

S < 10 20 40 80 140 240 400 660 1000 1500 2000 ~ (y) 

Q=0123 4 5 67 8 9 10 11 

(S = sum of absolute values of positive and negative deviations.) 

The Q indices are available from WDC-A for Geomagnetism. The 
availability information since 1957 is listed by station in the WDC-A 
Catalogue of Data on Solar-Terrestrial Physics. 
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SELECTED DAYS 


For many studies, a selection of magnetically quiet or disturbed 
days is necessary. For each month, the five most disturbed days, the 
five least disturbed days, and the 10 least disturbed days are selected 
by lAGA personnel using criteria described in Lincoln (1967) . 

Sources of selected days are the same as for Kp indices, except 
that selected days are not listed on the ESRO tape or in lAGA Bulletin 
NO. 18. 


EVENTS 


In addition to being used for the generation of indices, magneto- 
gram records are used to identify individual events. These events 
include magnetic storms, sudden impulses, polar magnetic substorms, 
micropulsations , and solar flare effects . 

The reader may wish to refer to the previous section on magneto- 
grams for availability information. Some geomagnetic activity infor- 
mation (especially on storms) is distributed on a real-time basis as 
part of the International Ursigram and World Days Service (lUWDS) 
effort (see discussion in Miscellaneous Information, Part VIII). 


Geomagnetic Storm 

A geomagnetic storm is a worldwide transient variation in the geo- 
magnetic field. In mid and low latitudes, four distinct phases in a 
typical storm can be distinguished. The storm sudden commencement 
(sc or ssc) phase is the nearly simultaneous worldwide increase in H, 
typically by a few tens of gammas. This state of ei^anced H may per- 
sist for a few tens of minutes and is called the initial phase. During 
the following main phase, H is depressed below its prestorm value typi- 
cally by several tens to a few hundred gammas for many hours. The 
recovery phase, or return to the prestorm state, usually requires a few 
days. 
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The sudden commencement and initial phases are manifestations of 
magnetospheric compression due to enhanced solar wind pressure, while 
the main and recovery phases are due to the great enhancement and 
gradual decay of the diamagnetic ring current found within the magneto- 
sphere. Further discussions of geomagnetic storms may be found in 
Sugiura and Heppner (1968) . 

Since 1961, the lAGA Bulletin No. 12 series has published for the 
year 4 years prior to publication a storm sudden commencement list 
using an "occurrence format." That is, for a given event, characterized 
by a given start time, reporting stations are separated into six classes 
according to whether the event in that station's magnetogram was: (1) 

very distinct; (2) fair, ordinary, but unmistakable; (3) doubtful; 

(4) not recorded despite satisfactory magnetogram records; (5) not 
discerned due to heavy disturbance; (6) missing because the magnetogram 
itself was missing. No event morphology information is given in these 
lists. 

The Three-monthly Bulletin on Geomagnetic Indices and Rapid Varia- 
tions of lUGG/IAGA lists the same information for the 3-month period 
about 8 months prior to publication. 

Until January 1967, the Journal of Geophysical Research published 
a monthly listing of storms, grouped by station, that had occurred 
about 6 months earlier. The begin and end times of storms, the types 
of sudden commencement, the magnetic amplitude of each field component 
during initial and main phase, and the maximiun associated K index were 
all presented. Since February 1967, this information, grouped by 
event, has appeared in ESSA's Solar -Geophysical Data for the month 2 
months prior to publication. Beginning April 1969, storm magnetograms 
for a few events for selected observatories are presented in Solar- 
Geophysical Data . 


Sudden Impulse 


Transient magnetic disturbances are similar to storm sudden com- 
mencements, but they tend to be of smaller amplitude and of less 
abruptness than the ssc's. These sudden impulses (s.i.) appear to 
result from sudden magnetospheric compressions and expansions associated 
with shock waves and other discontinuities in the interplanetary medivim 
(Burlaga, 1970), but they are not followed by main phase storms. In- 
complete lists of sudden impulses appear in the lUGG/IAGA Three-monthly 
Bulletin on Geomagnetic Indices and Rapid Variations , for the 3-month 
interval about 8 months prior to publication, and in lAGA Bulletin No. 

12 for the year 4 years prior to publication. The "occurrence format" 
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described in the preceding Geomagnetic Storm subsection is used in both 
sources. Three-month lists have also appeared in Solar -Geophysical 
Data, with varying time lags, since January 1966. 


Polar Magnetic Substorm 


A polar magnetic substorm (formerly called a magnetic bay) is a 
predominantly night-side variation, often superposed on a worldwide 
geomagnetic storm. The magnetic deviation, which may be positive 
(usually pre-midnight) or negative (usually post-midnight), can have 
an amplitude of hundreds of gammas at auroral latitudes and a signif- 
icantly reduced amplitude away from the auroral zones . Each event 
usually lasts between 1 and 6 hours. 

Recent reviews of polar magnetic substorms and of their relations 
to auroras and to magnetospheric dynamics are given in Hultqvist (1969) 
and in Feldstein (1969). 

Lists of magnetic bays are given in the lUGG/IAGA Three -monthly 
Bulletin on Geomagnetic Indices and Rapid Variations , for the 3 -month 
interval about 8 months prior to publication, and in lAGA Bulletin No. 
12 for the year 4 years prior to publication. The "occurrence format" 
described in the Geomagnetic Storm subsection is used in both sources. 
The type of bay (calm or stormy time, sharp or gradual beginning, 
associated pulsations or not) is also indicated. 


Micropulsations 

Small amplitudej quasi -periodic magnetic fluctuations are found 
occasionally in both ground and satellite magnetometer records. Occur- 
rence frequencies and amplitude depend on geomagnetic latitude and 
longitude, radial distance, and the fluctuation frequency. Pulsations 
lasting for considerable lengths of time are classed as Pci (Pc - 
continuous pulsation), Pc2, Pc3, Pc4, or Pc5 according to whether the 
fluctuation periods lie within the respective ranges: 0.2 to 5 seconds, 

5 to 10 seconds, 10 to 45 seconds, 45 to 150 seconds, 150 to 600 seconds. 
Typically, the smaller period waves have amplitudes of a fraction of a 
gamma, while the Pc5 wave may have amplitudes of hundreds of gammas, 
especially in the auroral zones. One class of Pci waves is referred 
to as pearls, or hydromagnetic emissions. 
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Review papers on the subject include those by Troitskaya (1964) 
and Saito (1969) . The relationship between magnetic pulsations and 
polar magnetic substorms is discussed in Gendrin (1970) . A recent 
paper by Liemohn (1969) , which contains a discussion of magnetospheric 
ELF signals (periodicities 0.002 to 5 sec), presents a 148-member list 
of stations making micropulsation observations. Station location, 
frequency range observed, and the name of an analyst are given. 

Lists of stations observing Pc4 and Pc5 micropulsations in their 
normal magnetograms and lists of stations observing micropulsations of 
all Pc classes in their rapid-run magnetograms are given in the 
lUGG/IAGA Three-monthly Bulletin on Geomagnetic Indices and Rapid 
Variations for the 3-month interval about 8 months prior to publication. 
Start and stop times for each occurrence are given. 

In addition to the continuous micropulsations, irregular pulsations 
are also observed. Trains of pulsations, Pt, with periods of 40 to 100 
seconds are closely associated with magnetic bays. Pt occurrence data 
are published in conjunction with magnetic bay data. 

Other irregular pulsations are classed as Pil and Pi2 according to 
whether the periods lie within the ranges 1 to 40 seconds and 40 to 150 
seconds, respectively. Pulsations of unusually large amplitude are 
denoted as Pg (giant pulsations) . Lists of stations observing Pi2 and 
Pg pulsations in their normal magnetograms and lists of stations ob- 
serving Pil, Pi2, and Pg pulsations in their rapid-run magnetograms 
are given in the lUGG/IAGA Three-monthly Bulletin on Geomagnetic Indices 
and Rapid Variations . Again"^ start and stop times for each occurrence 
are given . 

A survey of magnetospheric micropulsations observed by OGO 3 and 
OGO 5 can be found in Heppner et al. (1969). 


Solar Flare Effects 


Enhanced lower -ionospheric ionization resulting from solar flare 
X rays causes an enhancement of the dayside Sq current, which in turn 
results in geomagnetic variations. The occurrence of this effect is 
listed by observing station in the monthly letter issued by the Inter- 
national Service of Geomagnetic Indices, Royal Netherlands Meteorolog- 
ical Institute, DeBilt, The Netherlands. Quarter-yearly lists, with 
varying time lags, have appeared in Solar -Geophysical Data since 
January 1966. 
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TELLURIGRAMS 


Temporal changes in the geomagnetic field induce electric fields 
and currents in the earth's crust. For changes (pulsations) with 
periods between tenths of a second to tens of minutes, the induced 
telluric currents are measured directly at a number of stations, 
largely (75%) in the Soviet Union. A more detailed discussion of the 
concepts involved is given in Chapman and Bartels (1940) and in 
Troitskaya (1957); the latter reference also contains a description 
of experimental techniques . 

Each station generates a tellurigram daily; each tellurigram gives 
magnitude of current versus time. World Data Center A for Geomagnetism, 
Rockville, has a number of these tellurigrams on microfilm. The station 
list and the availability of the tellurigrams, given by station and by 
year since 1957, are listed in the WDC-A Catalogue of Data on Solar- 
Terrestrial Physics (June 1970) . 


GEOMAGNETIC FIELD MODELS 


For many reasons it is useful to have an analytical representation 
of the steady component of the geomagnetic field. In regions of no 
magnetic sources (i.e., no currents), a static magnetic field, may 
be represented by 

B = -VV 

where the scalar potential V is given by 
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Here, a is the earth's effective radius, Pj (cos 0) is a Schmidt normal- 
ized Legendre pol:^omial, and g™, h®, and B® are coefficients whose 
specification defines a field model . Different coordinate systems for 
0 and (|) have been used in differing models. 

There have been two basically different types of models generated 
in the past. The first type of model is based on data gathered at or 
near the earth's surface; in this type, sources exterior to the earth 
are neglected, and = ^ = 0 (all n and m) . The second type of field 
model takes into account external sources in the form of model current 
systems such as those associated with the magnetopause and the geo- 
magnetic tail. The earth's core field in this type model is taken to 
be dipolar so that gj ^ 0, while g® = h® = 0 (all other n and m) . 


Models Based on Near-Earth Data 


Models based on near-earth data are defined by the specification 
of a finite number of the coefficients g® and h^. The early history 
of these models is given in Chapman and Bartels (1940) and Vestine 
(1960) . References to more- recent representations (some of which 
present time derivatives of the coefficients) are given in Cain and 
Cain (1968) and Zmuda (1970). 

Within the framework of the international effort named the World 
Magnetic Survey, as outlined by Heppner (1963) and Zmuda (1970), input 
data for recent models have been collected by surface observatories 
(annual means of magnetogram records), ships, aircraft, and low-alti- 
tude satellites (e.g., OGO 2, OGO 4, 1964-83C, and Cosmos 49). 

As is detailed in Cain and Cain (1968), lAGA (1969), and Zmuda 
(1970) , some recent models have been appropriately weighted to yield 
the International Geomagnetic Reference Field [IGRF 1965.0] for epoch 
1965. The IGRF specifies g]^ and h® through n = m = 8 and was developed 
to provide a standard field model where the long duration of validity 
(1960 to 1969) outweighs the advantages of high accuracy. The 80 
non- zero IGRF (1965.0) coefficients and their first time derivatives 
are listed for geographic coordinates in lAGA (1969), Cain and Cain 
(1968), Zmuda (1970), and Fabicino and Peddie (1969) and are also 
shown in table 10. Also contained in Fabiano and Peddie are tables 
of the total magnetic field intensity and its annual change at every 
2® grid intersection of geographic latitude and longitude for the 
earth's surface. Mead (1970) has transformed the coefficients and 
their time derivatives to values applicable in geomagnetic dipolar 
coordinates (see the section on Coordinate Systems in Part VIII of 
this Handbook) . 
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TABLE 10 


Spherical Harmonic Coefficients — IGRF — 1965 


Main Field Secular Change 
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n 
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Y/yr 
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0.0 

2 

1 

2994 

-2006 

0.3 

- 11.8 

2 

2 

1567 

130 

- 1.6 

- 16.7 

3 

0 

1297 

0 

0.2 

0.0 

3 

1 

-2036 

-403 

- 10 . 8 . 

4.2 

3 

2 

1289 

242 

0.7 

0.7 

3 

3 

843 

-176 

- 3.8 

- 7.7 

4 

0 

958 

0 

- 0.7 

0.0 

4 

1 

805 ' 

149 

0.2 

- 0.1 

4 

2 

492 

-280 

- 3.0 

1.6 

4 

3 ' 

-392 

' 8 

- 0.1 

2.9 

4 

4 

256 

-265 

- 2.1 

- 4.2 

5 • 

0 

, -223 

0 

1.9 

0.0 

5 

1 

357 

16 

1.1 

2.3 

5 

2 

246 

125 

2.9 

■ 1.7 

5 

3 

-26 

-123 

0.6 

- 2.4 

5 

4 ' 

-161 

-107 

0.0 

0.8 

5 

5 

-51 

77 

1.3 

- 0.3 

6 

0 

47 

0 

- 0.1 

0.0 

6 

1 

60 

-14 

- 0.3 

- 0.9 

6 

2 

4 

106 

1.1 

- 0.4 

6 

3 

-229 

68 

1.9 

2.0 

6 

4 

3 

-32 

- 0.4 

- 1.1 

6 

5 

-4 

-10 

- 0.4 

0.1 

6 

6 

-112 

-13 

- 0.2 

0.9 

7 

0 

71 

0 

- 0.5 

0.0 

7 

1 

-54 

, -57 

- 0.3 

- 1.1 

7 

2 

0 

-27 

- 0.7 

0.3 

7 

3 

12 

-8 

- 0.5 

0.4 

7 

4 

-25 • 

9 

0.3 

0.2 

7 

5 

-9 

23 

0.0 

0.4 

7 

6 

13 

-19 

- 0.2 

0.2 

7 

7 

-2 

-17 

- 0.6 

0.3 
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TABLE 10 (continued) 


n m 


8 

8 

8 

8 

8 

8 

8 

8 

8 


0 

1 

2 

3 

4 

5 

6 

7 

8 


Main 

Field 

Secular 

Change 

ID 

n 

‘m 

^n 

n 

Y 

Y 

Y/yr 

Y/yr 

10 

0 

0.1 

0.0 

9 

3 

0.4 

0.1 

-3 

-13 

0.6 

-0.2 

-12 

5 

0.0 

-0.3 

-4 

-17 

0.0 

-0.2 

7 

4 

-0.1 

-0.3 

-5 

22 

0.3 

-0.4 

12 

-3 

-0.3 

-0.3 

6 

-16 

-0.5 

-0.3 
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Th.e POGO (10/68) and POGO (8/69) models generated by J. C. Cain 
from OGO 2 and OGO 4 data provide greater accuracy for the post-1965 
period than does the IGRF (1965.0) field (Cain and Sweeney, 1970). Of 
course all models become progressively less accurate as one moves in 
time away from the model epoch date. Future models will be generated 
as the need arises. It is suggested that the potential user of any 
model field consider his accuracy requirements in determining which 
of the available field models to use. 

IBM card decks containing the GSFC (12/66), POGO (10/68), POGO 
(8/69), and IGRF (1965.0) model coefficients for use with geodetic 
coordinates are available at NSSDC. A deck containing IGRF (1965.0) 
coefficients for use with geomagnetic dipolar coordinates is also 
available at NSSDC. 

The standard FORTRAN program used to generate the actual magnetic 
field at a fixed point in space from any of the field models is avail- 
able on cards at NSSDC. Also available is a fast, comprehensive 
FORTRAN routine (ALLMAG) that condenses six selected, time dependent 
models into one operational assembly, permitting a multiple successive 
selection of models and/or time periods during a single program exe- 
cution. This routine was developed by Stassinopoulos et al. (1970) 
and contains, in addition to four models mentioned above, the GSFC 
(9/65) model and that of Leaton et al . (1965). 

The field computation programs are compatible with the programs, 
also available at NSSDC, used to evaluate the B, L coordinates of a 
spatial point (refer to section on Coordinate Systems) and to trace 
out field lines. 

WDC-A, Rockville, has card decks for the 80 IGRF (1965.0) model 
coefficients and for the 168 World Magnetic Chart coefficients for 
1970. World Magnetic Chart coefficients are computed from extensive 
ground-based data by personnel at the United States Coast and Geodetic 
Survey (Rockville, Maryland) under the direction of E. B. Fabiano, and 
are used in the generation of the World Magnetic Charts issued every 
5 years by the Naval Oceanographic Office, Suitland, Maryland. WDC-A 
also has the program decks required to compute field values from either 
of these sets of coefficients. 


Models with External Sources 


As mentioned, there are field models that take into account current 
systems exterior to the earth. These models involve a specification 
of the coefficients a”* and b"" and are valid only at distances interior 
to the location of the current systems . Early versions of such models 
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assxime that the geomagnetic dipole axis is perpendicular to the inci- 
dent solar wind flow. Two early models that include the effects of 
currents located on the magnetospheric boundary are given in Mead (1964) 
and Midgley (1964). 

One of these (Mead, 1964) postulates a model magnetopause current 
system and computes nimerical values of the resultant magnetic field 
at 273 points within the magnetosphere. Least-squares fitting to the 
computed field values is used to determine the coefficients and B®. 

A parameter, r^,, which depends on the solar wind pressure and which 
scales the magnetosphere, is the subsolar geocentric magnetopause dis- 
tance. Mead presents two sets of numbers for the products A^ 

(B = 0 by symmetry, all n and m) for n ^ 6 and for n £ 2, respectively. 
Use of the coefficients with values of n up to 2 (only A° and A^ are 
non-zero) yields a fit with only 3% inaccuracy for geocentric distances 
less than 0.7 r^^. 

Using different techniques but the same general approach as Mead, 
Midgley (1964) also presents a model with coefficients having values of 
n up to 6. Midgley' s results are not significantly different from those 
of Mead. 

An early model in which a magnetic field due to a model geomagnetic 
tail current is superposed on the dipolar core field and magnetopause- 
current field is that of Williams and Mead (1965) . The postulated tail 
current lies in a semi-infinite sheet in the earth's equatorial plane, 
at distances between 10 and 40 earth radii at local midnight, and flows 
everywhere perpendicular to the midnight meridian. The field adjacent 
to the current sheet is 40 gammas. This model was generated to explain 
satellite (1963-38C) particle observations near the noon-midnight merid- 
ian plane. No Gaussian coefficients A^ and B® are presented. 

A more recent field model that includes the effect of magnetopause 
and magnetotail currents is that of Olson (1970) . This model is dis- 
tinguished from earlier models in that the geomagnetic dipole axis is 
tilted with respect to the incident solar wind; the angle between the 
geomagnetic dipole axis and the earth-sun vector is 90° - y. This model 
also incorporates a cylindrical geomagnetic tail in which current flows 
across the neutral sheet and returns across the two halves of the 
cylindrical boundary. The neutral sheet typically lies outside the 
solar magnetospheric equatorial plane and contains a current density 
that is a function of geocentric distance. 

To facilitate the computation of the magnetic fields resulting 
from the postulated current systems, the effects of magnetopause and 
magnetotail currents are considered separately. As such, two scalar 
potentials are involved and two sets of Gaussian coefficients are 
specified. Due to the linearity of all relevant equations, linear 
superposition of the fields (addition of coefficients for given n and 
m) is possible. 
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As in Mead's original model, there is a parameter, r^ef, that 
depends on the solar wind pressure and that scales the magnetosphere. 
For a given solar wind pressure, rj.g£ also depends on the tilt angle, 
p. Values of r^ef for eight equispaced values of y between 0“ and 
35° are presented for constant solar wind pressure. 

For solar magneto spheric coordinates and for solar magnetic coor- 
dinates separately and for four discrete values of y (0°, 10°, 20°, 
30°), numerical values of AJJ are given for the fields due to 

the magnetopause (0 <. n <. 6) and magnetotail (0 £ n ^ 4) current sys- 
tems. ^ = 0 for all n and m due to symmetry. Further, A^ ^ref^*^ 
is expanded as a power series in y: A^ ~ 

For the field due to magnetopause currents (0 £ n £ 6) and for the 
field due to magnetotail currents (0 £ n £ 4) and for each n and m, 
values of a^, and are presented for each coordinate system. 

IBM card decks containing the coefficients (A^ ^ref^^^) 
tain tilt angles and the fit coefficients (aj, a^, ap, as well as a 
program deck for the computation of the magnetic field at any point 
in space, are available at NSSDC. 
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INTRODUCTION 


The ionosphere is the ionized part of the earth's atmosphere that 
results mainly from the ionizing effect of solar electromagnetic radi- 
ation on the earth's atmosphere. For poleward latitudes, the ionizing 
effect of incident energetic particles is also significant. Tradition- 
ally, the following ionospheric regions and their approximate height 
ranges have been designated: D region (60-90 km) ; E region (90-150 km) ; 
FI region (150-250 km); and F2 region (above 250 km). Recent discus- 
sions of the ionosphere are those of Davies et al. (1968), Hines et al. 
(1965), and Rishbeth and Garriott (1969). 

Of principal interest for ionospheric studies are observations 
of electron density, ionic composition, and temperatures, and the 
dependence of these parameters on position and time. 

Much useful information about these quantities has been obtained 
indirectly by analysis of the effects of the ionospheric plasma on the 
propagation of electromagnetic waves. The presence of the plasma leads 
to wave group velocities less than that of light in a vacuum. Colli- 
sions between neutral and charged particles yield absorption of wave 
energy. The geomagnetic field splits a wave into two waves of differ- 
ent polarization (birefringence) ; these two waves generally propagate 
along different paths. 

Full wave reflection in the ordinary polarization mode occurs when 
a wave of a given frequency reaches that spatial position in the plasma 
at which the plasma frequency [fp, related to electron density N accord- 
ing to N(cm ^) = 1.25 X 10 ® fp^lHz)] is equal to the wave frequency. 
Partial reflections, such as those associated with sporadic E, may 
occur when a wave encounters thin layers of enhanced ionization, steep 
gradients in electron density, or localized scattering irregularities. 
For more complete discussions of absorption, geomagnetic effects, and 
reflection, see the previously cited general references. A recent 
paper contrasting ionospheric measuring techniques and indicating the 
relative advantages of each is that of Booker and Smith (1970) . 


lONOSONDES 


lonosondes utilize the radio wave-reflecting properties of the 
ionosphere. The product of the speed of light in a vacuum (c) and 



half the elapsed time (t) between vertical transmission of a single- 
frequency electromagnetic wave and' reception of the reflected wave at 
the transmitting location is defined as the virtual range (h^) for that 
frequency (f) : 


h' (f) = c t(f)/2 

A sweep -frequency ionogram is a plot of virtual range versus 
frequency and is recorded as nearly instantaneously as possible. A 
fixed -frequency ionogram is a plot of virtual range for a single fre- 
quency versus time. 

In the ionospheric plasma, there is normally a primary maximum in 
the electron density between 200 km and 400 km (F2 maximum) . Except 
for weak backscattering, waves with frequencies greater than the plasma 
frequency associated with this density maximum are not reflected but 
continue propagating through the atmosphere. Thus, the ionosphere is 
divided into two regions, bottomside and topside, relative to the layer 
of maximum density. These two regions can be studied using normal 
ionograms observed only below or above this layer, respectively. 

In addition, sweep- and -fixed -frequency ionograms are distinguished 
according to whether they are A-scan or B-scan ionograms. A-scan ionp- 
grams contain information that indicates the strength of the returned 
signal and hence information about the effects of absorption and partial 
transmission. A-scan ionograms are not produced on a routine basis 
because they are more cumbersome; i.e., they are longer records and 
more difficult to read. The common B-scan ionograms do not nonnally 
contain signal -amplitude information. 


Bottomside Soundings 

Sweep- Frequency Soundings 

A sample daytime bottomside B-scan sweep -frequency ionogram is 
shown in figure 21. (The ordinate is usually called "virtual height" 
instead of "virtual range" for bottomside ionograms.) Such ionograms 
have been obtained with the ionosonde techniques and equipment described 
in Wright, Knecht, and Davies (1957) and by the network of ground-based 
stations listed and illustrated in the WDC-A Catalogue of Data on Solar - 
Terrestrial Physics. A given station typically generates a sweep-fre- 
quency ionogram in about 2 minutes. Five ionograms are normally taken 
per hour, two on the hour, and one at each quarter hour. Microfilmed 
copies of these ionograms and information on time coverage of data of 
individual observatories are available from WDC-A for Upper Atmosphere 
Geophysics, Boulder, Coloradp. 
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NIGHTTIME 


Bottomside nighttime 'B' scan ionogram from Wallops Island, Va. The bottomside 
enhanced regions are normally absent at night. The two hop traces result from 
pulses reflected twice from the ionosphere and once from the ground. 
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DAYTIME 

Bottomside daytime 'B' scan ionogram from Wallops Island, Va. The E, F^, and 
F 2 critical frequencies are indicated. An 'O' (ordinary) and 'X' (extraordinary) 
mode of propagation can usually be identified. 

Figure 21. Sample Bottomside B-Scan Sweep-Frequency lonograms 
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There are some difficulties in interpreting the bottomside B-scan 
sweep-frequency ionograms, and a considerable portion of the literature 
is concerned with interpretation of particular features noted on the 
ionograms. (See Wright, 1967, and Piggott and Rawer, 1961). In most 
cases, some useful quantities may be obtained directly from the shape 
of the virtual height versus frequency curve. These quantities include 
the critical frequencies foF2, foFl, foE, and foEs (ordinary wave fre- 
quencies just large enough to permit transmission through the F2, FI, 
and E regions; a slightly different definition applies to the partially 
reflecting Es level) , the blanketing frequency fbEs (the lowest ordi- 
nary wave frequency at which echo traces are observed on the ionogram) , 
and the virtual heights h^F, h'^Es, and h'E associated with local 
electron density enhancements or maxima. MUF (Maximum Usable Frequency) 
factors, which are useful for oblique propagation computations, can be 
immediately derived from the critical frequencies (Piggott and Rawer, 
1961) . Collections illustrating various types of ionograms have been 
made by Wright and Knecht (1957) and by Shapley (1970). 

F-Plots - Many of the characteristics digitized (scaled) from 
bottomside sweep -frequency ionograms are presented in a graphical for- 
mat called f-plots. For a given station for a given day, the reflection 
heights for a number of characteristic frequencies are presented on a 
single plot. The heights for each characteristic frequency are typically 
given four times per hour. The different characteristic heights are 
distinguished by different symbols . 

Information related to classification of the Es trace is presented 
on a separate small grid usually located below the main plot. A full 
description of this classification scheme and a more complete discussion 
of the f-plot are found in Piggott and Rawer (1961), Wright et al. 

(1957), and Piggott and Bossy (1963). The f-plots are available on 
microfilm from WDC-A for Upper Atmosphere Geophysics, Boulder. Avail- 
ability information, by time and station, is given in the WDC-A Catalogue 
of Data on Solar-Terrestrial Physics . 

Hourly Values and Monthly Siimmary Tabulations - Scaled ionosonde 
data are available in tabular form. Each table presents data taken at 
one station; data points are usually read from the ionograms taken 
nearest the hour. An hourly value table presents data for one charac- 
teristic. The table presents the data in columns for each of the 24 
hours with rows for each day of the month. At the bottom of each 
column, monthly summaries of median and quartile values and the number 
of hourly values from which these values were obtained are often pre- 
sented. These monthly summaries are usually available separately in 
tables using columns for each hour of the day and rows across the tables 
for the various characteristics . Qualifying and descriptive letters 
used in the hourly value and monthly stnumary axe defined and discussed 
in Piggott and Rawer (1961) and are often contained in brief form 
in the introductions to published tables of data. 
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These digitized ionospheric data tables are available for more than 
200 soTinding stations. For recent years, tables for over 100 stations 
are available. The most complete and convenient source for these tables 
(in both hard-copy and microfilmed forms) is WDC-A for Upper Atmosphere 
Geophysics, Boulder. Data on magnetic tape for over 60 stations are 
also available from the same source. Detailed availability information 
is given in the WDC-A Catalogue of Data on Solar -Terrestrial Physics . 
WDC-A personnel can often be of assistance in obtaining data that are 
not currently available in their files. 

Electron Density Profiles - Profiles of electron density (N) versus 
real (geocentric) height can be calculated from ionograms using the 
equation defining virtual height [h‘'(f)]. 


h'(f) = I n'[N(h), f, B(h), 0(h) v]dh 

J o 


where h" and hj^ are virtual and real heights of reflection of the wave 
with frequency (f) . The group index of refraction is designated by n'; 
N(h) is the electron density, B(h) is magnetic field magnitude, 0(h) is 
the angle between the magnetic field direction and the wave propagation 
direction, and v is the collision frequency. When the changes in mag- 
netic field strength along the propagation path are small and the 
collision frequency is small (true above the D region), the equation 
can be simplified to 


h"(f) 



n'[N(h), f]dh. 


A vertical propagation path is usually assumed. 

There are several techniques for solving this equation to obtain 
the N(h) profile from observed values of h'. The frequently used 
lamination technique involves a step-by-step procedure. It is assumed 
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that N(h) is known for h i h^. For a wave frequency fj, reflected 
from a known electron density level at a small but unknown distance 
Ah above h^, the virtual height is 

fhi rhj^-Ah 

h"(f2) = n'[f^, N(h)]dh + n'[f2» N(h)]dh. 

Jo J hj 

The quantity h" (f^) is read off the ionogram, and the first integral is 
known because N(h) is known for ^hj. in the second integral, both 
N(hj) and N(h^ + Ah) are known, although Ah is not itself known. By 
assuming a functional form for the variation of N(h) with h in the small 
height interval h^ to h^ + Ah, Ah may be determined as that value yielding 
the required equality. This completes a given step. 

Complicating effects in this approach are: (1) low level absorp- 

tion, which results in having to assume a. value for hj for the first 
step; (2) an E region electron density maximum; and (3) determination 
of h" at the absolute density maximrnn where the h'‘(f) curve becomes 
asymptotic to the vertical. A much more complete discussion of these 
effects and of details describing other reduction techniques can be 
found in "Special Issue on Analysis of lonograms for Electron Density 
Profiles," Radio Science , volume 2 (new series) #10, October 1967. 

An extensive bibliography on profile calculation is included in this 
issue. 

In addition to the generation of N(h) profiles from individual 
ionograms, a composite N(h) profile can be generated from a "median 
ionogram." Twenty-four median ionograms are generated (one for each 
hour) for a given station each month. 

For a given N(h) profile, a number of parameters have been defined 
(Wright, 1967, and Piggott and Rawer, 1961). These include: 

1 . h (height of the maximum) - real height of the F layer 

maximum assuming a parabolic profile near the maximum 
(same as the more commonly used h^, or HMAX) . 

2. q (quarter thickness) - half the height difference 
between h and the two points at which the electron 
density ii three-quarters of that at h^. 

3. SCAT - same as q ; also called characteristic thickness. 


108 



4. HMIN, HMAX - true heights from which the minimum and maximum 
frequencies are reflected. 

5. NMAX - electron density at h^. 

6. SHMAX - total electron content (TEC) downward from h^,. 

7 . SHINF - TEC from the ground to infinity, assuming a 
Chapman layer density profile above the density maximimi. 

8. RATIO - ratio of TEC to SHMAX. 


Hard-copy tabular data, which contain ionogram station and time 
information, electron densities, virtual and true heights, total co- 
lumnar content below given heights, and some of the parameters listed 
above, are available at WDC-A, Boulder. The WDC-A Catalogue of Data 
on Solar-Terrestrial Physics includes the station list, data availa- 
bility, and sample formats. These tables have been prepared by the 
Institute for Telecommunication Sciences and Aeronomy (ITSA) , NOAA, 
from ionograms (individual and median) of 41 stations for various 
periods of time since 1956. The format of the tables is given in a 
paper by Wright and McKinnis (undated). This paper is available from 
the authors or from NSSDC. 


Oblique Soundings 

The oblique incidence sweep -frequency ionospheric sounding tech- 
nique uses the same principle of operation as the vertical incidence 
sounder (Davies, 1964). The primary difference is that the transmitter 
and receiver are located at different places so that the signal is 
reflected at an angle. The most common use of these sounders is to 
operationally investigate optimiom conditions for radio propagation. 
According to Rawer (1969) , no attempts are known to have been made 
to compute electron density versus height profiles from oblique inci- 
dence ionograms. These oblique incidence soundings do not adequately 
substitute for the vertical soundings. A very limited number of oblique 
incidence soundings are available through the WDC-A for Upper Atmosphere 
Geophysics. This same WDC-A subcenter can assist in the identification 
of currently operating stations which may have observational records 
that can be made available. 
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Fixed- Frequency Soundings 


Fixed -frequency bottomside sounders (also called virtual height 
recorders) are in operation at some ground ionospheric observatories. 
Recordings are not generally available except through the individual 
stations. These recordings normally operate continuously and are used 
to monitor ionospheric changes . 


Topside Soundings 


Sweep- Frequency Soundings 

A sample topside sweep -frequency ionogram is shown in figure 22. 
Such ionograms have been obtained by equipment and techniques described 
in Franklin and Maclean (1969) and Franklin et al. (1969). These tech- 
niques and the equipment were used on the Alouette 1 (launched 9/62) , 
Alouette 2 (launched 11/65), and ISIS 1 (launched 1/69) satellites. 

All are near -polar satellites and were operating as of January 1971. 
Observations have been taken in the vicinity of as many as 22 telemetry 
stations, and a limited number of ionograms from ISIS 1 have been tape- 
recorded at other locations. The availability of observations depends 
on many factors, including: 

1. Location of telemetry station, 

2. Satellite altitude and other orbit characteristics. 

3. Satellite operating conditions (power supply, etc.). 

4. Satellite scheduling at telemetry stations. 

For telemetry reception with elevation angles ^ 15°, the observing 
range of a satellite at a height of 1000 km is within a radius of 2100 
km from the station. For some stations, elevations of 0° are not un- 
common, with a resulting observing range of nearly 3400 km for satel- 
lites at a height of 1000 km (Franklin et al . , 1969) . Of about 24 
stations which have received the telemetered data, all but a few lie 
near the 90th geographic meridian; these stations are fairly well 
distributed in latitude. The satellite sounders have generated one 
ionogram every 20 to 30 seconds for up to 9 hours per day . For Alouette 
1, in a circular orbit at about 1000 km, this means typically one iono- 
gram (20 seconds per frame) for each degree of latitude along the orbit 
when in range of, and scheduled by, a telemetry station. For ISIS 1 
and Alouette 2, with an orbit perigee/apogee of about 500/3000 km, it 
means typically one ionogram (30 seconds per frame) for each 2° of 
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Figure 22, Sample Topside Sweep - Frequency lonogram 



latitude near perigee and one for each 1.1® of latitude near apogee. 
Microfilmed copies of these ionograms are available from NSSDC approxi- 
mately 2 years after observation. These ionograms can be obtained also 
through the WDC-A for Rockets and Satellites or the WDC-A for Upper 
Atmosphere Geophysics. Current listings of specific times and places 
for processed observations (i.e., prepared ionograms) are also main- 
tained at NSSDC. As with bottomside ionograms, interpretation of top- 
side ionograms requires experience, skill, and training. (Hagg et al., 

1969) . 

The appearance of the ionograms from each of the three sweep- 
frequency topside sounders is different. Sweep rates are all the same 
above 2 MHz (i.e., about 1 MHz/sec), but below 2 MHz, the rate for 
Alouette 2 is 0.13 MHz/sec and for ISIS 1 is 0.31 MHz/sec. This causes 
the frequency marker intervals below 2 kHz on the ionograms to be 
stretched out on the two later satellites, with a resulting change in 
appearance of the ionograms . The frequency ranges are also different 
for each satellite: Alouette 1, 0.5 to 12 MHz; Alouette 2, 0.12 to 

14.5 MHz; and ISIS 1, 0.1 to 20 MHz. An added feature on the ISIS 1 
ionograms is a short record of fixed -frequency sounding normally ap- 
pearing prior to each sweep-frequency ionogram. (See figure 23.) No 
enhanced electron density regions, except for the F2 maximum, can be 
identified regularly on topside ionograms. 

Electron Density Profiles - The principles used in obtaining electron 
density profiles from topside ionograms are the same as those used for 
bottomside ionograms. Reduction of topside ionograms has been done primarily 
by the lamination method described by Jackson (1969) or by some modifi- 
cation of this technique. Important differences which must be taken 
into account are (1) that errors due to non-vertical propagation cannot 
be assumed to be small when the sounder range starts to exceed 1000 km 
(bottomside sounder range seldom exceeds 500 km) and (2) that the X trace, 
in contrast to the 0 trace for bottomside reductions, is normally used 
to make reductions from topside ionograms. Use of the X trace also 
makes the computation more complex since its use requires consideration 
of the magnetic field and its variations along the propagation path. In 
addition to Jackson (1969), other primary references treating topside 
reduction are. Jackson (1967) , Fitzenreiter and Blumle (1964) , and Thomas 
et al. (1965). 

Of the approximately 2.5 million satellite ionograms that have been 
produced to date by the two Alouette satellites, it is estimated that 
about 5% to 7% will be reduced. At present, reductions are being made 
as needed to support research projects of the major ISIS project par- 
ticipants including the Communications Research Center (CRC) , Ottawa; 
the Radio and Space Research Station, Slough, Bucks, England; NOAA, 

Boulder; NASA-Goddard Space Flight Center, Greenbelt, Maryland; and 
NASA-Ames Research Center, Moffett Field, California. A limited 
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number of these reductions prepared by the CRC and the Radio and Space 
Research Station are available in hard copy and on tape from NSSDC or 
directly from the agency (see references 2, 3, and 28). The NASA- 
Ames Research Center facility has reduced satellite ionograms on a 
mass production basis, and a quantity of Alouette 1 profiles are avail- 
able in hard copy (Chan et al . , 1966; Thomas et al., 1966; and Colin 
and Chan, 1967) or on tape from NSSDC. About 8000 Alouette 2 ionograms 
have also been reduced by the Ames Research Center, and it is expected 
that these will also be available on tape soon. Since October 1967, 
a Japanese telemetry station has been taking observations at Kashima, 
Japan. Their first reductions for October to December 1967 are now 
available in hard-copy form, and continuation of routine reduction of 
Kashima ionograms is expected (Ref. 12) . 

Scaled Digital Data - Alouette Synoptic (Alosyn) data have been 
tabulated from a large number of the Alouette 1 ionograms recorded at 
11 of the primary telemetry stations. They contain satellite position 
information and related geophysical data such as zenith angle and elec- 
tron gyrofrequency based on a model geomagnetic field. Listed with 
each satellite position are observed values of plasma frequency at the 
sounder at Fj penetration and any observed sporadic E. For topside 
sounding, the X trace is normally used in scaling, but both ordinary 
(0) and extraordinary (X) trace observed values, as well as 0 values 
computed from the X values, are included in the Alosyn data. Ground 
echo occurrence and intensity have also been listed. Alouette 1 Alosyn 
data are available from launch (September 29, 1962) through January 1965 
and for a few months in 1967 and 1968. It is anticipated that additional 
listings will be prepared by CRC. Additional information and/or copies 
of these data are available on tape, microfilm, or hard copy through the 
experimenter (CRC) or NSSDC. 

Fixed-Frequency Soundings 

A typical fixed -frequency topside ionogram is shown in figure 23. 
These ionograms have been obtained by equipment and techniques described 
by Russell and Zimmer (1969). There have been two satellite-borne fixed- 
frequency sounders, one on Explorer 20 and the other on ISIS 1. Explorer 
20 sampled each of six frequencies sequentially, and, in processing, 
records of each frequency were placed on separate ionograms. About 1000 
hours of soundings made during 1964 and 1965 near telemetry stations are 
available on microfilm. ISIS 1, still operational, sounds for 3 to 5 
seconds at one of six possible frequencies prior to each sweep -frequency 
sounding. These are recorded ahead of the sweep -frequency ionograms on 
microfilm and have been observed near telemetry stations since January 
1969. These data have been used to study plasma resonances and small- 
scale ionospheric structure. No digital scaled data are available. 
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Figure 23. Typical Fixed -Frequency Topside lonogram 


INCOHERENT BACKSCATTER RADAR OBSERVATIONS (THOMSON SCATTER) 


For radio waves transmitted through the ionosphere with frequencies 
much larger (50 to 1300 MHz) than the plasma frequency (up to 15 MHz) , 
very small-scale fluctuations of electron density in the plasma will 
scatter a portion of the signal back toward the transmitter. These 
fluctuations are caused by the random thermal motions of the electrons 
(Evans, 1969, and Davies, 1965). The height at which the scattering 
occurs can be determined by either a pulse technique or by the use of 
two intersecting cylindrical radar beams. The height resolution of a 
backscattered pulse or beam varies from 4 to 35 km at a height of 100 
km to about 300 km at altitudes above 1000 km, depending on the equip- 
ment used and the geographical location of the station. 

From analysis of the backscattered signals , electron density pro- 
files, electron and ion temperature profiles, and ion composition and 
ionospheric drift can be determined. The equipment used for the obser- 
vations determines the vertical extent of the backscatter profiles and 
the resolution, as well as influencing the procedures used to convert 
the scattered return to measures of meaningful ionospheric parameters. 

Radar installations that can make these observations are complex 
and expensive to operate; consequently, these observations are not 
scheduled on a routine basis. Rather, they are scheduled on an "as 
needed" basis , and the equipment is shared for other uses . The 
available data can be obtained by inquiry to the observatory or to the 
WDC-A for Upper Atmosphere Geophysics . The seven stations taking inco- 
herent backscatter radar observations are located in the U.S., Canada, 
France, and England. The stations are listed and described in detail 
by Evans (1969) . 


IONOSPHERIC ABSORPTION OF RADIO WAVES 


Absorption of radio waves occurs when electrons responding to the' 
wave fields collide with and transfer energy to the neutral particles. 
A study of ionospheric absorption yields information on both the D and 
E regions, where both electrons and neutral particles exist in suffi- 
cient quantity to produce significant absorption. Knowledge of D 
region characteristics is unattainable from normal sweep-frequency 
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ionograins because absorption prevents the return of a detectable re- 
flected signal. In practice, absorption characteristics are studied 
primarily in connection with radio communications applications. 

Three techniques of observing ionospheric radio wave absorption, 
identified as Al, A2, and A3, have been designated by the International 
Council of Scientific Unions. They are (Al) , radio waves transmitted 
vertically and received at the same station (vertical ionosondes), (A2) , 
waves received from extraterrestrial sources by riometers , and (A3) , 
waves transmitted obliquely and received by spatially separated sta- 
tions. Discussions of ionospheric absorption theory, techniques, and 
results are found in Geisweid et al. (1963, 1968), Piggott et al. 

(1957), Allcock (1954), and Mitra (1970). 


Pulse Reflection Method (Al) 

An amplitude versus time oscilloscope (A scan) display of a single- 
frequency vertical ionosonde pulse is used to measure the amplitude (I) 
of an echo pulse. Comparison of this value with the amplitude that 
would have been returned in the absence of collision losses (lo) provides 
a quantitative scale with which to measure absorption. This ratio is 
called the apparent reflection coefficient (p) 



The absorption loss (L) is measured in decibels . 

L = -20 logp 

The reference value lo is determined by nighttime observations, usually 
a mean of values over night observations for several different dates. 
The noon observation of L is most frequently available and is generally 
considered representative of absorption for that day. 

It is usually difficult to determine the relative contributions 
to L arising from the D and E regions . 

Data are available from the WDC-A for Upper Atmosphere Geophysics 
for 37 stations, 14 of which have records that cover more than 5 years. 
Data for a given station and month appear on a single page. Local 
observing time (usually noon) and values of absorption for each day 
and observing frequency (usually one value near 2.2 MHz and at least 
one other value) are listed. Monthly median values for each frequency 
are also listed. For further details, see the WDC-A Catalogue of Data 
on Solar -Terrestrial Physics . Symbolism used in recording these Al 
observations has been standardized (Piggott et al., 1957). 
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Riometer Method (A2) 


A riometer (Relative Ionospheric Opacity Meter) is a detector that 
measures the absorptivity properties of the ionosphere. This direc- 
tional detector is oriented in a fixed direction relative to the earth 
(usually vertical) and measures the incident cosmic radio noise at a 
single frequency between 8 MHz and 120 MHz (usually 20, 30, or 60 MHz). 
During the course of each sidereal day, the riometer receives noise 
from all celestial longitudes ($) . On days of negligible ionospheric 
disturbance, the riometer signal power will be Po($). Due to the gen- 
eral constancy of cosmic radio noise, Po($) will be constant except 
for a large-scale, small-amplitude, daily solar effect which regresses 
through the observations at the rate of 1 day per year (4 minutes earlier 
each day) . Deviations of observed riometer signal power (P) from Po 
are primarily attributable to ionospheric variations particularly in 
the D and E regions . This method is particularly useful during radio 
blackout (high absorption) since it involves the use of higher frequency 
signals that are not completely absorbed. More detailed discussions 
of riometers are found in Little (1959) and Little and Leinbach (1957) . 

Roll chart recordings of riometer output current versus time obtained 
at over 60 stations are available in hard-copy and microfilm forms from 
the WDC-A for Upper Atmosphere Geophysics. Availability information, 
by station and time period, is found in the WDC-A Catalogue of Data on 
Solar-Terrestrial Physics . 

Absorption (A) in decibels is defined as A = -10 log (P/Po) . Tables 
of the maximum absorption reading for each hour are also available in 
hard-copy and microfilm formats at WDC-A, Boulder. No conversions are 
routinely made from absorption data to total electron content or changes 
in content. 


Continuous Wave (CW) Field Strength Method (A3) 

Method A3 consists of comparing signal strengths of radio signals 
obliquely reflected from the ionosphere during times when absorption is 
present with similarly obtained signal strengths for times when absorp- 
tion is absent (Geisweid et al., 1963, 1968). Criteria for observing 
practices depend upon many factors . The normal operating range between 
the transmitter and receiver is greater than 200 km in order to avoid 
any contribution to the received signal from the ground wave, and the 
distance is normally less than 500 km so that the reflection from one 
level (E region) is much stronger than from any other level. Signal 
strength observations from night sporadic E (or adjusted values from 
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the F region) reflections are used for reference "no absorption" values, 
against which the day values are normalized. The absorption is given 
by 

L = -(cos a) 20 log 

b 1 
n n 

where a is the incident angle, and Ejj are the received day and night 
signal strengths, and l^ and 1^ represent the signal path lengths for 
day and night observations. The absorption occurs in the D and E re- 
gions; again, separation of D and E effects is very difficult. As with 
the A1 method, relatively small changes in electron density can be de- 
duced from these observations, but with large enhancements of electron 
density, the signals may be entirely absorbed and no reflections de- 
tected. Data are available for 19 receiving stations. Unsealed data 

are available from 10 stations, all for periods of over 6 years. 

Digitized values are available from the other nine stations. These 
are tabulated values for day of month versus time of day on each sheet. 

Six of these nine stations have records available for periods of over 

5 years. 


"F min" as a Measure of Absorption 


Another measure of absorption that is sometimes used is the f min 
value scaled from the ionogram (Piggott et al., 1957). Caution must be 
used when interpreting f min as absorption since changes in recorder 
gain and variation in noise level also affect this value. The f min 
value is the lowest frequency at which combined recorder gain, noise, 
and absorption permit echoes of the ionosonde transmitter to be returned 
and detected on the ionosonde receiver. These data are readily avail- 
able from WDC-A, Boulder, for most ionosonde stations as one of the 
scaled parameters for which both daily values and monthly medians are 
available (see page 106) . 


SATELLITE BEACON EXPERIMENTS 


Changes in radio frequency signal characteristics in the range of 
15 to 150 MHz result when a radio signal traverses the ionosphere. 

The signal changes can be studied to deduce ionospheric characteristics 
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and their variations (Mass, 1963, and Garriott et al., 1970). Several 
arrangements have been used to measure these signal variations, includ- 
ing reflection of ground transmissions from the moon or a satellite, 
transmission from the ground to a satellite, and the most commonly used 
arrangement, transmission from a satellite beacon to the ground. A 
synchronous altitude (approximately 35,800 km) satellite transmission 
is convenient for the study of time changes of ionospheric conditions 
at any given location. Transmissions from other satellites can be used 
to obtain wide geographical coverage, but numerous ground stations are 
required for obtaining wide coverage. 

The signal characteristics affected by the ionosphere are signal 
strength, frequency, direction of arrival, and polarization angle. 

The observations of variations in these characteristics are interpreted 
in terras of ionospheric total electron content, scintillations, and 
irregularities. Determinations of the changes in time and space of 
these ionospheric features, as well as the processes related to these 
changes, are the objectives of these satellite beacon studies. The 
results of these studies are useful in the design and operation of 
radio communications and tracking equipment. 

Ionospheric effects on radio signals increase rapidly with decreas- 
ing radio frequency. However, signals with frequencies lower than the 
F layer critical frequency will not penetrate the ionosphere but will 
be reflected. Since F layer critical frequencies range from 5 to 15 
MHz (depending primarily on time of day, location, and solar cycle 
phase), optimum frequencies for beacon experiments are around 20 MHz, 
with a higher frequency sometimes used concurrently for comparison. 
Transmissions that have not been specifically designed for use in 
ionospheric study have often been successfully employed; one example 
is the Syncom tracking beacon, which transmitted at 136 MHz. 

One major difficulty with all beacon observations is that the 
propagation path does not pass perpendicularly through the ionosphere. 
Even if the satellite were at zenith, tilts in the ionospheric layer, 
multipath propagation, refraction, etc., would make it difficult to 
normalize observations to some standard path. Nevertheless, some nor- 
malizations or assumptions must be made in order to make comparisons 
between sequential observations. It is common practice to use a simple 
geometric adjustment in allowing for the incident angle of the wave, to 
assume a horizontally stratified ionosphere, and to limit observations 
to zenith angles of less than 40“ in order to avoid gross errors . 
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Total Electron Content (TEC) 


TEC is the total number of electrons in a column of unit cross 
section extending from the transmitter to the receiver. TEC is calcu- 
lated from observations of the rotation of a polarized signal (Faraday 
rotations) transmitted from the satellite. TEC can also be computed 
from observations of the dispersive Doppler frequency shift of the 
beacons due to the ionosphere (Mass , 1963) . There are several specific 
techniques for obtaining TEC by using Faraday rotation and/or Doppler 
shift observations. Estimates of accuracy for these data range from 
2% up to 30%, depending on the technique and on the validity of assump- 
tions for the particular observational conditions. Digitized tabula- 
tions of TEC data taken in the vicinity of several stations for various 
time periods between 1960 and 1970 may be obtained from the WDC-A for 
Upper Atmosphere Geophysics, and its catalog gives information on the 
availability of TEC data. Yuen and Rpelofs (1970) contains a large 
number of diurnal TEC profiles using the synchronous altitude ATS 1 and 
Syncom 3 satellite located at the equator over the Pacific Ocean. The 
data, which cover the period 1964 through 1968, are considered valid for 
a sub-ionospheric point a few hundred kilometers south of Hawaii. 


Scintillation 


Scintillations are rapid fluctuations in the amplitude, frequency, 
polarization, or direction of an observed signal. Scintillation of 
beacon signals is caused by small-scale variations of electron content 
that occur along the propagation path. These irregularities in the 
structure of the ionosphere are usually considered to extend less than 
1 km. The scintillations are commonly observed as noise-like structures 
(rapid signal strength variation) on strip chart recordings of satel- 
lite tracking records or Faraday rotation recordings. Tracking re- 
cords of selected equatorial and polar stations are available for a few 
satellites from WDC-A. The intensity and the start and stop times of 
scintillations occurring on records of signal polarization changes 
(Faraday rotation) are sometimes noted in the tables of TEC. 


Irregularities 


Ionospheric irregularities are generally considered to represent 
middle-scale variations of total electron content, i.e., those with 
horizontal dimensions from 1 to 1000 km. Irregularities are sometimes 
identified by rapid changes in TEC which show strong horizontal elec- 
tron density gradients in the ionosphere. The large areas in the 
ionosphere in which scintillation occurs are also referred to as 
irregularities . 
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Joint Satellite Studies Group - An informal organization of scien- 
tists who are particularly interested in using beacon equipment and 
data has been formed and is called the Joint Satellite Studies Group. 

If more detailed information about the availability of beacon data is 
desired, it can be obtained through the chairman of this group. Dr. Jules 
Aarons, AFCRL, Bedford, Massachusetts 01730. 


IONOSPHERIC DRIFTS 


Ionospheric motion is a complex phenomenon. At ionospheric heights 
of 100 to 500 km, less than 1% of the particles are ionized. Ionospheric 
motions, then, may be expected to respond to neutral motions as well as 
to electrical fields and currents, magnetic fields, diffusive processes, 
etc. Ionospheric motions of various scales of magnitude have been ob- 
served. As for most ionospheric observations, remote sensing techniques 
are used also for observing ionospheric drift. Some identifiable iono- 
spheric feature is tracked, and the motion of that feature is then as- 
sumed to represent the "drift" motion of the ionosphere. 

Several papers in the literature compare observations from differ- 
ent observing techniques and attempt to give physical reasons for dif- 
ferences in results. There are five well-known drift observing tech- 
niques. A considerable amount of the data that are readily available 
have resulted primarily from the fading method described below. IGY 
data (1957 to 1959) from over 20 stations are available at the WDC-A 
for Upper Atmosphere Geophysics. In subsequent years, about 10 stations 
have routinely provided data to WDC-A. Detailed availability of data 
is given in the WDC-A for Upper Atmosphere Geophysics data catalog. 

Rawer (1965) contains summary data of fading method observations from 
10 stations for the years 1957 to 1959 . The drift tracking methods , 
along with comments on the procedures and availability of data, are 
outlined below. For a more detailed survey on ionospheric drifts, see 
Beynon and Brown (1957), Briggs (1957), Huxley and Greenshow (1957), 

Rawer (1968, 1969), Spencer (1955), and Spizzichino (1968). 


Fading Method (Dl) 


Radio pulses are transmitted upward, and their downcoming reflec- 
tions are monitored by three receivers spaced some distance (about one 
wavelength) apart. Pulsed transmissions at fixed frequencies between 
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2 and 10 MHz are most frequently used. Signal fading that occurs along 
the propagation paths will appear on the amplitude records of the re- 
ceived signals. When similar fading features are noted on all three 
records of the received signal, it is assumed that that particular 
signal results from the same ionospheric feature. The differences in 
time of reception of the commonly returned signal occurrence at the 
three receivers can be interpreted as twice the velocity of the iono- 
spheric motion near the reflection height. It is assumed that (1) the 
cause of fading is due to small-scale variations in electron density 
near reflection height, (2) these variations are conservative over short 
time periods, and (3) time changes in these small-scale variations rep- 
resent general horizontal motion of the ionosphere. 

The altitude of the observed drift is a function of both electron 
density and frequency; thus, for a given frequency, the height of drift 
information will show considerable diurnal variation and occasional 
shorter term storm variation. Data are limited to the bottomside iono- 
sphere (80 to 350 km) . 

Data listings for each station consist of observation time; speed 
of the irregularity, v (one half the observed ground speed) (Wright, 
1968); azimuth toward which the drift moves, ()); and component velocities, 
+Vjj and +Vy. The listed values are a composite of values obtained from 
many individual observations and thus represent a mean motion over a 
period of minutes. There are several methods for converting observations 
to drift motions (Briggs, 1957; Rawer, 1969, p. 123-25; and Pfister and 
Bibl, 1968) , but all result in data listings of similar form. Most of 
the hard-copy data readily available in WDCs are derived from this 
fading method. A large amount of data from Breisach, Germany, has 
been summarized recently (Hamischmacher and Rawer, 1968) . 


Radio Star Scintillation Method (D3) 


Radio signals from radio stars are received at a selected frequency 
between 30 and 100 MHz. Three receivers are spaced at about 4 km or less, 
and signal strength fluctuations that are similar on all three receiver 
records are presumed to be caused by the same ionospheric feature. In 
a manner similar to the fading method, the time delays can be converted 
to drift motions given in speed, V, and azimuth, <j>. The altitude of 
drift has been associated with the F layer (Little and Maxwell, 1951), 
but it is possible for the altitude range affecting these observations 
to extend over a considerable range of both the bottom and top sides of 
the ionosphere. The observation itself cannot resolve the drift altitude 
accurately. Observations by this technique have not been extensively 
made, and no records are available at WDCs. Observations were made 
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in the early 1950 's at Jodrell Bank, University of Manchester; the 
University of Cambridge; and Mayaguez, Puerto Rico (Olivers, 1961; 
Dueno, 1961; Hewish, 1952; and Maxwell and Little, 1952.) 


Meteor Trail Method (D2) 


The ionized trails created by collisions of meteors with neutral 
particles between 80 and 100 km are tracked, and the tracking data are 
converted to speed and direction values . Two methods have been used, 
one a continuous wave technique and the other a coherent -pulse tech- 
nique. Both techniques involve similar principles: a radio frequency 

signal is directed toward a selected area of the sky, and the portion 
of the signal that is scattered back from the trail is "beat" with a 
known signal; these beat characteristics are used to determine range 
and radial motion of the trail. 

Varying the transmitted signal Characteristics and observing in 
two different areas of the sky at a zenith angle of about 45® provides 
adequate additional information to obtain horizontal drift motions at 
relatively accurate determined altitudes. The smaller vertical com- 
ponents of motion are sometimes computed. For the continuous wave 
technique, the trail must persist for 0.5 to 3 seconds, and the trails 
must occur frequently enough to allow an average of four to five obser- 
vations per hour to be made. For the pulsed observations, trails per- 
sisting for less than 0.1 second can be used, and hundreds of such 
meteor trails per hour normally are available for making observations. 
No extensive files of meteor trail data are readily available, but 
obseirvations have been made at the following stations (Bradley et al., 
1970) : 


Australia, Adelaide, University of Adelaide 

Australia, Newcastle, University of Newcastle 

Canada, Ottawa, National Research Council 

Czechoslovakia, Prague, Ondrejov Astronomical 
Observatory 

East Germany (GDR) , Kiihlungborn, Institut fur 
lonospharen Forschung 

France, Garchy, National Center for Telecommunications 
Studies (CNET) 
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Sweden, Lund, Lund Observatory 

U. K., Sheffield, University of Sheffield 

U. K., Malvern, Royal Radar Establishment 

U. S. A., Alaska, College, University of Alaska 

U. S. A., California, Stanford, Stanford University 

U. S. S. R., Kazan, Kazan University 


Characteristic Reflection Method (D4) 


Records of unique reflection patterns from ionospheric soundings 
can be recognized as the same pattern at different receiving station 
locations. When time differences between station passage of these ' 
patterns can be determined, an ionospheric drift velocity vector can 
be computed. The most useful arrangement is for a single station to 
receive pulsed signals that were sent at the same sounding frequency 
by three or more transmitting stations. Continuous records of signal 
transit time versus time are kept for each of the three stations. Many 
kinds of patterns that are returned can be recognized as originating 
from the identical ionospheric feature. An ionospheric drift can be 
calculated from the time differences between pattern occurrence and 
distances between transmitters and receivers . Patterns from sweep- 
frequency ionograms spaced at greater distances have also been used 
in a similar way. No long periods of such drift observations have 
been prepared using this technique. 


Chemical Release Methods (D5) 


Both luminous and illuminated (by twilight) chemicals that are 
released from rocket vehicles can be used to observe drift velocities. 
An identifiable part of the cloud is optically tracked, or sequential 
photographs are taken. Chemicals may be released either as a trail or 
in discrete blobs. The altitudes are known from rocket trajectory in- 
formation. Some of these wind data for mid latitudes have recently 
been summarized (Rosenberg, 1968) . 
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IONOSPHERIC EVENTS 


There are a variety of observable effects accompanying ionospheric 
events that are commonly referred to by their abbreviations (Lincoln, 
1968; Warwick, 1963; Hultqvist, 1963; and Solar- Geophysical Data , 1970). 
The events are related to enhancements of electromagnetic radiation, 
energetic particle flux, and plasma flow (and the resulting geomagnetic 
storm), .associated with solar activity. The observable effects result 
from enhanced ionospheric electron density, mainly in the D region; 
these enhancements are due to photo -detachment and/or collisions. 


PCA 


Polar Cap Absorption refers to D region electron enhancement over 
the entire polar cap and is caused by incoming solar protons and parti- 
cles with energies of 1 to 100 Mev. A PCA event begins 15 minutes to 
several hours after a flare occurrence, and it may persist for 1 to 6 
days. PCA was first identified in 1956 from VHP scatter data, vertical 
sounder ionogramsv-and riometer records. An early criterion for PCA 
was shortwave blackout (i.e., absorption below 3 to 4 MHz on vertical 
sounder records), which persisted for over 3 hours and covered a con- 
siderable polar geographic area. More recent definitions utilize a 
given absorption loss. 

There is no standard definition for PCA; thus, any particular list 
of events will reflect definitions that were adopted for its preparation. 
Lists available from the literature are noted on page 193 of Hultqvist 
(1963) . Routine listings by the WDC-A for Upper Atmosphere Geophysics 
are expected soon; they will likely include UT, time of maximum, start 
and stop, technique, frequencies, locations, remarks, and observations 
of related phenomena. An index of auroral blackout (Ip) is described 
below. 


AZA 


Auroral Zone Absorption, also known as "auroral blackout" or 
"auroral absorption" (AA) , is delayed 20 to 40 hours from the associated 
solar disturbance (not necessarily a flare) and is caused primarily by 
precipitating electrons with low energies between 50 and 500 kev. As 
with PCA, high frequencies (3 to 30 MHz) are strongly absorbed. These 
occurrences normally persist for periods ranging from 1 to 6 hours, and 
the blackout zone is generally concentric with, and a few degrees equa- 
torward of, the auroral oval. Since identification of AZA depends upon 
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information other than that contained in any one particular observational 
record (i.e., geographical distribution of absorption, delay time from 
disturb^ce, etc.)» these events are most frequently included in listings 
of PCA events . They are sometimes referred to as a type of PCA and are 
presently not identified separately in any known listings . Frequencies 
absorbed are the same as for PCA; only the geographical area covered, 
the time relations to flare occurrence, and species and energies of 
incident particles are different. 

PCA and AZA blackout indices. Ip and la, have been discussed by 
Hakura et al . (1967). These daily blackout indices provide a measure 
of the temporal extent of blackout for a given day. The index values 
run from 0 to 9 and are related to the number of hours (H) of the day 
during which blackout occurred, as shown in the following table. 


Ip or la 

0 

1 

2 

3 

4 

5 

6 

7 

8 
9 


Hp or Ha 

0 . 0 - 0. 4 
0.5-0. 9 
1 . 0 - 1. 5 

1 . 6 - 2. 5 

2.6- 3. 5 

3. 6- 5. 5 

5.6- 8. 5 

8.6- 13.0 

13.1- 20.0 

20.1- 24.0 


These indices are tabulated for the period 1957 to 1965 in Annals 
of the IQSY , volume 2, and in the "Abbreviated Calendar Records" from 
1964 to the current issue. (See Part VII of this Handbook . ) 


SID 


The Sudden Ionospheric Disturbance is usually attributed to solar 
flare electromagnetic radiation incident on the ionosphere. SIDs are 
observed on the sunlit hemisphere and occur almost simultaneously with 
visual flare observations. SIDs and subsequent recovery of the iono- 
sphere have a time duration somewhat longer than flare duration; gen- 
erally from minutes to an hour, with rise more rapid than decay. 

Several of the various SID observations are routinely identified 
and digitized. Samples of several types of original records as they 
appear before digitizing are nicely illustrated by Lincoln (1968). The 
types of SID are alphabetically listed below with brief descriptions of 
each. The SCNA, SEA, SES, SFD, and SWF listings are available from the 
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WDC-A for Upper Atmosphere Geophysics and have also been published in 
issues of Solar-Geophysical Data since 1963. Bursts and SFE are also 
listed in Solar-Geophysical Data . Information available in the data 
listings includes, in addition to that noted for PCA, the "definiteness" 
rated on a subjective scale. 

Burst - A solar noise burst at riometer frequencies is a sudden 
increase of signal strength on the riometer records . This is caused 
by solar radio noise emissions at riometer frequencies . This noise 
is affected by and can mask the ionospheric absorption that is present. 
The bursts are of such short duration, however, that absorption changes 
during the burst can be estimated. In addition to riometer observations, 
solar noise is routinely observed at much higher frequencies than those 
used for riometer receivers. 


Crochet - See SFE. 

SCNA - Sudden Cosmic Noise Absorption, also referred to as CNA 
(Cosmic Noise Absorption) Type I (see below), is characterized on 
riometer recordings by a sudden daytime absorption increase of several 
db within a few minutes. Decay is somewhat slower than rise, and the 
riometer record is perturbed for a period of minutes to over an hour. 

SCNA involves absorption outside of auroral regions. 

SEA - Sudden Enhancement of Atmospherics is observed as an increase 
in signal strength on wideband equipment operated to detect electromag- 
netic emissions from lightning. This equipment operates in the VLF 
(10 to 50 kHz) range with the most commonly used frequency near 27 kHz. 

SES - Sudden Enhancement of Signal is observed on VLF frequencies 
(15 to 50 kHz) . These observations are nearly identical to SEA except 
that the receivers are narrow-band receivers designed to pick up man-made 
VLF transmissions. As with SEA, signal strength increase is the SID 
indicator. 

SFD - Sudden Frequency Deviation, caused by ionospheric disturbances, 
is observed by beating the received radio signal with a reference fre- 
quency that is near the received frequency. This technique allows small 
frequency changes to be observed. SFD is observed in the high-frequency 
range, and the reference signal is obtained locally from a stable oscil- 
lator. Occurrence times and deviations in frequency are listed in data 
presentations . 
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SFE - Solar Flare Effect, or crochet, is observed as a small hook 
on magnetometer records. It is caused by the magnetic field response 
to increased current flow in the E region due to electron enhancement 
induced by flare X rays. 

SPA - Sudden Phase Anomaly is observed in the same manner as SFD, 
except that VLF frequencies are used and the reference signal is the 
ground wave. For observing SPA on transmitters more distant than 
possible for ground wave propagation (about 200 km) , the ground wave 
reference signal is provided by telephone, or a local reference signal 
source is used as in SFD observations . 

SWF - Short Wave Fadeout is observed from signal strength records 
of any shortwave (3 to 30 MHz) receiver. Signals from sweep -frequency 
ionosondes (vertical or oblique incidence) may be completely absorbed. 
SWF is categorized as gradual (G) , slow (SI), or sudden (S) , depending 
primarily on how rapidly the signal loss occurs. 


CNA 


In addition to SCNA discussed under SID, there are two other types 
of CNA. Data on these latter types are not generally available but are 
mentioned here to broaden perspective. Cosmic Noise Absorption is 
observed on riometers that operate between 15 and 60 MHz, but most often 
near 25 MHz. Absorption is normalized to low nighttime absorption oc- 
currences for that station. CNA can be divided into SCA (sudden com- 
mencement absorption), CNA Type I, also called SCNA (sudden cosmic 
noise absorption), and CNA Type II. CNA Type II is local nighttime 
absorption associated with active aurora. The absorption may persist 
for hours with irregular variation including individual peaks as high 
as 8 to 10 db persisting for several minutes. This is sometimes called 
auroral absorption, but is not the same as AZA, which is not a local 
phenomenon . 

SCA (Sudden Commencement Absorption) is characterized by sudden 
onset, occurrence in zones centered on the auroral zone (58° to 75° 
geomagnetic) , and persistence for less than 1 hour. The SCA is similar 
to SCNA in rise time and duration. Types I and II CNA are somewhat 
local and irregularly distributed in longitude in contrast to the zonal 
occurrence of SCA. SCAs have been associated with geomagnetic storm 
sudden commencement X-ray fluxes . 
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VLF EMISSIONS, WHISTLERS 


Natural VLF emissions are recognized as structured noises in the 
0.2- to 40-kHz radio band (Helliwell, 1965). They are normally observed 
by use of a long antenna connected to an audio amplifier. The signals 
are usually tape-recorded and studied either by ear, using a speaker or 
earphones, or visually by use of a spectrogram. The aural data -consist 
of multi -frequency or single-frequency noises, whistles, chirps, hiss, 
chorus, and many combinations of the structured noises which are pos- 
sible. The spectrogram graphically depicts frequency on an ordinate 
versus time on the abscissa, with signal intensity indicated qualita- 
tively by the relative brightness of the plotted signals . Spectrograms 
are normally available on 4- x 12-1/2-inch sheets of facsimile paper 
(Kay Electric Sonograph) or rolls of 35-mm paper or film (Raytheon 
Rayspan) . Data summaries have been prepared from some observations 
which consist of counts of selected noises (e.g., whistlers) by time 
of day and day, based on one short observation repeated hourly. 

Most frequently, the noises are studied to determine their origin 
and their propagation characteristics. Such study is of value in under- 
standing the structure of the propagating medium, along with its tem- 
poral and spatial variations. It has been determined that whistlers 
result from frequency dispersion of lightning-generated emissions. The 
origin of other noises has not been so specifically identified, but 
these are generally attributed to origins located in ionospheric and 
magnetospheric regions. Propagation is strongly influenced by the 
magnetic field in the ionosphere and magnetosphere. 

Since 1957, a VLF technique has been used to observe propagation 
characteristics (transmission loss and group delay) of pulsed, fixed- 
frequency emissions (Helliwell, 1965). With this technique, the strength 
of the received signal versus time is observed. A closely related 
technique is to transmit a frequency modulated pattern instead of a 
pulse and to compare the frequency spectrum of the transmitted and 
received wave. Both methods require a man-made transmission of known 
characteristics and a receiver at a considerable distance from the 
transmitter. 

Time coverage of VLF wideband observations is iiregiilar and very 
limited from 1951 to 1957. Subsequent to 1957, a limited volume of 
VLF data (spectrograms and digital summaries on hard copy or microfilm) 
are available from the WDC-A for l^per Atmosphere Geophysics and are 
indexed in its catalog. More complete data for both wideband and 
fixed-frequency observations at other times may be available at the 
"observing station" or at the sponsoring institution (e.g., Stanford 
University has sponsored Antarctic observations). A list of VLF 
groimd transmitters has been prepared by Watt (1967) . 
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ATMOSPHERIC RADIO NOISE 


Atmospheric radio noise is the term applied to electromagnetic 
disturbances at radio frequencies, which occur naturally in the atmo- 
sphere. Lightning strokes associated with thunderstorms represent the 
major source of this noise below the ionospheric electron density 
maximum. 

Discussions on atmospheric radio noise can be found in Watt (1967) , 
CCIR Report 322 (1964) , and Hartz (1969) . 


Ground-Based Observations 


During 15 -minute intervals each hour, and at each of a series of 
frequencies between 0.013 and 20.0 MHz, a given station measures the 
effective omnidirectional antenna noise in decibels above a reference 
value (kTpb) . (See the 54-member station list in the WDC-A data 
catalog.) Here, k is Boltzmann’s constant (1.38 x 10"^® ergs/deg). 

Tq is a reference temperature (288°K) , and b is the bandwidth in Hz. 
The average over the 15 -minute period of observation is recorded as 
the hourly value. At the end of each month, 24 monthly median values 
of noise power, one for each hour, are determined. Upper and lower 
decile values of noise power are also determined for each hour. 

Hard-copy tables of monthly values determined by individual sta- 
tions are available at WDC-A for Upper Atmosphere Geophysics, Boulder, 
Colorado. Parameters presented include the monthly median values and 
the ratios of the upper and lower decile values to the median noise 
power in decibels . 


Satellite-Based Observations 


Observations of radio noise between .1 and 20 MHz have been made 
by the ionospheric satellites Alouette 1, Alouette 2, and ISIS 1. 

The noise observations consist of the AGC voltages measured by the 
receiver of the sounding equipment and represent noises originating 
in the ionosphere and above, which cannot be ordinarily observed at 
the ground due to the shielding effect of the ionosphere. These data 
are presently available only through the experimenter. 
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INTRODUCTION 


This section presents a discussion of the instrumentation and 
techniques used for measuring upper level temperature, winds, pressure, 
and ozone concentrations and identifies the availability and sources 
of data obtained from these measurements. The availability of satellite 
cloud photographs is also discussed. The last section of Part VI de- 
scribes some of the available global atmospheric models and the assump- 
tions upon which they are based. Specific models and their sources are 
given for both the lower atmosphere and the upper atmosphere. 


UPPER LEVEL MEASUREMENTS OF TEMPERATURE, WINDS, AND PRESSURE 


Observations of temperature, winds, and pressure are routinely 
taken in the first 250,000 feet of the atmosphere. 

Temperature is a measure of the molecular kinetic energy of the 
air. Near the surface of the earth, it is measured by a thermometer; 
vq)per air temperatures are measured by balloon-borne instruments 
(radiosondes) and rocket-bome instruments (rockets ondes) . The temper- 
ature measuring instrument on a radiosonde is a thermistor, which is 
a device having an electrical resistance that varies markedly with 
tenperature. Rocketsonde temperature instrumentation consists of a 
modified 1680-MHz radiosonde transmitter unit and a bead thermistor 
sensor. 

Wind is air in motion relative to the surface of the earth. Since 
the vertical component of wind is usually relatively small, the term 
refers to the horizontal component unless specified otherwise. Surface 
wind speed or force is measured by anemometers, and the direction by 
wind vanes. Upper level winds are measured by tracking pilot balloons, 
balloon-borne rawin or rawinsondes,* and rocketsondes . These upper 
level wind measuring devices are tracked visually with radio direction 
finders, or by radar. 


*Rawin -- A method of determining upper level wind speeds and direc- 
tions by tracking a balloon-borne radar target, responder, or radio- 
sonde transmitter with either radar or a radio direction finder. 
Rawinsonde -- A method of determining upper level temperature, pres- 
sure, relative humidity, and wind speed and direction by tracking a 
radiosonde with a radar or radio direction finder. 
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Atmospheric pressure is the pressure exerted by the atmosphere 
as a consequence of gravitational attraction exerted upon the column 
of air lying directly above the point in question. It is measured by 
many varieties of barometers and is expressed in several unit systems. 
The most common unit used is the millibar (1 millibar equals 1000 dynes 
per cm^) . 

Balloon-borne instruments are capable of making observations only 
in the lower 100,000 feet (30 km) of the atmosphere. Meteorological 
rockets routinely take observations in the lower 250,000 feet (76 km), 
especially in that portion inaccessible to balloons. See reference 
16 for a complete description of meteorological rocket instrumentation. 


Meteorological Rocket Network Data 

The Meteorological Rocket Network (MRN) was formed in October 1959 
to combine the efforts of a number of organizations engaged in upper 
atmospheric studies. The data are presented in the series: "Data 
Report of the Meteorological Rocket Network Firings," Document 109-62, 
of the Inter-Range Instrumentation Group (IRIG) of the Range Commands 
Council and later the Meteorological Working Group (MWG) of the IRIG. 

The Meteorological Rocket Network data reports were first published in 
the fall of 1959. These data reports began on a quarterly basis but 
were changed to a monthly publication with the September 1962 data. 

Each report contains the wind and temperature measurements from all 
available network firings for that particular month, along with com- 
puted values for pressure, density, and speed of sound. It also in- 
cludes the rawinsonde observations nearest the time (within 6 hours) of 
rocket firing. All stations attempt to launch rockets concurrently and 
at regular intervals. In 1959, data from only five stations were in- 
cluded in the data reports, but by 1966, data from a greater number of 
stations were included (see table 11). These data are available from 
October 1959 through February 1966. Copies of any of the 54 volumes of 
the "Data Report of the Meteorological Rocket Network Firings" can be 
obtained by writing to the Defense Documentation Center, Cameron Station, 
Alexandria, Virginia 22314. 

Since 1966, the "Data Report of Meteorological Rocket Network Firings" 
has been issued monthly by the World Data Center A (WDC-A) for Meteorology 
at the National Climatic Center (formerly National Weather Records Center) , 
Asheville, North Carolina 28801. These reports are basically the same 
as the previously issued IRIG data reports, but they also include data 
from other countries. The first issue of the WDC-A data report began 
with January 1964 data. Thus the 1964 to 1966 issues of the Meteorolog- 
ical Rocket Network data reports are essentially a reprint of the data 
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TABLE 11 


METEOROLOGICAL ROCKET NETWORK STATIONS 


Station 

Lat 

Long. 

Start 

Date 

Stop 

Date 

Antigua AAFB, BWI 

17“09'N 

61°47'W 

May 

1963 



Ascension Island AFB 

07°59'S 

14°25'W 

Aug. 

1962 



Barking Sands (Kauai) , Hawaii 

22°03'N 

159°47'W 

Oct . 

1960 



Cape Kennedy, Florida 

28°27'N 

80°32*W 

Mar. 

1960 



Eglin AFB, Florida 

30°23'N 

86°42'W 

Sept. 

1960 



Eniwetok, Marshall Islands 

11°26'N 

162°20'E 

Mar. 

1961 

June 

1969 

Fort Churchill, Canada 

58®44'N 

93®49'W 

Oct. 

1959 



Fort Greely, Alaska 

64“00'N 

145°44'W 

Dec. 

1959 



. Fort Sherman, Canal Zone 

09°20'N 

79°59'W 

Jan. 

1964 



Grand Turk Island AAFB, BWI 

21°27'N 

71°09'W 

Sept. 

1963 

Dec. 

1966 

Green River, Utah 

38®56'N 

110°04'W 

Jan. 

1964 



Harp, Seawell, West Indies 

. 13"06'N 

59°37'W 

Jan. 

1964 



Holloman AFB, New Mexico 

32“51'N 

106°06'W 

Nov .- 

1959 

June 

1961 

Kindley AFB, Bermuda 

32°21'N 

64°39'W 

Dec. 

1961 



Kwajalein, Marshall Islands 

08°44'N 

167°44'E 

Oct . 

1963 



McMurdo Sound, Antarctica 

77'’51'N 

166°39'E 

Mar. 

1962 

Oct. 

1963 

Point Barrow, Alaska 

71°21'N 

156°59'W 

July 

1961 

Aug. 

1961 

Point Mugu, California 

34°07'N 

119°07'W 

Oct . 

1959 



Primrose Lake, Alberta, Canada 

54°45'N 

110°03'W 

Jan. 

1964 



San Salvador Island AFB 

24°04'N 

74®31'W 

Sept . 

1963 

Dec. 

1964 

Thule AFB, Greenland 

76"33'N 

68°49’W 

Jan. 

1964 



Tonopah Range , Nevada 

38°00'N 

116°30'W 

Mar. 

1960 



Wallops Island, Virginia 

37°50*N 

75°29'W 

Nov. 

1959 



West Geirinish, Scotland 
White Sands Missile Range, 

57“21'N 

07°22'W 





New Mexico 

32°23'N 

106°29'W 

Oct . 

1959 



Yuma Proving Ground, Arizona 

32“52'N 

114°19'W 

Jan. 

1964 
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report of the IRIG Document 109-62, with a few more rocket stations 
reporting. Meteorological rocket data from January 1964 to the present 
time are available at the WDC-A for Meteorology. These data are on 
magnetic tape and hard copy. Subscriptions to the WDC-A "Data Report 
of the Meteorological Rocket Network Firings" should be requested from 
the Superintendent of Documents, Government Printing Office. Washington 
D.C. 20402. ® ’ 

Beginning with 1969 data, the data publication entitled "Meteoro- 
logical Rocketsonde Network" has been changed to "High Altitude Meteoro- 
logical Data." 


EXAMETNET Data 


Upper level wind data from rocketsondes are also available for the 
period 1966 to 1968 from the Experimental InterAmerican Meteorological 
Rocket Network (EXAMETNET) . This is a cooperative program among the 
national space organizations of Argentina, Brazil, and the U.S. The 
EXAMETNET meteorological rocket launchings and data dissemination are 
conducted synoptically from launch sites at Chamical, Argentina (Lat 
30"22' S, Long. 66°17' W) , Natal, Brazil (Lat 05°55' S, Long. 35°10' W) , 
and Wallops Island, Virginia (Lat 37“51' N, Long. 75°29' W) . These 
data are available in the EXAMETNET Data Report Series, Annual Report 
for 1966 (NASA SP-175), 1967 (NASA SP-176) , and 1968 (NASA SP-231), 
respectively. The reports can be obtained from the National Technical 
Information Service, Springfield, Virginia 22151, The 1969 EXAMETNET 
data are to be included in the new World Data Center A "High Altitude 
Meteorological Data" publication. 


NOAA Stratosphere Data 

A series of weekly constant-pressure charts for the upper strato- 
sphere, based on rocketsonde and very high level rawinsonde data ob- 
tained over North America and adjacent ocean areas, is prepared. These 
high level data consist of a series of charts for 5, 2, and 0.4 mb 
(corresponding to approximately 36, 42, and 55 km, respectively). Wind, 
temperature, and pressure values are included. These charts are avail- 
able for the years 1964 through 1967 as Environmental Science Services 
Administration (ESSA) Technical Reports WB-2, 3, 9, and 12 (see refer- 
ences 19, 20, 21, and 22). Copies of these documents can be obtained by 
writing to the Superintendent of Documents, U.S. Government Printing 
Office. 
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upper Air Network Data 

Wind, temperature, relative humidity, and pressure data measured by 
rawinsondes (or rocketsondes from which wind data are not available) are 
obtained from the Upper Air Network Stations, which are reasonably well 
distributed over the land masses of the Eastern and Western Hemispheres. 
The soundings are usually taken twice each day at 0000 GMT and 1200 GMT. 
These data are available at National Climatic Center, Asheville. Data 
listed daily by station are available on microfiche. The data are 
tabulated by average monthly values year by year. Synoptic charts of 
plotted data are available for sea level and the 850-, 700-, 500-, 

300-, 200-, 100-, 50-, and 30-mb levels. 

Wind data are also obtained from pibals (pilot balloons). A pibal 
is a small balloon whose ascent is followed by a theodolite in order 
to obtain data for the computation of the speed and direction of winds 
in the upper air. The pibals are released as many as four times each 
day but usually at 0600 GMT and 1800 GMT, which are off -times for the 
soundings. Most of the Upper Air Network Stations release pibals as do 
many of the principal surface weather observing stations. These data 
are also available at the National Climatic Center. 


SATELLITE CLOUD PHOTOGRAPHS 


All usable cloud photographs taken by meteorological satellites are 
available from the National Climatic Center, Asheville. 

A complete listing of the TIROS and ESSA cloud photographs are con- 
tained in the Key to Meteorological Records Documentation (KMRD) series, 
which is made up of catalogs of meteorological satellite television 
cloud photography. This series has been established to provide ^idance 
to personnel making use of the cloud photographs. It also describes how 
copies of these photographs may be obtained from the National Climatic 
Center. The catalogs are published by the National Oceanic and Atmo- 
spheric Administration and are available from the Superintendent of 
Documents, U.S. Government Printing Office. The first of these catalogs 
is the "Catalog of Meteorological Satellite Data - Tiros 1 Television 
Cloud Photography," KMRD Series No. 5.31. (See table 12.) 
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TABLE 12 


CATALOGS OF METEOROLOGICAL SATELLITE DATA 
TELEVISION CLOUD PHOTOGRAPHY 


5.31 

TIROS 

I 



1961 

Apr. 

1, 

1960 

- 

June 

15, 

1960 

5.32 

TIROS 

2 



1963 

Nov. 

23, 

1960 

- 

Sept. 

27, 

1961 

5.33 

TIROS 

3 



1962 

July 

12, 

1961 

- 

Jan. 

23, 

1962 

5.34 

TIROS 

4 



1963 

Feb. 

8, 

1962 

- 

June 

18, 

1962 

5.35 

TIROS 

5 



1964 

June 

19, 

1962 

- 

May 

14, 

1963 

5.36 

TIROS 

6 



1964 

Sept. 

18, 

1962 

- 

Oct. 

21, 

1963 

5.37 

TIROS 

7 



1965 

June 

19, 

1963 

- 

Dec. 

31, 

1963 

5.37 

TIROS 

7 



1965 

Jan. 

1, 

1964 

- 

June 

30, 

1964 

5.37 

TIROS 

7 



1965 

July 

1, 

1964 

- 

Dec. 

30, 

1964 

5.37 

TIROS 

7 



1966 

Jan. 

1, 

1965 

- 

Dec. 

31, 

1965 

5.38 

TIROS 

8 



1965 

Dec. 

21, 

1963 

- 

June 

30, 

1964 

5.38 

TIROS 

8 



1965 

July 

1, 

1964 

- 

Dec. 

31, 

1964 

5.38 

TIROS 

8 



1966 

Jan. 

1, 

1965 

- 

Aug. 

31, 

1965 

5.39 

TIROS 

9 



1966 

Jan. 

22, 

1965 

- 

Apr. 

30, 

1965 

5.39 

TIROS 

9 



1967 

May 

1, 

1965 

- 

July 

26, 

1965 

5.310 

TIROS 

10 


1967 

July 

2, 

1965 

- 

Sept. 

30, 

1965 

5.311 

ESSA 

1 



1966 

Feb. 

3, 

1966 

- 

Mar. 

31, 

1966 

5.311 

ESSA 

1 



1968 

Apr. 

1, 

1966 

- 

Oct. 

6, 

1966 

5.313 

ESSA 

3 



1968 

Oct. 

4, 

1966 

- 

Jan. 

1, 

1967 

5.313 

ESSA 

3 



1967 

Jan. 

1, 

1967 

- 

Mar. 

31, 

1967 

5.314 

ESSA 

3 

8 

5 

1968 

Apr. 

1, 

1967 

- 

June 

30, 

1967 

5.315 

ESSA 

3 

§ 

5 

1968 

July 

1, 

1967 

- 

Sept. 

30, 

1967 

5.316 

ESSA 

3 


5 

1969 

Oct. 

1, 

1967 


Dec. 

31, 

1967 

5.317 

ESSA 

3 

§ 

5 

1969 

Jan. 

1, 

1968 

- 

Mar. 

31, 

1968 

5.318 

ESSA 

3 

8 

5 

1970 

Apr. 

1, 

1968 

- 

June 

30, 

1968 

5.319 

ESSA 

3, 

C 

5, 8 7 

1970 

July 

1, 

1968 

- 

Sept. 

30, 

1968 
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The Nimbus television cloud photographs are listed and described 
in the follovdng NASA/GSFC documents: 


"Nimbus I Users' Catalog: AVCS and APT, August 28 through 

September 22, 1964" 


"The Nimbus II 
1966" 

"The Nimbus II 
1966" 

"The Nimbus II 
August, 1966 

"The Nimbus II 
30 September 


Data Catalog, Vol. 1, 15 May through 30 June 


Data Catalog, Vol. 2, 1 July through 31 July 


Data Catalog, Vol. 3, 1 August through 31 
(Orbits 1035 - 1447)" 

Data Catalog, Vol. 4, 1 September through 
, 1966 (Orbits 1448 - 1846) " 


> 


"The Nimbus II Data Catalog, Vol. 5, 1 October through 15 
November, 1966 (Orbits 1847 - 2458)" 

"The Nimbus III Data Catalog, Vol. 1, Part 1, 14 April through 
31 May, 1969 (Orbits 109 - 639)" 

"The Nimbus III Data Catalog, Vol. 2, June 1969 (Orbits 
640 - 1041)" 

"The Nimbus III Data Catalog, Vol. 3, July 1969 (Orbits 
1042 - 1457)" 

"The Nimbus III Data Catalog, Vol. 4, August 1969 (Orbits 
1458 - 1872)" 

"The Nimbus III Data Catalog, Vol. 5, 1 September to 
31 December, 1969 (Orbits 1873 - 3508)" 


Copies of these catalogs can be obtained from NSSDC. Additional Nimbus 
III and Nimbus IV data catalogs are currently being compiled and, when 
printed, will also be available from NSSDC. 


Cloud photographs are also available from the Applications Tech- 
nology Satellites (ATS) spin scan camera system. These data are listed 
and 'described in the following volumes of the "Meteorological Data 
Catalog for the Applications Technology Satellites": 

Vol. I - "Users' Guide, ATS I, January 1 through June 30, 1967" 

Vol. II - "ATS I Data Catalog, ATS II Summary, and ATS III Users' 
Guide and Data Catalog, November 7, 1967 through January 31, 1968" 
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Vol. Ill - "ATS I Data Catalog and ATS III Data Catalog, February 
1 through December 31, 1968" 

Vol. IV - "ATS I Data Catalog and ATS III Data Catalog, January 1 
through July 31, 1969" 

Vol. V - "ATS I Data Catalog, August 1, 1969 through May 25, 1970" 

Copies of these catalogs can be obtained from the Nimbus -ATS Data Utili- 
zation Center, Code 450, GSFC, Greenbelt, Maryland 20771. 

When the desired cloud photographs have been identified from these 
catalogs, they may be obtained from the National Climatic Center, 
Asheville, North Carolina 28801. 


OZONE DATA 


Ozone, O 3 , a nearly colorless gaseous form of oxygen, is found in 
trace quantities in the earth's atmosphere with a maximum concentration 
occurring between 20 and 30 km above the surface. Ozone's absorption 
of solar ultraviolet light makes possible the presence of life as we 
know it on earth. Extensive discussions of ozone and ozone measurement 
techniques can be found in Vasay (1965) . 

Parameters measured include total ozone content and vertical dis- 
tribution versus location and time. Total ozone content is the thickness, 
reduced to standard pressure and temperature conditions, of the ozone 
contained in a vertical cylindrical column of unit cross section above 
a specific location. It has been recently expressed in units of milli- 
centimeters . Ozone concentration for a given pressure level can be 
expressed as the ozone partial pressure and is in units of micromilli- 
bars. 


Most ozone observations are taken by balloon-borne ozonesondes. 
The data are commonly presented in a plotted form called ozonagrams on 
which ozone content is given versus altitude (figure 24) . Isopleths 
of ozone mixing ratios (micrograms of ozone per gram of air) can also 
be drawn on the diagram. 

Ozonagrams can be obtained from a number of sources. Data from a 
network of six stations using chemiluminescent -type ozonesondes taken 
from January 1963 to January 1966 have been published in "Ozonesonde 
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TABLE 13 


OZONE NETWORK 



Station 

Lat (“N) 

Long. ("W) 


Albrook Field, Canal Zone 

9.0 

79.6 

Cape Kennedy, Florida 

28.2 

80.6 

Goose Bay, Labrador (Canadian Meteorological 
Branch) 

53.3 

60.4 

Grand Turk AAFB, BWI 

21.5 

71.1 

L. G. Hanscom Field, Bedford, Massachusetts 

42.5 

71.3 

Wallops Island, Virginia 

37.9 

75.5 


Observations over North America," AFCRL 64-30, Volumes I to IV. (See 
Hering, 1964, and Bering and Borden, 1964, 1965, and 1967.) The par- 
ticipating stations are listed in table 13. At the beginning of each 
volume, a listing is given of the days that ozonesondes were flown at 
the various stations. 

Ozonagrams from 1962 to 1966 are available for a number of stations 
operated or supported by NOAA in Boulder, Colorado, in the Western 
Hemisphere and in Antarctica. (See Komhyr and Sticksel, 1967, and 
Dutsch, 1966.) During this period a major effort was made to gather 
ozone data during the time of the International Years of the Quiet Sun, 
1964 to 1965. These Western Hemisphere ozonagrams are published in 
Komhyr and Sticksel (1967). Antarctic data appear in Komhyr and 
Grass (1968). The stations are listed in table 14. Three different 
types of ozonesondes were used during this period: Regener chemilumi- 

nescent, Brewer-M5st electrochemical, and carbon-iodine instruments. 
Ozonagrams after 1966 are not available from NOAA at this time. 

The data from the Regener chemiluminescent soundings for the IQSY 
were extracted at 1 -minute intervals and reduced by computer. These 
data include ozone partial pressure at specified levels and total ozone 
content, and can be obtained from the National Climatic Center, Asheville. 
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TABLE 14 


NOAA OZONESONDE NETWORK 


Station 

Lat. (“N) 

Long. (°W) 

Albrook Field, Canal Zone (AWS) 

9.0 

79.6 

Colorado State diversity. Fort Collins 

40.6 

105.1 

Fairbanks, Alaska (USWB) 

64.3 

147.9 

Florida State University, Tallahassee 

30.4 

84.3 

Fort Churchill, Manitoba (Canad. Met. Br.) 

58.8 

94.1 

Goose Bay, Labrador (Canad. Met. Br.) 

53.3 

60.4 

Grand Turk AFB, BWI 

21.5 

71.1 

Green Bay, Wisconsin (USWB) 

44.5 

88.1 

L. G. Hanscom Field, Bedford, Massachusetts 

42.5 

71.3 

Point Mugu, California 

34.1 

119.1 

Thule AFB, Greenland (AWS) 

76.5 

68.8 

University of New Mexico, Albuquerque 

35.0 

106.6 

University of Washington, Seattle 

47.4 

122.3 

University of Wisconsin 

43.1 

89.4 


Two years of ozonagrams taken over Boulder, Colorado, are available 
in Dutsch (1966). Data coverage is from August 1963 to August 1965. 
Ozone measurements are presented in several forms other than ozonagrams: 
tables of ozone concentration at standard levels for single flights in 
monthly coti 5 >ilation, tables of monthly mean values of ozone concentra- 
tion at standard levels for each year, and monthly mean values of ozone 
concentration at standard levels compiled from the 2-year period. These 
measurements are taken over Boulder, and all are available in Dutsch 
(1966). 

Other statistical summaries of ozone data are available. The 
3-year 1963 to 1966 average vertical ozone distribution for each season 
for each station listed in table 15 is given in Bering and Borden (1967). 
Average values of ozone density, partial pressure, and mixing ratio are 
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tabulated for 2-km intervals. The number of observations and average 
temperature and pressure values are shown for each season and altitude 
Supplementary ozone data (coded values of total amount of ozone tem- 

levelrfn applicable to standard and significant pressure 

levels for all soundings) appear in Ozone Data for the World, 

bimonthly by the Meteorological Branch, Canadian Department o/lrans- 

available f?om 1965 ?hrS 
National Climatic Center, Asheville. Similar Japanese 

L of Japan, published monthly 

by the Jap^ Meteorologi^i Agency, Tokyo. Diti“are from 1947, 1948 

from^?h^ ®lso available 

from the National Climatic Center in Asheville. 


TABLE 15 

LOCATION OF STATIONS 
PRESENTING SEASONAL AVERAGES 


Station 

Lat (“) 

Long. (°) 

WMO* Region 

Amundsen -Scott , Antarctica 

90.0 S 

24.8 S 


Byrd, Antarctica 

80.0 S 

119.5 W 


Canton Island, South Pacific 

2.8 S 

171. .7 W 

V 

USNS Eltanin, South Pacific 




and South Atlantic Oceans 

Variable 

Variable 

III and IV 

Fairbanks, Alaska 

64.8 N 

147.9 W 

IV 

Hallett, Antarctica 

72.3 S 

170.2 E 


Hilo, Hawaii 

19.7 N 

155.1 W 

V 

La Paz, Bolivia 

16.5 S 

68.2 W 

III 

Puerto. Montt, Chile 

41.5 S 

73.0 W 

III 

Sterling, Virginia 

39.0’ S 

77.5 W 

IV 

Wilkes, Antarctica 

66.3 S 

.110.. 6 E 



*WM0: World Meteorological Organization 
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MODELS OF THE EARTH'S NEUTRAL ATMOSPHERE 


Two distinct regimes can be identified in the neutral atmosphere. 
The lower atmosphere, or the homosphere, is that region characterized 
by mixing of the components; chemical composition is nearly constant. 
The heterosphere overlaps the upper portion of the homosphere and is 
characterized by diffusive equilibrium in which chemical composition 
depends on height. The transition altitude is between 90 and 120 km. 
Above about 90 km, dissociation of molecular oxygen due to extreme 
ultraviolet (EUV) radiation is significant. 

Construction of model atmospheres is complicated by the temporal 
variability of the atmosphere. The homosphere exhibits diurnal and 
seasonal variations. The characteristics of the upper atmosphere (the 
heterosphere) change with the solar cycle, daily solar activity, geo- 
magnetic activity, local time, and geographic latitude. Seasonal and 
semiannual variations are also exhibited. 

It has proven convenient in the past to generate separate atmo- 
spheric models for the regions below and above the homosphere-hetero- 
sphere transition layer. They are likewise considered separately here. 


The Lower Atmosphere 

One of the earlier models for the lower atmosphere was the "U.S. 
Standard Atmosphere, 1962," adopted by the United States Committee on 
Extension to the Standard Atmosphere (COESA) . This model contains the 
International Civil Aviation Organization (ICAO) standard atmosphere up 
to 20 km, the proposed extension iq) to 32 km, and tables and data up to 
700 km. The tables depict an idealized, moisture-free, middle-latitude 
atmosphere. Seasonal and solar activity variations are neglected. 

"U.S. Standard Atmosphere Supplements, 1966" was prepared by COESA 
depicting variations of the models listed above. The document provides 
tabulations of atmospheric parameters from the surface up to 1000 km, 
based upon latitude and seasonal variations for the lower atmosphere. 
Both publications were sponsored jointly by the Environmental Science 
Services Administration (now called National Oceanic and Atmospheric 
Administration), the National Aeronautics and Space Administration, 
and the U. S. Air Force. They can be obtained from the Superintendent 
of Documents, U.S. Government Printing Office. The Handbook of Geo- 
physics and Space Environments (Valley, 1965) discusses these publi- 
cations and contains tables reprinted from them. 
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Another early model was represented in the COSPAR International 
Reference Atmosphere 1961 (CIRA 1961) (see Kallmann -Bijl et al., 1961). 
The tables appearing in CIRA 1961 depict average conditions represented 
by single profiles of temperature, pressure, and density for the 20- 
to 100-km region. These tables were revised in CIRA 1965 to account 
for latitude and seasonal variations. 

One of the latest models of the lower atmosphere appears in a 
publication by Groves (1970), which was prepared as part of a COSPAR 
panel project. The publication contains tables that present atmospheric 
temperature, pressure, and density from 25 to 110 km and their seasonal 
and latitudinal variations. Some of the problems of atmospheric models 
are discussed, and many comparisons of the model with actual data are 
presented. 


The Upper Atmosphere 

Two excellent review papers, those of Moe (1969) and of Blum et al. 
(1970), discuss the physics of the upper atmosphere. They also review 
the various atmospheric models that have been developed in the past and 
discuss the limitations or shortcomings of each. An extensive bibliog- 
raphy is included with each paper. 

Some of the most significant early references are Nicolet (1961), 
"U.S. Standard Atmosphere, 1962," Harris and Priester (1962a), Jacchia 
(1965), CIRA 1965, and the "U.S. Standard Atmosphere Supplement, 1966," 
which incorporated the Jacchia (1965) model atmosphere. 

Weidner et al. (1969) present a brief review paper that discusses 
the models of Harris and Priester (1962b) and of Jacchia (1965) and some 
of the limitations of each. Weidner also presents a static diffusion 
model that is a computerized version of Jacchia's model developed at the 
NASA Marshall Space Flight Center. The model varies with solar condi- 
tion and location and can be predicted for particular times and loca- 
tions. The model provides atmospheric density, chemical composition, 
temperature, molecular mass, and density scale height between 120 and 
1000 km. 

Jacchia (1970) has submitted to COSPAR (Leningrad, May 1970) a 
provisional model atmosphere for altitude range 90 to 2500 km. This 
model, with modifications in some numerical values and for altitudes 
above 110 km, will appear as the ( CIRA 1971) model. In constructing 
this model, Jacchia assumes constant temperature and density values 
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at 90 km. He also chooses a temperature profile T = T(h; T^) , which 
yields results compatible with the density data of satellite drag 
analysis. In the equation, h is height and is the temperature of 
the isothermal exosphere.* 

To obtain total mass density values, Jacchia integrates the equation 
of hydrostatic equilibrium up to 105 km, the assumed upper limit of the 
mixing regime. The mean molecular mass variation, which is non-zero 
due to oxygen dissociation, is described by a sixth degree polynomial. 

To obtain densities as functions of height in the region of assumed 
diffusive equilibrium above 105 km, Jacchia integrates the diffusion 
equation, which is the equation of hydrostatic equilibrium for each 
species, with an additional thermal diffusion term used for helium. 
Special treatment is given to the hydrogen component. 

The resultant tables of atmospheric parameters are presented for 
exospheric temperatures (T^) between 600“K and 2000°K, at 100°K inter- 
vals. Temperature, mass density, mean molecular weight, pressure scale 
height, and number densities of molecular nitrogen and oxygen and of 
atomic hydrogen (above 500 km), helium, argon, and oxygen are given at 
95 discrete altitudes between 90 km and 2500 km. Separate tables pre- 
sent total mass densities only as functions of altitude for exospheric 
temperatures between 600°K and 2000“K, in 50®K steps. Since good 
observational data do not exist for altitudes above 1100 km, tabular 
data above this height must be considered as unconfirmed extrapolation. 

Jacchia (1970) includes equations and tables that permit the user 
of the model to choose the appropriate value of the exospheric tempera- 
ture for given values of the solar activity level (as characterized by 
the 10.7-cm flux level), the geographic latitude and local time, the 
geomagnetic activity level (characterized by Kp or Ap) , and the phase 
of the year. Seasonal and latitudinal variations of the mass density 
in the lower thermosphere and of helium are treated separately. 


*The alternative to choosing a temperature profile en5)irically is to 
assume the presence of heat input and transport mechanisms and to 
solve an energy conservation equation.. This approach has been fol- 
lowed by others. See, for example, Harris and Priester (1962a). 
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MODELS OF MAGNETOSPHERICALLY TRAPPED PARTICLES 


Due to the large-scale geometry of the geomagnetic field, energetic 
charged particles may be constrained for long periods of time with a 
combined motion consisting of gyration about field lines, latitudinal 
bounce along field lines, and azimuthal drift around the earth. The 
trapped particle populations, often referred to as the Van Allen belts, 
have been studied using data from numerous spacecraft since 1958. 

Recent discussions of these particles include those by Hess (1968) and 
by Williams (1970) . 

Trapped particle data gathered by various spacecraft have been 
analyzed, synthesized and presented in a format useful to space systems 
analysts and space scientists. The resultant model environments are 
distinguished from each other by particle species, spatial and temporal 
coverage, and energy range. 

Within a given model, the integral omnidirectional flux, J(>E, B, 
L), is given as a product of the flux, F(B, L) , above a certain energy, 
Ej, times a spectral function, N(>E, E^; B, L) . For proton models, the 
spectral function is represented over a limited energy interval by* an 
exponential function, N = N(,exp[(E-Ej)/E (B, L>] , or- by a power law, 

N = N (E/Ej)"PCB, L) , Thus, a proton moael includes specification of 
the functions F(B, L) and E(, (B, L) or P(B, L) . On the other- hand; for 
electron models, the spectral function N is given as a tabular fuhct4on 
of E and L, with no B dependence. For each environment, the daily 
accumulated fluxes above certain energy levels and within certain energy 
intervals are also given for an assumed spacecraft in various orbits. 

The exception to this discussion is the synchronous -altitude elec- 
tron environment, AES, in which the mean differential and integral 
fluxes (for solar minimum) are explicitly presented for given values 
of B/Bo, local time, and energy. Probability distributions for devia- 
tions from mean values are also given. 

The following model electron environments are currently available: 


Env. 

E, (Mev) 

Energy (Mev) 

Spatial 

Temporal 

. EPOCH 

Name 

1 

Range 

Range 

Range 

AEl 

0.5 

>0.5 

1.2<L<3.0 

1962-1963 

7/63 

AE2 

0.5 

O-oo 

1.2<L<6.2 

1962-1964 

8/64 

E68 

0.5 

O-oo 

1.2<L<6.2 

1962-1964 

1968 

AES 


0.01-5 

L=6.6 

1959-1965 

1964 
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The following model proton environments are currently available; 


Env. 

Ej(Mev) 

Energy 

Spatial 

Temporal 


Name 


Range (Mev) 

Range 

Range 

Remarks 

API 

34 

30-50 

1.2<L<2.8 

1958-1963 


AP2 

15 

15-30 

1.2<L<3.0 

1958-1963 


APS 

50 

>50 

1.2<L<2.8 

1958-1963 


AP4 

4 

4-15 

1.2<L<4.2 

1962-1963 


APS 

0.4 

0.1-4* 

1.2<L<6.6 

1961-1965 


AP6 

4 

4-30 

1.2<L<4.0 

1962-1965 

Supersedes AP2 § 4 

AP7 

50 

>50 

1.15<L<3.0 

1961-1966 

Supersedes AP3 


New inner and outer zone electron environments are currently being 
prepared. Updating of model environments will proceed as warranted by 
new data as they become available. 

Except for E68, the model environments have been published in the 
NASA SP-3024 series. Models of the Trapped Radiation Environ ment. The 
following volumes are currently available from NSSDC: 


Volume 

1: 

Inner Zone Protons and Electrons (Vette, 1967) 
(Includes model environments AEl, API, AP2. APS 
AP4) 

Volume 

2: 

Inner and Outer Zone Electrons (Vette, Lucero, and 
Wright, 1967) 

(Includes model environment AE2) 

Volume 

3: 

Electrons at Synchronous Altitudes (Vette and Lucero. 
1968) 

(Includes model environment AES) 

Volume 

4: 

Low Energy Protons (King, 1967) 
(Includes model environment APS) 

Volume 

5: 

Inner Belt Protons (Lavine and Vette, 1969) 
(Includes model environment AP6) 

Volume 

6: 

High Energy Protons (Lavine and Vette, 1970) 
(Includes model environment AP7) 


Using API, AP2 , APS, AP4, and AE2, Stassinopoulos (1970) has gener- 
ated isoflux contours for integral electron and proton fluxes . These 
contours in geographic latitude and longitude are given at discrete 
altitudes in 100-km increments between 200 and 1000 km. 

*The lower energy limit is L dependent, being 0.1 at L>3.5 and 1.0 at 
L<1.8. 
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Except for AES, NSSDC can supply a BCD tape or card decks containing 
distribution functions F(B, L) and spectral functions N(E, L) for elec- 
trons and Eo(B, L) or P(B, L) for protons. NSSDC can also supply punched- 
card copies of the TRECO program, with which a user can generate the 
orbital integrated fluxes mentioned earlier. 


SOLAR WIND 


The solar wind is the interplanetary extension of the solar corona 
in its state of dynamical expansion. Parker (1969) and the references 
contained therein provide detailed information on the solar wind. 

Using multichannel observations (in energy and in angle) of posi- 
tive ion fluxes, macroscopic bulk flow parameters such as flow velocity 
and particle density are derived. A secondary derived parameter is the 
time required for a steady-state solar corotating plasma beam to rotate 
from the spacecraft to the earth; this time is referred to as the coro- 
tation delay time. 

In Solar -Geophysical Data, the bulk velocity and corotation delay 
times as provided by both the NASA/AMES and MIT plasma groups from their 
Pioneer 6 and 7 data are given for the month prior to the publication 
date. Proton densities are also provided by the MIT group from Pioneer 
6 data. In addition, bulk velocity and proton density values are pro- 
vided by the Los Alamos plasma group from its Vela 3 and 4 data. 


AIRGLOW 


Airglow is a faint luminescence of the earth’s upper atmosphere. 

The emitted electromagnetic radiation results mainly from chemical 
reactions of upper atmospheric constituents, although ionic recombina- 
tion and bombardment by energetic particles are other mechanisms that 
can cause airglow. The main emissions include the atomic oxygen for- 
bidden green line (5577 A) and forbidden red line (6300 A) , the sodium 
doublet (5890 A and 5896 A), and the hydroxyl (OH) band system. For a 
more complete discussion of airglow, see Chamberlain (1961). More infor- 
mation on photometric observations of airglow can be found in Roach (1957). 
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A photometer is used to measure the intensity of the various lines 
and bands in the airglow spectrum. Photometric data have been obtained 
at several stations on a regular basis. The World Data Center A Cata- 
logue of Data on Solar -Terrestrial Physics lists 41 stations active 
during at least part of the time since 1959, according to the wavelength 
at which measurements were made. The data from some stations are held 
at the WDC-A for Upper Atmosphere Geophysics, and the catalog indicates 
the completeness of these data. Zenith intensities in rayleighs and 
ratios of intensities at equal angular distances (usually 75 , although 
70® and 67® are common) north and south from the zenith are the two types 
of data held at WDC-A. Observations are made for each hour on the hour 
(UT) The ratio data have advantages in that they are independent of 
the absolute calibration of the photometer, and they give the relative 
intensities at widely separated geographic locations. Most of these 
data are available in printed tabular form. However, some data obtained 
during the IGY and IQSY are available on punched cards. 

Photometric data, zenith and ratio, from 19 stations for the IQSY 
are contained in WDC-A report UAG-1, July 1968. The report, entitled 
"IQSY Night Airglow Data," was prepared by L. L. Smith, F. E. Roach, 
and J. M. McKennan. 

In some instances, the data from a particular station exist or are 
presumed to exist but are not held at WDC-A. The catalog indicates this, 
and, in these cases, WDC-A will assist in obtaining information concerning 
the completeness and availability of data from these stations. 


AURORA 


An aurora involves emitted electromagnetic radiation that results 
primarily from the interaction of energetic, extra-atmospheric particles 
with the neutral gases of the upper atmosphere. These emissions extend 
from the ultraviolet to the infrared. The portion of the emitted elec- 
tromagnetic radiation that is visible to the human eye is referred to 
as the visual aurora. For an in-depth treatment of the aurora, see 
Chamberlain (1961) . 


The World Data Center A for Upper Atmosphere Geophysics in Boulder 
is a repository for visual auroral data. Photographs obtained from 
all-sky cameras constitute the principal type of available data. These 
photographs are taken throughout all clear nights, and at least 12 
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exposures per hour are made. The photographs show the whole sky as 
reflected in a convex mirror, and all but the faintest aurora are re- 
corded. The data are then put on 16- or 35-mm film in rolls 100 feet 


All -sky camera photographs have been taken on a regular basis 
Camera stations are located mainly in the northern and southern aurcral 
regions (60 N to 90“N, 60°S to 90“S), but a few are in the subauroral 
regions (45 N to 60 N, 45 S to 60®S) . (All references are to geomag- 
netic latitude.) Locations of stations that have been involved in this 
t^e of observation for any substantial length of time during the period 
since 1957 are listed in the WDC-A Catalogue of Data oh Sblaf- Terrestri al 
^pics. In instances where data from a particular station are held at 
WDC-A, the time period and completeness of the data are indicated. The 
catalog also lists cases for which data exist or are presumed to exist 
but are not held at WDC-A. In these cases, WDC-A will assist both in 
determining the availability of data and in obtaining it. 

In addition to microfilm, all -sky camera data can be obtained in 
the form of the ascaplot. The ascaplot contains information such as 
the part of the sky where the display is located, the auroral intensity 
(weak medium, strong), and the presence of cloudiness. Information is 
recorded on the plot every half hour. Only a limited amount of data 
essentially the IGY and IQSY data, is held at WDC-A in this form. 


Visual observations form the basis of another type of data. Infor- 
mation is collected from a network of observers. This information is 
then summarized and can be put into several different forms. 


One form, which indicates the occurrence or nonoccurrence of auroras 
without any reference to form or brightness, is known as a visoplot. 

With geomagnetic latitude as ordinate and universal time as abscissa 
various symbols are used to indicate whether, at a certain location *an 
auroria is observed or is estimated to be overhead, to the north, or*to 
the- south. This is done for each of eight 45° geomagnetic longitude 
sectors. Magnetic activity index values, Kp, are given for 3-hour 
periods . Northern and Southern Hemisphere visoplots for the IGY period 

Annals of the International Geophysical 


Another form, that shows the auroral geometry and the geomagnetic 
latitude and longitude of a display at each hour of universal time is 
the synoptic auroral map on polar projection. By the use of various 
symbols, these maps show the location of veils, patches, arcs, rays 
and bands. Also given is the value of Kp for that particular ’time 
These maps are available on microfilm from WDC-A for certain time 
periods. Northern Hemisphere maps are available for 1964, 1965, 1967, 
and 1968. Southern Hemisphere maps can be obtained for the years 
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1964 and 1965. A similar form for displaying auroral data from the U.S. 
and Canada is the North American map. This type of map is not on polar 
projection, but it does indicate the locations and forms of observed ^ 
aurora. A large number of these maps are held at WDC-A, and the avail- 
able time periods can be obtained upon inquiry. Still another form in 
which some visual data are recorded is punched cards. Since 1957, 
various modifications have been made in the system that is used to 
convert data from observers' reports to cards, and, consequently, a 
wide variety of basic data cards are at WDC-A. Information recorded on 
punched cards may include time, point of observation, location of the 
display in the sky, duration, motion, forms, form sequence, and color. 

These three types of data, i.e., visoplots, synoptic maps, and 
punched cards, can be obtained from WDC-A for certain periods. The 
available time periods are listed in the WDC-A catalog. 

The all -sky camera and visual observation data are supplemented by 
several other groups of data. A limited amount of radar echo and patrol 
spectrograph data is available. Radio aurora is a term used for the 
ionization associated with the aurora and is responsible for certain 
types of radio reflections. Some radar echo data are available for the 
period 1957 to 1962.. The data are on 16-ram film and consist of pictures 
of radar range echoes taken from the face of a cathode-ray tube. A 
patrol spectrograph is a low-resolution measurement that- can be used 
for verifying the presence or absence of spectral lines and bands. Some 
patrol spectrograph data are available for the period 1957 to 1960. 
Spectrograms have been recorded on 16-mm film. The WDC-A catalog lists 
the stations that recorded these two types of data and the amount of 
data available. Most of these same patrol spectrograph data can be 
found in tabular form in volumes 25 and 40 of the Annais of the Iriter- 
national Geophysical Year. The information from the. spectrograms was 
punched on cards, and the card data were used to compile the tables in 
these volumes. In addition, both volumes, contain a description of the 
instrumentation and ‘of the data reduction procedures. 


SOLAR GEOPHYSICAL CALENDAR. RECORDS, ABBREVIATED CALENDAR RECORDS, 
AND CONDENSED- CALENDAR RECORDS 


Solar Geophysical Calendar Records contain the daily highlights of 
solar geophysical activity' for the period 1957 to 1964. They consist 
of a collection, on a daily basis', of data related, to geomagnetism, 
solar activity, ionosphere^ aurora, and cosmic rays. Charts, mainly 
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containing indices and brief textual discussions comprise the Calendar 
Records. The Calendar Records are useful in that they may readily sug- 
gest correlations between various phenomena; however, it is recommended 
that original data be used for more detailed studies. 

Calendar Records for the IGY and for the ICC (International Geo- 
physical Cooperation - 1959) are found in the Annals of the Intema- 
tional Geophysical Year , volume 16 (pp. 1-159 iHd' 201-300, ' respectively) 

Since 1964, Abbreviated Calendar Records have been generated. 

somewhat smaller quantity of tabulated data than 
e earlier Calendar Records, but, at least through June 30, 1967 they 
contain ejq)anded textual discussion of daily highlights. Since July 1, 
1967, most data have been presented in chart form, with very limited 
textual discussion. ^ 


Data currently included in Abbreviated Calendar Records are as 
follows. For each month, a calcium plage map is given; for each day, 
the times (start and duration), importance, and associated calcium 
plage region of flares greater than If, start times of solar noise 
bursts in four frequency ranges, occurrence of spectral types II or 

X-W bursts 

^ geomagnetic indices Ap and Kp, ionospheric indices I, Ip, 

and la, 10-cm solar flux value, sunspot number Rz, and central meridian 
passage information related to calcium plages are given. When appro- 
priate, textual information relates to the geomagnetic character of the 
day (quiet or disturbed), auroral displays, general ionospheric condi- 
tions , Forbush decreases, polar cap absorption, noctilucent clouds 
coronal emission, and solar proton flux increases. 


Some data contained in Abbreviated Calendar Records are based on 
provisional data reports available at WDC-A for Upper Atmosphere Geo- 
physics. Thus, It IS recommended that the Records be used as a guide 
to more definitive data that can be obtained from WDC-A or found in 
standard data publications. 


For the period J^uary 1964 to December 1966, Abbreviated Calendar 
Records can be found in selected issues of the IQSY Notes For the 
period January 1966 to December 1967, these Records 5?TTound in WDC-A 
for Upper Atmosphere Geophysics Report UAG-4. (See Data Sources, Part 
IX, tor a description of the document and procedures for acquiring it ) 
For the periods January 1967 to April 1967, May 1967 to September 1967 
November 1967 to April 1968, May 1968 to June 1968, and July 1968 to ’ 
November 1968, Abbreviated Calendar Records are found iri STP Notes, 

No. 1, 2, 3, 5, and 7, respectively. (See Data Sources, Part IX, for 
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discussion of STP Notes.) For the period December 1968 to the present, 
these Records are found in the Solar-Geophysical Data (Comprehensive 
Reports) , with a 7 -month lag. 

An even more succinct presentation of solar -terrestrial data is 
the Condensed Calendar Record, the first of which appeared for the 
period December 1968 through September 1969 in STP Notes, No. 7 (May 
1970). One line gives the following data for one day: the sunspot 

niimber Rz, the 10-cm solar flux value; the geomagnetic Ap index; the 
ionospheric I, Ip, and la indices; the number of flares of importance 
one or above; the number of X-ray and SID events; occurrence of type 
IV radio emission; a proton event index; and comments. A monthly 
summary of two to five lines is also given. Future Condensed Calendar 
Records will also be published in the STP Notes. 


SOLAR-TERRESTRIAL ACTIVITY CHARTS 


As an outgrowth of the First General Assembly of the Inter -Union 
Commission on Solar Terrestrial Physics (London, January 1969), T. 
Obayashi and others at the National Committee on Solar -Terrestrial 
Physics, Science Council of Japan (Tokyo), are generating solar -ter- 
restrial activity charts on three time scales: yearly charts, 27-day 
charts, and special event charts. Each chart displays individual 
time plots related to: solar radio and X-ray emission, solar flare 

occurrence, solar proton flux, solar wind velocity and density, inter- 
planetary magnetic field intensity and direction, neutron monitor data, 
and geomagnetism (Dst and Ap) . On the yearly charts, most parameters 
are presented as daily averages. On the 27-day charts and special 
event charts, greater time resolution is presented. A sample yearly 
chart for the year 1968 appears in STP Notes, No. 7. Information spec- 
ifying the publication in which these charts will routinely appear will 
be given in a subsequent issue of the STP Notes . 


JOVIAN RADIO EMISSIONS 


Except for the sun, the planet Jupiter is the only source of intense 
sporadic radio emissions that can be observed from the earth. These 
Jovian emissions occur at wavelengths greater than 7.5 meters (a frequency 


165 



of 40 MHz). These emissions are observed with either fixed -frequency or 
sweep-frequency receivers such as those described by Kraus (1966) . A 
conq)rehensive review of Jovian radio astronomy is given by Warwick (1967) . 
Jovian radio emission data include relative intensity, burstiness, and 
frequency range of the emissions. 

Data from the University of Colorado 7.6- to 41-MHz sweep -frequency 
receiver, covering the period January 1960 to September 8, 1968, are 
available in hard-copy form from WDC-A at Boulder. Florida State Uni- 
versity data from single -frequency or multichannel receivers operating 
between 18 and 38 MHz and covering the 1961 to 1964 period are available 
in tabular form on microfilm and as oscillograph tracings on microfilm 
at NSSDC. Similar data from the University of Florida and from a co- 
operative station in Santiago, Chile, for the period 1957 to 1962, are 
also available at NSSDC. Goddard Space Flight Center (GSFC) data taken 
at 16.7 and 22.2 MHz between November 11, 1965, and February 28, 1966, 
are available in hard-copy tabular form at NSSDC. 

Data are still being taken by the University of Florida and the 
GSFC groups and are published on an irregular basis. 
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COORDINATE SYSTEMS 


A number of coordinate systems that have proven useful in the 
analysis of space science data are discussed in this part. Some co- 
ordinate systems are better suited for certain studies than are others . 
The various systems are distinguished by the orientation of their axes. 
Except for the B, L and a, B coordinate systems, all systems are right 
handed and orthogonal; as such, the directions of the Y axes are not 
explicitly given. Coordinate systems used principally in spacecraft 
trajectory analyses or similar studies are not included in this dis- 
cussion. 


Geographic Coordinate System 

In the geographic coordinate system, the positive Z axis is aligned 
northward along the earth's spin axis. The positive X axis lies in the 
Greenwich meridian and intersects the earth's surface in the Atlantic 
Ocean, south of Ghana, Africa. 

The geocentric latitude of a point, P, is the angle measured at 
the earth's center between the equatorial plane and a radius to P. The 
geodetic latitude of P (sometimes called the geographic latitude) is 
the angle between the equatorial plane and the line which passes through 
P and which is normal to the oblate ellipsoid whose surface most nearly 
approximates the mean sea-level surface of the earth (i.e., the geoid) . 
The maximum difference between geodetic and geocentric latitudes at 
the earth's surface is 11.6 minutes of arc and occurs near 45° latitude. 
The astronomical latitude of P is the angle between the equatorial 
plane ^d the plumb-bob vertical direction at P. The difference between 
geodetic and astronomical latitudes is almost everywhere less than 1 
minute of arc. 


Equatorial Coordinate System 

In the equatorial coordinate system, the positive Z axis is per- 
pendicular to the earth's equatorial plane (positive north), while the 
X axis points from the earth to the sun at the moment of the vernal 
equinox (or, equivalently, to the point of Aries) . This is an inertial 
system and is sometimes referred to as the geocentric celestial inertial 
coordinate system. Due to the slow precession of the earth's equatorial 
plane, an epoch date should be defined for the sake of completeness. 
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Right ascension is defined as the angular distance east of the 
vernal equinox as measured along the earth’s equator; declination is 
the angular distance north (positive) or south of the earth’s equatorial 
plane. Right ascension and declination are often used in specifying 
vector directions (as of magnetic fields) in an inertial coordinate 
system. 


Geomagnetic Dipole Coordinate System 

In the geomagnetic dipole coordinate system, the positive Z axis 
is aligned northward along the geocentric geomagnetic dipole axis. 

When the field model defined by the g^, g\, and hj terms of the IGRF 
1965.0 model is used (a common choice), this positive Z axis inter- ^ 
sects the earth’s surface at 11.44® geocentric colatitude and 290.24 
geographic east longitude (extreme western Greenland). The positive 
X axis lies in the plane formed by the geographic and geomagnetic Z 
axes and intersects the earth in southeastern Peru. Analogous to the 
conventional local time system, geomagnetic local time is defined by 
replacing the geographic axis with the geomagnetic axis and defining 
geomagnetic noon for a location P as the time when the subsolar point 
lies on the geomagnetic meridian of P (Sugiura and Heppner, 1968) . 

The transformation between the geocentric and geomagnetic coordi- 
nate systems is given in a spherical representation as follows. If 0 
and $ are geomagnetic colatitude and east longitude, respectively, 
while 0 and <)> are geocentric colatitude and east longitude, respectively, 
then 


cos 0 = cos e cos(11.44®) + sin 0 sin(11.44®) cos (4>-290.24®) 

sin $ = sin 0 sin((|»-290.44®)/sin 0 

cos ♦ = cos 0 cos 11.44® - cos (|)/(sin 0 sin 11.44®) 

A Fortran program to transform positions and field components from 
geocentric to geomagnetic dipole coordinates and vice versa has been 
generated by G. D. Mead (1970). Program decks are available at NSSDC. 


Solar Magnetic Coordinate System 

In the solar magnetic coordinate system, the positive Z axis is 
aligned northward along the geocentric geomagnetic dipole axis while 
the X axis lies in the plane formed by the Z axis and the earth-sun 
line; X is positive in the sunward hemisphere. Solar magnetic coordi- 
nates are recommended for those magnetospheric studies in which the 
main geomagnetic field is a dominant factor. 
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Solar Magnetospheric Coordinate System 


^ In the solar magnetospheric coordinate system, the positive X 
^is points toward the sun while the Z axis lies in the plane formed 
by the X axis and the geomagnetic dipole axis; Z is positive in the 
Northern Hemisphere. This coordinate system is recommended for those 
magnetospheric studies involving the geomagnetic tail and for others 
in vdiich the overall magnetosphere configuration is a dominant factor. 


Solar Ecliptic Coordinate System 

In the sol^ ecliptic coordinate system, the Z axis is perpendic- 
ular to the ecliptic plane (positive north) and the X axis is directed 
toward the sun. In terms of the equatorial coordinate system pre- 
viously discussed, the Z axis of the solar ecliptic coordinate system 
has a right ascension of -90“ and a declination of 66.56“ 


The sol^ ecliptic coordinate system is recommended for studies 
of solar-wind-related interplanetary phenomena. The solar ecliptic 
coordinate system should not be confused with the heliocentric 
ecliptic coordinate system used in some spacecraft trajectory anal- 
yses. In the latter system, the Z axis is also normal to the 
ecliptic plane, but the X axis lies along the earth-sun line defined 
at the tune of the vernal equinox. 


B, L Coordinates , Invariant Geomagnetic Coordinates 

To properly organize magnetospheric trapped particle data, 
Mcllwain (1961, 1966) introduced the B, L coordinate system. B is 
the magnetic field strength at the point of interest, and L is ex- 
plained as follows. In a dipole field, the equatorial crossing 
dist^ce, Rq, of a field line passing through a spatial point of 
interest is functionally related to the values, at that point, of B and 
I. Here, I is the integral (or longitudinal) invariant. Thus, in a 
dipole field, Rq = F(B,I). For a nondipolar field, L is defined by 
L = F(B,I). B and I are evaluated at the point of interest using the 
nondipolar field, and the function F is taken to be the same as for the 
dipole field. 

The advantage in the use of the B, L coordinates is that a parti- 
cle conserving its first two adiabatic invariants mirrors at the same 
B, L values as it drifts around the earth. Thus, in the absence of 
significant azimuthal gradients in particle populations such as those 
resulting from the South Atlantic magnetic anomaly, B, L coordinates 
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organize perpendicular flux data very well. To the extent that the 
geomagnetic field approximates a dipole field, the B, L coordinates are 
useful in organizing omnidirectional fluxes. This is the case for L 
values out to about 5. At L values greater than about 5, the geomag- 
netic field is sufficiently nondipolar that each field line is charac- 
terized by a range of L values, and the B, L coordinates are somewhat 
deficient in organizing omnidirectional fluxes. 

The INVAR program deck, generated by Mcllwain, used for calculat- 
ing B, L values at a given point in the earth's vicinity for any given 
geomagnetic field model is available from NSSDC. A second Fortran 
program, generated by Kluge (1970) computes B and L values faster 
and more accurately than does INVAR. This deck will become available 
at NSSDC shortly. A coefficient deck for a field model, such as those 
discussed in the geomagnetic field models section in Part IV of this 
Handbook , are used with these B, L programs. 

Using the Hendricks and Cain 1960.0 geomagnetic field model, as 
projected to 1965.0, Stassinopoulos (1970) has generated contours of 
constant B and L. These contours in geographic latitude and longi- 
tude are given at discrete altitudes in 100-km increments between 0 km 
and 3000 km. Similar plots, based on the GSFC (12/66) geomagnetic 
field as updated to 1970.0, are found in Wiley and Barish (1970) for 21 
altitudes between 0 km and 10,000 km. (This reference also contains 
contours of constant conjugate altitude difference, conjugate latitude 
and longitude, field dip angle and declination, and horizontal field 
intensity.) 

Useful auxiliary quantities have been defined to complement B, L 
coordinates. Variables R and X (a radius and a latitude) have been 
defined in terms of B and L by assuming dipolar relations (Roberts, 
1964). The quantity x [x=(1-Bq/B) ^/^] , where Bo is the equatorial 
value of B, is the cosine of the equatorial pitch angle of a particle 
mirroring at the point of interest (Roberts, 1965). Of the five vari- 
ables, B, L, R, X, and x, only two are independent. Vette and Porjes 
(1966) have generated 17 tables showing the relation of one of the 
five variables to two others; the Jensen and Cain (epoch 1960) model 
magnetic field was used. Copies of this document are available at 
NSSDC. 

Another pair of useful variables is the invariant geomagnetic 
latitude and longitude, A and $, respectively. A is defined simply by 
cos^A = 1/L (see Sugiura and Heppner, 1968); A is a field line earth- 
intersection latitude derived from the actual L value using a dipolar 
relation. The invariant geomagnetic longitude $ cannot be defined as 
simply as the invariant geomagnetic latitude A because the choices of 
a polar axis, of an equatorial plane, and of a prime meridian involve 
a series of arbitrary decisions. 
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Evans et al. (1969) present a system of invariant geomagnetic 
coordinates (not to be confused with the geomagnetic coordinate systems 
discussed previously) based on the International Geomagnetic Reference 
Field [IGRF 1965.0 discussed previously] as updated to 1969.75. This 
document, which can be obtained from the authors (see page 196), describes 
the difficulties encountered in constructing such a system and presents 
reasons for the decisions that lead to their definition of the invariant 
geomagnetic longitude $. 

They present, at each of six altitudes (0, 100, 300, 600, 1000, 
and 3000 km) and in each hemisphere, a polar azimuthal equidistant 
projection plot containing concentric constant-geographic-latitude (X) 
circles (40® S. ^ 5. 90®; AX = 5®) and radial constant -geographic -longi- 
tude (4>) lines (0® ^ <j> <. 360®; A4» = 30®). Superposed on these are 
lines of constant invariant geomagnetic latitude (AA = 5®) and lines of 
constant invariant geomagnetic longitude (A$ = 10®). They also present, 
at each of the six altitudes identified above, a cylindrical equidistant 
projection plot containing a Cartesian grid of constant geographic 
latitude (-90® _< X _< 90®; AX = 30®) and longitude (0° < (p < 360®; A(f> = 30®) 

lines. Superposed on these lines are lines of constant invariant geo- 
magnetic latitude (AA = 10°) and longitude (A4 = 10®) lines. Finally, 
they present for each of the same six altitudes numerical listings of 
calculated values of geographic latitude and longitude for specified 
values of invariant geomagnetic latitude (-90® ^A < 90®; AA = 5®) and 
longitude (0° <. $ <. 360°; A$ = 10°). 

NSSDC has available for full scale reproduction the master copies 
of the graphs found in Evans et al . (1969) . 


g, B Coordinate System 

The g, 8 coordinates, also called Euler Potentials, are constant 
on a given field line and satisfy B fr . t) = 0‘(£> t) x ^ 3(r, t) . 

A third coordinate, usually designated by s, provides a measure of dis- 
tance along a field line from some reference point. The g, 8 coordi- 
nates are not unique in that two g, 8 pairs related by a Jacobian of 
unity are equally valid. Some of the advantages of this system are 
indicated in Stem (1968); the relation of g, 8 to B, L is also given 
in this reference. 
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Heliocentric Corotating Coordinate System 


For studies involving the evolution of solar surface features, a 
coordinate system in which these features are stationary is advanta- 
geous. An approximation to such a system was developed by Carrington 
and is discussed in The American Ephemeris and Nautical Almanac (1968, 
p. 494). In this system, the positive Z axis points northward along 
the solar spin axis. Heliographic longitudes on the solar surface are 
measured from the solar meridian which passed through the ascending 
node of the solar equator on the ecliptic on January 1, 1854, Green- 
wich mean noon; the longitudinal angle increases in the direction of 
solar rotation. A sidereal period of rotation of 25.38 mean solar 
days is assumed for this coordinate system. Note that the differential 
character of solar rotation is neglected. 


INTERNATIONAL URSIGRAM AND WORLD DAYS SERVICE 


Data Exchange and Forecasts 


A system for the real-time collection, interchange, and distri- 
bution of solar and geophysical data, alerts, and forecasts has been 
established under the direction of the International Ursigram and World 
Days Service (lUWDS) . A thorough description of this system, along 
with a list of contributing observatories and the operating procedures 
of each of 11 Regional Warning Centers (RWCs), is given in the "Inter- 
national Ursigram and World Days Service Synoptic Codes for Solar and 
Geophysical Data." 

Briefly, individual observatories report data to the most con- 
venient RWC. These data may involve sunspots, calcium plages, solar 
flares, solar X rays, solar radio emission, coronal line emission, solar 
particle events, solar wind parameters (Pioneer satellite data), iono- 
spheric critical frequencies, sudden ionospheric disturbances, ionospheric 
absorption, geomagnetic activity, auroral observations, or cosmic rays. 
These data are exchanged among RWCs on a scheduled basis once or twice 
daily; each RWC then forwards the data to users in its region. 

When an observatory records a major event (optical or radio flare, 
cosmic-ray event, geomagnetic storm), the pertinent information is 
transmitted to the most convenient RWC, which immediately issues an 
ADALERT (ADvance ALERT) to all other RWCs. 
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The Regional Warning Center for the Western Hemisphere, located 
at the National Oceanic and Atmospheric Administration (NOAA), Boulder, 
Colorado, serves as the World Warning Agency (WWA) . Once each day, 
this RWC issues to all other RWCs a GEOALERT (GEOphysical ALERT) con- 
taining information on the commencement, continuation, or termination 
of geomagnetic storms, cosmic-ray events, and stratospheric warnings, 
as well as forecasts of flares, solar proton events, and geomagnetic 
storms. The contents of the GEOALERTs are based on the data provided 
by the worldwide network of RWCs. 

The data handled in the lUWDS-sponsored activity also go to the 
World Data Centers (WDCs) . The data catalogs from the various WDCs 
are referenced in the discipline-oriented parts of this Handbook and 
are discussed further in Part IX. This system will probably be of 
interest primarily to those data users who require essentially real- 
time data. 


International Geophysical Calendar 

Each year, a set of "world days" is selected by lUWDS and is pre- 
sented on the International Geophysical Calendar. The selections pro- 
vide various kinds of samples (3 days per month, 14 contiguous days 
about every third month, etc.) for intensified measurement of the time 
variation of geophysical phenomena. The calendar also indicates pre- 
dictable events such as meteor showers and eclipses and any special 
intervals selected for intensified cooperative observations such as 
for the Proton Flare Project. 

The International Geophysical Calendars for the period 1957 to 
1969 are contained in Shapley and Lincoln (1969) . This is a useful 
document in that for some disciplines the data held by World Data 
Centers is increased for some of the world days marked on the calendar. 
Current calendars appear annually in the WDC-A Catalogue of Data on 
Solar-Terrestrial Physics and in the January 1 issue of the Journal of 
Geophysical Research. 


ROCKET DATA 


Rocket data may be conveniently categorized into three types as 
follows. 
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Synoptic Observations 

Synoptic observations are rocket observations that involve re- 
peated time -scheduled observations of the parameters by very similar 
sensors from a relatively unchanging set of launch sites. At present, 
this category consists primarily of relatively small meteorological 
rockets with apogees of 60 to 100 km. Profiles of wind, pressure, 
temperature, and/or density are obtained. Synoptic observations are 
routinely reduced (usually at the observing site) and then collected 
by the National Climatic Center, Asheville, North Carolina. 


Research Rockets 

Research rockets are launched with a scientific payload aboard 
and obtain data that are used for the study of some aspect of the 
natural environment. The discipline areas include astronomy, solar 
physics, aeronomy, and energetic and thermal particles. By international 
agreement, research rocket launches are routinely reported to World 
Data Centers. A machine-sensible file, which contains selected infor- 
mation on these launches, is maintained at NSSDC. The information in 
this file is taken from the launch reports and from open literature. 

From this listing, launches can be selected by location, by time of 
day or date, by type of experiment, etc. The data must be obtained 
directly from the experimenter, but NSSDC can provide the experimenter's 
name and current address in most cases. A current copy of this listing 
or a uniquely sorted version of this listing can be obtained from NSSDC. 
Also available is a library file that lists publications that contain 
or have used rocket data. A very small number of engineering and 
development rockets are included in this file when the rocket type 
under development or testing is being prepared especially for subsequent 
use as a research rocket. 


Engineering and Development Rockets 

Engineering and development rockets are instrumented to obtain 
rocket performance data. These data are intended to improve or eval- 
uate rocket design. Occasionally, research payloads are included on 
this type of rocket. In such a case, a WDC attempts to include infor- 
mation about the flight in this file along with the research rockets. 
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PART IX - SOURCES 


This part contains brief descriptions of the principal sources of 
correlative data. Both data centers and publications are discxissed. 

A single alphabetically ordered list of documents containing correla- 
tive data or information on the availability of such data is also given. 
Sufficient detail is given to enable the reader to obtain the document. 
In general, no attempt has been made to include the publications of 
individual observatories either in this part or in the discipline- 
oriented parts of this Handbook . More extensive lists of data sources 
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can be found in: Report UAG-4 of the World Data Center A (WDC-A) for 

Upper Atmosphexe Geophysics, 1969 ("Abbreviated Calendar Record 1966- 
1967") ; Annals of the IQSY, volume 1, 1968; and the AFCRL Geophysics 
and Space Data Bulletin, volume IV, No. 2, 1967. 


lUCSTP, STP NOTES 


.The Inter-Union Commission on Solar Terrestrial Physics (lUCSTP) 
was established in 1966 to provide a single international organization 
with a comprehensive interest in solar-terrestrial physics. lUCSTP is 
composed . of five preexisting ICSU (International Council of Scientific 
Unions) organizations, each .interested in certain aspects of solar- 
terrestrial physics . These organizations are the International Union 
of Astronomy (I AU) , the.- International Union of Geodesy and Geophysics 
, (IUGG),,-;the' International Union of Radio Science (URSI) , the Interna- 
tional" Union of Pure’ and Applied Physics (lUPAP) , and the Committee on 
. Space Research (COSPAR) . ' 

Since April 1968, TUCSTP has published the STP Notes to provide 
solar-terrestrial physicists at the working level with information on 
relevant research that is being performed or proposed throughout the 
world. Recommendations of- lUCSTP and its individual working groups 
are published, as are data -.synopses in the form of "Abbreviated 
Calendar Records." (See Part VII.) STP Notes are published three to 
four times a year and are available through the lUCSTP Secretariat, 
c/o National Academy of Sciences, 2101 Constitution Avenue, N.W. , 
Washington, D. C. 20418. . 


WORLD DATA CENTERS 


World Data Centers (WDC) conduct international exchange of geo- 
physical observations in accordance with the principles set forth by 
the International Council of Scientific Unions. First established for 
the International Geophysical Year (IGY) , 1957 to 1958, they provide 
the means by which data gathered at individual observatories are made 
available to the international scientific community. 
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Three major centers exist: WDC-A in the U.S.A., WDC-B in the 
U.S.S.R. and WDC-C in Western Europe, Australia, and Japan. Each 
WDC is made up of a number of subcenters, each of which is devoted 
to one or more specific geophysical disciplines. The geographical 
location of each subcenter is usually different than that of the other 
subcenters. Certain of the subcenters have recently been designated 
as Solar Terrestrial Physics World Data Centers (STP-WDC) by the Inter- 
Union Commission on Solar Terrestrial 'Physics . 

Overall responsibility for WDC-A (coordination office andnine 
subcenters) resides with the U.S. National Academy of 'Sciences through 
the Geophysics Research Board (GRB) and its Committee on Data Exchange 
and Data Centers; The WDC-A subcenters related to the phenomena dis- 
cussed in this Handbook are : 

WDC-A for Upper Atmosphere Geophysics 
National Oceanic and Atmospheric Administration 
Boulder, Colorado, U.S. A. .80302 
Telephone: (303) 447-1000, Ext. '3181 

WDC-A for Geomagnetism, Seismology, and Gravity 
Environmental Data Service, NOAA 
Rockville, Maryland, U.S. A. 20852 
Telephone: (301) 496-8160 

WDC-A for Meteorology 
National Climatic Center, NOAA 
Asheville, North Carolina, U.S. A. 28801 
Telephone: (704) 254-0236 

WDC-A for Rockets and Satellites 
Goddard Space Flight Center, Code 601 
Greenbelt, Maryland, U.S. A. 20771 
Telephone: (301) 982-6695 

More extensive discussions of the WDC system can be found in IQSY 
Instruction Manual No. 6 or in the IQSY Annals, volume 1, "Geophysical 
Measurements, Techniques, Observational Schedules, and Treatment of 
Data." These references contain the "Guide to International Data Ex- 
change through the World Data Centers." Likewise STP Notes No. 6, 

"Guide for International Exchange of Data in Solar-Terrestrial Physics," 
contains a full discussion of the STP-WDCs. 

Individual WDC subcenters periodically issue catalogs listing their 
data holdings; any special publications involving data syntheses or 
special data compilations are also listed in the catalogs. Detailed 
information on the method of ordering data and the associated costs is 
presented. The catalogs are available upon request to the appropriate 
WDC subcenter. 
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Of principal interest to most users of this Handbook will be the 
WDC-A Catalogue Of Data on Solar -Terrestrial Physics , which has been 
issued annually during recent years. This data catalog lists the 
holdings of the WDC-A for Upper Atmosphere Geophysics, the WDC-A for 
Geomagnetism, and the WDC-A for Rockets and Satellites. The June 1970 
issue of the catalog is Report UAG-11 of the WDC-A for Upper Atmosphere 
Geophysics. The annual series of UAG documents can be obtained by 
subscription (see Data Compilations, p. 194), or a single copy of the 
catalog can be obtained from the Superintendent of Documents or from 
WDC-A at a cost of $1.50. 

For the convenience of the reader, the contents of the data catalog 
(1970) are given here in table 16. The numbering scheme shown not only 
reflects the organization of the catalog but also shows the discipline 
classification system as proposed by lUCSTP. 

In addition to the data listed, this WDC-A catalog contains an 
alphabetically ordered, cross -discipline station list that gives ob- 
serving program, station latitude and longitude (geographic and geomag- 
netic) , magnetic dip angle, L value, open and close dates, and, in 
appropriate cases, cosmic ray cutoff rigidity and altitude. 

Volume 7 of the Annals of the IQSY, "Sources and Availability of 
IQSY Data," contains lists of observing stations, rockets, satellites, 
and space probes active during part or all of the 1964 to 1965 period. 
This source also lists data available at the World Data Centers for the 
period 1957 to 1965. 


SOLAR-GEOPHYSICAL DATA 


The Institutes for Environmental Research (lER/NOAA) publish Solar- 
Geophysical Data (SGD) monthly in two parts. The first part is Solar- 
G eophysical Data (Prompt Reports), which contains data for the first 
and second months prior to publication; the second part is Solar -Ge ophys - 
ical Data (Comprehensive Reports) , which contains data for the sixth 
and seventh months prior to publication along with other miscellaneous 
information. A descriptive text is published annually. Prior to July 
1969, data for the first, second, and sixth months prior to publication 
appeared in the single monthly issue of Solar-Geophysical Data. Between 
November 1955 and November 1965, this series was the CRPL-FB-Solar Geo- 
physical Data, issued by the Central Radio Propagation Laboratory of the 
National Bureau of Standards. 


184 



TABLE 16 


CONTENTS OF THE WORLD DATA CENTER A 
CATALOGUE OF DATA ON SOLAR -TERRESTRIAL PHYSICS 


MASTER STATION LIST 

SOLAR AND INTERPLANETARY PHENOMENA 

A. I Sunspot Positions, Areas, and Classification 

A. 2 Sunspot Numbers 

A. 3 Solar Magnetic Fields 

A. 4 H-a Observations (other than flares) 

A. 5 Calcium Plages - Positions, Areas, Maximum Intensities 

A. 6 Combined and Special Optical Observations 

A. 7 Optical Observations of the Corona 

A. 8 Total Radio Flux Measurements 

A, 9 Radio Maps of Solar Disk 

A. 10 Radio East-West Scans of Solar Disk 

A. 11 Solar X-Ray and UV Background Levels 

A. 12 Energetic Solar Protons 

A. 13 Solar Wind 

A. 14 Comet Tails, Interplanetary Scintillations, Zodiacal Light 

A. 15 Sporadic Radio Emissions from Jupiter 

IONOSPHERIC PHENOMENA 

B. l Ionosphere Vertical Soundings 

B.2 Topside-Vertical Incidence Soundings and Satellite Probe Data 

B.3 Incoherent Scatter Soundings 

B.4 Oblique Incidence Soundings 

B.5 Ionospheric or Aeronomical Rockets 

B.6 Total Electron Content - Satellite Beacons 

B.7 Absorption - Method A 1 (pulse echo) 

B.8 Absorption - Method A 2 (riometer) 

B.9 Absorption - Method A 3 (CW field-strength) 

B.IO Ionospheric Drifts 

B.ll Ionospheric Scintillations from Satellite Beacons 
B,12 Ionospheric Back- and Forward-Scatter 
B.13 Whistlers and VLF Emissions 
B.14 Atmospheric Radio Noise 
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TABLE 16 (continued) 


FLARE -ASSOCIATED EVENTS 

C.l H-a Flares 

C.2 Solar Magnetic Field in Active Regions and Their Short-Term 
Changes 

C.3 Solar Radio Events, Fixed Frequency 
C.4 Solar Radio Spectrograms of Events 

C.5 Solar X-Ray Observations 

C.6 Sudden Ionospheric Disturbances - Ground-Based Observations 
C.7 Solar Protons - Direct Measurement 

C.8 Solar Protons - Riometer 

C.9 Solar Protons - Ionospheric Vertical Incidence Soundings 

C.IO Solar Protons and Electrons - VHF Forward Scatter 
C.ll Solar Protons - Other Types of Measurements 
C.12 Solar, Ionospheric, or Aeronomical Rockets Launched During 
an Event 

C. 13 Cosmic-Ray Ground Level Increases 
GEOMAGNETIC PHENOMENA 

D. l Standard and Rapid Run Measurements 

D.2 Magnetospheric Micropulsation Phenomena 

D. 3 Space Magnetism 

AURORA 

E. l All -Sky Cameras and Others 

E. 2 Satellite Measurements (particle precipitation and photometers) 

COSMIC RAYS 

F. l Neutron Monitors and Supermonitors 

F.2 Ionization Chambers 

F.3 Meson Telescope (cubical, crossed, narrow angle, and wide angle) 

F.4 Balloon Measurements 

F.5 Aircraft and Ship Measurements 

F. 6 Satellite Measurements 

AIRGLOW 

G. l Ground-Based Photometers 

G. 2 Satellite Measurements (photometers) 

MISCELLANY 

H. l Nocti lucent Clouds 

H.2 Meteorological Rockets 
H.3 Ozone 
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The Solar-Geophysical Data series is presently edited by Miss Hope 
J. Leighton tinder the supervision of Miss J. Virginia Lincoln of the 
Aeronomy and Space Data Center, NOAA, Boulder. Copies can be obtained 
on a data exchange basis from the World Data Center A for Upper Atmos- 
phere Geophysics. Copies can be purchased from the Superintendent of 
Documents, Government Printing Office, Washington, D.C. 20402. The 
annual subscription cost is $30.50 for both parts or $15.50 for one 
part. (Add $4.00 per part for foreign mailing.) The single issue price 
is $1.25 per part. 

Table 17, based on the index of Issue No. 304, December 1969, indi- 
cates the volume in which data of a certain type and for a given month 
can be found. Issue No. 303, for example, which was published in Novem- 
ber 1969, contains data for October 1969. 


GEOPHYSICS AND SPACE DATA BULLETIN 


The Air Force Cambridge Research Laboratories (AFCRL) , Bedford, 
Massachusetts, issues the Geophysics and Space Data Bulletin on a quar- 
terly basis. The bulletin contains AFCRL -sponsored data which are 
collected from sites throughout the world. 

The first four issues, encompassing 1964 data, constitute Volume 
I, No. 1-4. Volume II, with few exceptions, contains quarterly data 
for 1965; Volume III, 1966; Volume IV, 1967; and Volume V, 1968. Each 
of four appendices have been published in various issues as follows: 


Auxiliary Data Sources: Volume I, No. 1, 2; II, No. 1; 

III, No. 2; IV, No. 2. 

Trapped Radiation Model Environments: Volume II, No. 4; 

IV, No. 1. 

Annual Summary of Activity Cycle No. 20 Solar Flares : 

Volume IV, No. 3. 

Cumulative Index : Volume IV, No. 1, No. 4. 

Program descriptions for all data sections are updated annually or 
as required and are presented in Issue No. 1 of each volume. Each new 
section is described the first time it appears . 
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TABLE 17 


DATA PUBLISHED IN SOLAR-GEOPHYSICAL DATA 

1969 

Apr. May June July Aug. Sep. Oct. Nov. 


American Relative Sunspot 


Numbers R^ 

Ztlrich Provisional Relative 

297 

298 

299 

Sunspot Numbers R 2 
Zurich Final Sunspot 

297 

298 

299 

Numbers R 2 

2800 MHz - Daily Values of 

— 

“ “ — 

“ “ “ 

Solar Flux (ARO-Ottawa) 
2800 MHz - Daily Values of 
Adjusted Solar Flux (ARO- 

297 

298 

299 

Ottawa) 

8800, 4995, 2695, 1415, 606 
MHz Adjusted Solar Flux 

297 

298 

299 

(AFCRL) 

Mt. Wilson Magnetic Charac- 

297 

298 

299 

teristics of Sunspots 

298 

299 

300 

Mt. Wilson Magnetograms 
Ha (Sacramento-Peak or 

298 

299 

300 

Boulder) 

Sunspots (Boulder or 

298 

299 

300 

Sacramento Peak) 

Calcium Plage Maps (McMath- 

298 

299 

300 

Hulbert or Catania) 
Calcium Plage and Sunspot 

298 

299 

300 

Region 

298 

299 

300 

Coronal Line Emission 

298 

299 

300 

Optical Observations Flares 
Optical Observations Flares 
(Including Standardized 

297 

298 

299 

Data) 

302 

303 

304 

Flare Patrol Observations 

297 

298 

299 

Flare Patrol Observations 
Solar X-Ray Spectroheliograms 

302 

303 

304 

(OSO 5) 

Solar X-Ray Radiation 



" 

(Explorer 37) (Graphs) 
Solar X-Ray Radiation 

302 

303 

304 

(Explorer 37) 

Solar X-Ray Radiation 

298 

299 

300 

(Explorers 33 S 35) 

303 

304 

304 

Solar Protons (Explorer 34) 

303 

303 

304 


300 

301 

302 

303 

304 

300 

301 

302 

303 

304 

300 

301 

302 

303 

304 

300 

301 

302 

303 

304 

300 

301 

302 

303 

304 

301 

302 

303 

304 



301 

302 

303 

304 

— 

301 

302 

303 

304 

— 

301 

302 

303 

304 

— 

301 

302 

303 

304 

— 

301 

302 

303 

304 



301 

302 

303 

304 

— 

300 

301 

302 

303 

304 

300 

301 

302 

303 

304 

301 

302 

303 

304 

— 

301 

302 

303 

304 
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TABLE 17 (Continued) 


1969 



Apr. 

May 

June 

July 

Aug. 

Sep. 

Oct. 

Nov, 

Solar Proton Events 









(Riometer) ' 

Cosmic-Ray Protons (Pioneers 
6 8 7) 

Cosmic-Ray Protons (Pioneers 

297 

298 

299 

300 

301 

302 

303 

304 

8 8 9) 

Solar Wind (Pioneers 687) 

297 

298 

299 

300 

301 

302 

303 

304 

Solar Wind (Vela 2, 3, 8 4) 
Magnetic Field (Pioneer 9) 
Solar Radio Waves - Outstand- 

297 

298 

299 

300 

301 

302 

303 

304 

ing Occurrences - Worldwide 
Sudden Ionospheric Distur- 

302 

303 

304 

* 

“ “ 


"■ 


bance (SWF-SCNA-SEA-SPA- 
SES-SFD) 

297 

298 

299 

300 

301 

302 

303 

304 

Sudden Ionospheric Distur- 









bance (Addenda) 

302 

303 

304 

— 

— 

— 



— 


10700, 2700, 9,60, 328 MHz - 
Outstanding Occurrences 
(Pennsylvania State Uni- 
versity) 297 298 304 

15400, 8800, 4995, 2695, 

1415, and 606 MHz Out- 


standing Occurrences 
(AFCRL, Sagamore Hill) 297 

7000 MHz - Outstanding 

Occurrences (Sao Paulo) 297 

2800 MHz - Outstanding 

Occurrences (ARO- Ottawa) 297 

486 Mlz - Outstanding 
Occurrences (Washington 
State University) 297 

408 MHz - Outstanding 

Occurrences (San Miguel) 297 

408 MHz - Interferometric 

Observation (Nan^ay) 297 

169 MHz - Interferometric 
Observation (Nancay) 297 

184 MHz- Outstanding 
Occurrences (Boulder) 297 

18 MHz T Bursts (Boulder) 297 

18 MHz - Bursts (McMath) 297 

30-1000 MHz - (Weissenau, 

G.F.R.) 298 


298 304 

298 304 

298 304 


298 

304 

— 

— 

— 

— 

— 

298 

299 

300 

301 

302 

303 

304 

298 

299 

300 

301 

302 

303 

304 

298 

304 



— 

— 

— 

— 

298 

304 

— 

— 

— 

• 

— 

298 

304 

— 

— 

— 

— 

— 

299 

300 

301 

302 

303 

304 

— 
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TABLE 17 (Continued) 


1969 



Apr. 

May 

June 

July 

Aug. 

Sep. 

Oct. 

Nov 

10-580 MHz - (Fort Davis) 

298 

299 

300 

301 

302 

303 

304 

^ M » 

10-210 MHz - (Culgoora)* 

298 

299 

300 

301 

302 

303 

304 

— 

7.6-80 MHz - (University 

of Colorado)** 

298 

299 

300 

301 

302 

303 

304 

— 

19-41 MHz - (AFCRL, 

Sagamore Hill) 

298 

299 

300 

301 

302 

303 

304 

— 

9.1-cm Spectroheliograms 

(Stanford) 

298 

299 

300 

301 

302 

303 

304 

— 

21-cm Spectroheliograms 
(Fleurs) 

298 

299 

300 

301 

302 

303 

304 


10.7-cm Solar Scans 

(Ottawa-ARO) 

297 

298 

299 

300 

301 

302 

303 

304 

21-cm Solar Scans (Fleurs) 

297 

298 

299 

300 

301 

302 

303 

304 

43-cm Solar Scans (Fleurs) 

297 

298 

299 

300 

301 

302 

303 

304 

Cosmic- Ray Neutron Counts 

(Deep River) 

298 

299 

300 

301 

302 

303 

304 

— 

Cosmic- Ray Neutron Counts 

(Alert) 

298 

299 

300 

301 

302 

303 

304 

— 

Cosmic-Ray Neutron Counts 

(Churchill) 

298 

299 

300 

301 

302 

303 

304 

— 

Cosmic- Ray Neutron Counts 

(Climax) 

298 

299 

300 

301 

302 

303 

304 

— 

Cosmic-Ray Neutron Counts 

(Dallas) 

298 

299 

300 

301 

302 

303 

304 

— 

Geomagnetic Indices Ci, Cp, 
Kp, Ap - Selected Days 

298 

299 

300 

301 

302 

303 

304 

... 

Principal Magnetic Storms 

298 

299 

300 

301 

302 

303 

304 

— 

Reduced Magnetograms 

302 

303 

304 

— 

— 

— 

— 

— 

Sudden Commencement and 

Solar Flare Effects 

303 

303 

303 

304 

304 

304 

— 

— 

27-Day Chart of Kp Indices 

for Year 

— 

— 

— 

— 

— 

— 

— 

--- 

27-Day Chart of Cg for Year 

— 

— 

— 

— 

— 

— 

— 

— 


High Latitude Quality 


Figures and Forecasts 

298 

299 

300 

301 

302 

303 

304 

High Latitude Comparison 

Graphs 

298 

299 

300 

301 

302 

303 

304 

Graphs of Transmission 

Frequency Range 

298 

299 

300 

301 

302 

303 

304 

lUWDS Alert Decisions 

297 

298 

299 

300 

301 

302 

303 

Abbreviated Calendar Record 
*8-2000 MHz beginning August 

303 

1969. 

304 







**7.6-41 MHz beginning September 1968. 
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The bulletin lists the data by category and by station name and 
date. The following categories of data are included: 

• Total Field Intensity Magnetometer 

• Cosmic-Ray Neutron Intensity Monitor 

• Extremely Low Frequency Noise 

• Riometer 

Hourly Absorption Values 

• Solar Optical Patrol Observations 
Sunspot 

Chromospheric Flares 
Flare Film Record 

Chromospheric Activity Other Than Flares 
White-Light Visual Simspot Drawing 
Spectrographic Coronal Observations 

Sunspot Group Analysis (White-Light Photographic Patrol) 
Solar Event Atlas 

• Solar Radio Emission 
Spectral Observations 
Radio Noise Measurements 
Radio Map 

• Vertical Incidence Ionospheric Soundings 
Absorption Loss 

Median Values 

Median Curves and Hourly Values 

Copies of the Geophysics and Space Data Bulletin may be requested from: 

Miss Anne L. Carrigan (editor) 

Space Physics Laboratory 

Air Force Cambridge Research Laboratories 

L. G. Hanscom Field 

Bedford, Massachusetts 01730 
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lAGA BULLETINS NO. 12 AND NO. 18 


The International Association of Geomagnetism and Aeronomy (lAGA) 
of the International Union of Geodesy and Geophysics (lUGG) publishes 
bulletins that contain, in some cases, transactions of meetings and, in 
other cases, listings of geomagnetic data. Copies of the bulletins can 
be purchased from the lUGG Publication Office, 39ter, Rue Gay-Lussac, 
Paris (V) , France . 

In 1948, Bulletin No. 12 was issued to present compilations of the 
geomagnetic indices, C and K, that were received from participating 
observatories from 1940 to 1946. Since the publication of Bulletin No. 
12, annual supplementary bulletins (No. 12a to 12-£) have been published 
toT present data that were collected between 1947 and 1957. ( Bulletin 
No. 12d is an exception.) In each of subsequent years, two bulletins 
have been published, e.g., 12ml and 12m2. Typically, data for the 
third year prior to publication are presented. 

An example of the contents of the I AGA Bulletin No. 12 series, 
supplements 12ol and 12o2, is given below: 

(12ol) 


List of observatories 
Observatories reporting K-indices 
Three-hour -range indices K, 1960 
Magnetic character-figures C, 1960 
International character-figures Ci, 1960 
Magnetically selected days 

International character-figures Ci, 1900-1960 
Planetary indices Kp, equivalent ranges ap, daily average 
ranges Ap, daily character-figures Cp, 1960 
Frequencies of Kp-indices, 1960 
Monthly averages of Ap and Cp, 1960 
Very quiet intervals 
27-day recurrence for Kp, 1960 

References to tables and diagrams for Kp, ap, and Cp 
(12o2) 

List of observatories 

Storm sudden commencement, ssc, 1960 

Bays and pulsational disturbances, b, bs, bp, bps, pt, pg, 1960 
Pulsational disturbances, 1960 (pt and pg) not associated 
with bays 
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Sudden impulses, s.i., 1960 
Minor disturbances, 1960 
Solar-flare effects 
Doubtful solar-flare effects, 1960 
Rejected solar-flare effects, 1960 

lAGA Bulletin No. 18 , "Geomagnetic Planetary Indices Kp, Ap, and 
Cp, 1932 to 1961" contains, in addition to the data indicated by its 
title, values of ap for the periods 1932 to 1936 and 1952 to 1958, and 
27-day recurrence diagrams for Kp for the period 1932 to 1961. These 
diagrams are also referred to as Bartel's musical diagrams. 


QUARTERLY BULLETIN ON SOLAR ACTIVITY 


The Quarterly Bulletin on Solar Activity (QBSA) is published for 
the International Astronomical Union by the Eidgen,. Sternwarte in 
Zurich, Switzerland, .with financial assistance from UNESCO. Since 
only a limited number of the bulletins are printed, only one copy is 
sent to a single institution. The data contained in an issue nominally 
were taken 1-1/2 years prior to publication. The QBSA is divided into 
four parts: Part 1, "Simspots," written in English, contains sunspot 

numbers and areas; Part .2, "Eruptions Chromospheriqiies Brillantes" 
(solar flares), in French, contains flare listings, observing hours, 
and a table of active regions; Part 3, "Intensities de la Couronne 
Solaire" (coronal emissions), in French, contains listings of coronal 
line emission intensities.; and Part 4, "Solar Radio Emissions," in 
English, contains tables of flux values, a list of distinctive events, 
maps showing the location of active regions, and a list of the spectral 
class of activity. These data have been more complete than those pub- 
lished in Solar-Geophysical Data. For instance, QBSA has contained 
worldwide radio data, whereas, prior to mid-1969, SGD contained only 
Western Hemisphere radio data. The . Quarterly Bulletins for January 
1917 to June 1968 are available on four rolls of microfilm from WDC-A, 
Boulder. 
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DATA COMPILATIONS FOR INDIVIDUAL EVENTS 


A number of data compilations related to individual, outstanding 
solar -terrestrial events are available in published form. 

• Solar Event of May 23, 1967 

"Data on Solar Event of May 23, 1967, and its Geophysical Effects," 
WDC-A for Upper Atmosphere Geophysics Report UAG-5. 

• October 24-November 6, 1968, Solar Activity 

"Data on Solar Geophysical Activity October 24-November 6, 1968," 
WDC-A for Upper Atmosphere Geophysics Report UAG-8, two parts. 

• Cosmic-Ray Event, November 18, 1968 

"Data on Cosmic-Ray Event of November 18, 1968, and Associated 
Phenomena," WDC-A for Upper Atmosphere Geophysics Report UAG-9. 

These reports were compiled by Miss J. Virginia Lincoln and contain 
discussions by principal investigators in addition to the actual 
data. Published by the Environmental Data Service of NOAA at 
Asheville, North Carolina, the documents can be purchased from the 
Superintendent of Documents, Government Printing Office, Washington, 
D.C. 20402. The single issue price varies: for UAG-5 the price is 

$0.65; for UAG-8, $1.75; and for UAG-9, $0.55. The annual subscrip- 
tion rate for the l^per Atmosphere Geophysics series is $9.00 in 
the U.S. and $11.50 for subscribers outside the U.S. The number 
of issues varies from year to year. 

• Solar Flare of July 7, 1966 

The solar flare of July 7, 1966, and its geophysical effects is 
the subject of Volume 3 of the Annals of the IQSY . Both the solar 
flare data and discussions by the principal investigators are pre- 
sented. 

• Solar Event of February 25, 1969 

Another potentially useful document is "Ground Based Observations 
for the Solar Event of February 25, 1969," compiled by K. G. 
-Lenhart, European Space Operations Center, Darmstadt, Federal 
Republic of Germany. Solar, ionospheric, geomagnetic, and neu- 
tron monitor data from several observatories are presented. 
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MISCELLANEOUS DOCUMENTS 


Infrequent references are made in this Handbook to publications 
that may not be readily available to most readers. This section lists 
these publications along with the publisher and/or center of distribu- 
tion. The publications are listed here alphabetically by title. 

Aerdldgical Data of Japan, published by the Japan Meteorological Agency, 
Tokyo, available from the National Climatic Center, Asheville, North 
Carolina 28801. 

"Data Report of the Meteorological Rocket Network Firings," published 
by IRIG. Copies issued prior to 1966 are available from the Defense 
Documentation Center, Cameron Station, Alexandria, Virginia 22314. 
Subscriptions for current editions are available from the Superinten- 
dent of Documents, Government Printing Office, Washington, D.C. 20402. 

Environmental Science Services Administration Technical Reports WB-2, 

3, 9, 12 (1964-1967), available from the Superintendent of Documents, 
U.S. Government Printing Office, Washington, D.C. 20402. 

EXAMETNET Data Report Series, Annual Report for 1966 (NASA SP-175), An- 
nual Report for 1967 (NASA SP-176), and Annual Report for 1968 (NASA 
SP-231), available from the National Technical Information Service, 
Springfield, Virginia 22151. 

"Format of Electron Density Profile Tabulations," prepared by J. W. 
Wright and D. McKinnis, available from Institute for Telecommunication 
Sciences and Aeronomy, NOAA, Boulder, Colorado 80302. 

"Handbook of Geophysics and Space Environments," edited by Shea L. 
Valley, available from Air Force Cambridge Research Laboratories, L. G. 
Hanscom Field, Bedford, Massachusetts 01730. 

"International Geophysical Calendars 1957-1969," published by WDC-A for 
Upper Atmosphere Geophysics, available from the Superintendent of Docu- 
ments, Government Printing Office, Washington, D.C. 20402 ($.30). 

International Ursigram and World Days Service Synoptic Codes for Solar 
and Geophysical Data, available from Miss J. V. Lincoln, Deputy Secre- 
tary, lUWDS Steering Committee, NOAA, Boulder, Colorado 80302. 

Key to Meteorological Records (KMRD) Series, published by the National 
Oceanic and Atmospheric Administration, available from the Superinten- 
dent of Documents, Government Printing Office, Washington, D.C. 20402. 
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"North Polar, South Polar, World Maps and Tables of Invariant Magnetic 
Coordinates for Six Altitudes: 0, 100, 200, 600, 1000, and 3000 km," 
published by DASA, available from the authors, J. E. Evans, L. L. 
Newkirk, . and B. M. McCormac, Lockheed Palo-Alto Research -Laboratory, 
3251 Hanover Street, Palo Alto, California 94304. - 

Ozone Data for the World, published bimonthly by the Meteorological 
Branch, Canadian Department of Transport, Toronto, Ontario, available 
from the National Climatic Center, Asheville, North Carolina 28801. 

Photographic Journal of the Sun, M. Cimino, editor, published by ' 
Osservatorio Astronomico di Roma, Italia. 

"Preliminary Report and Forecast of- Solar-Geophysical Activity," pub- 
lished by the National Oceanic and Atmospheric Administration, Boulder, 
Colorado 80302 . : : 

"Royal Observatoiy Bulletins," published in Photoheliographic Results, 
available from WDC-A for Upper Atmosphere Geophysics, Boulder, Colorado 
80302. 

"The American Ephemeris and Nautical Almanac" published by U.S. Naval 
Observatory, available from the Superintendent of Documents i Government 
Printing Office, Washington, D.C. 20402. ($4.25). 

"Three-monthly Bulletin on Geomagnetic Indices and Rapid Variations," 
International Service of Geomagnetic Indices, c/o Royal Netherlands 
Meteorological Institute, DeBilt, The Netherlands. 

USSR Academy of Sciences Solar Data , available from Astronomical Obser- 
vatory of the U.S.S.R. Academy of Sciences, Leningrad M-140, U.S.S.R. 

"U.S. Standard Atmosphere, 1962" and "U.S. Atmosphere Supplements," 
published jointly by the Environmental Science Services Administration, 
the National Aeronautics and Space Administration, and the U.S. Air 
Force, available from the Superintendent of Documents, Government 
Printing Office, Washington, D.C. 20402. 
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LIST OF ABBREVUTIONS 


A 

AAFB 

ADALERT 

AFB 

AFCRL 

ATS 

BCD 

Bev 

C5GS 

CIRA 

cm 

COESA 

COSPAR 

Dst 

ELF 

env 

ESRO 

ESSA 

ETAC 

euv 

EXAMETNET 

GEOALERT 

GMT 

GRB 

GSFC 

HAO 

hr 

lAGA 

lAU 

IBM, 

ICAO 

ICSU 

IGRF 

IGY 

IQSY 

IRIG 

lUCSTP 

lUGG 

lUPAP 

lUWDS 

K 

kHz 

km 

KMRD 

mb 

Mev 


angstrom 

Auxiliary Air Force Base . 

Advance Alert 
Air Force Base 

Air Force Cambridge Research Laboratories 
Applications Technology Satellite 
binary coded decimal 
billion electron volts 

Coast and Geodetic Survey > 

COSPAR International Reference Atmosphere 

centimeter 

Commission on Extension to the Standard Atmosphere 
Committee on Space Research 

axisymmetric component of geomagnetic disturbance field 

extremely low frequency . . 

environment 

European Space Research Organization 
See NOAA \ 

Environmental Technical Applications Center 
extreme ultraviolet radiation 

Experimental InterAmerican Meteorological Rocket Network 

Geomagnetic Alert 

Greenwich mean time 

Geophysics Research Board 

Goddard Space Flight Center 

High Altitude Observatory 

hour 

International Association of Geomagnetism and Aeronomy 
International Astronomical Union 
International Business Machines 
International Civil Aviation Organization 
International Council of Scientific Unions 
International Geomagnetic Reference Field 
International Geophysical Year 
International Years of the Quiet Sun. 

Inter- Range Instrumentation Group . 

Inter-Union Commission on Solar-Terrestrial Physics 

International Union of Geodesy and Geophysics 

International Union of Pure and Applied Physics 

International Ursigram and World Days Service 

Kelvin 

kilohertz 

kilometer 

Key to Meteorological Records Documentation 
millibar 

million electron volts 
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MHz 

min 

MIT 

mm 

MRN 

MSC 

MSFC 

MWG 

NASA 

NOAA 

NSSDC 

OGO 

Pc 

Pg 

Pi 

pibal 

POGO 

Pt 

QBSA 

rms 

RWC 

sc, ssc 
SGD 
s .i . 

SID 

Sp 

SP 

SPME 

STP-WDC 

UNESCO 

URSI 

USAF 

USCGS 

UT 

VHF 

VLF 

vs 

WDC-A 

W.M.S. 

WWA 


megahertz 

minute 

Massachusetts Institute of Technology 
millimeter 

Meteorological Rocket Network 
Manned Spacecraft Center 
Marshall Space Flight Center 
Meteorological Working Group 
National Aeronautics and Space Administration 
National Oceanic and Atmospheric Administration, 
Environmental Science Services Administration 
North Pole 

National Space Science Data Center 
Orbiting Geophysical Observatory 
continuous pulsation 
giant pulsation 
irregular pulsation 
pilot balloon 

Polar Orbiting Geophysical Observatory 
trains of pulsation 
Quarterly Bulletin on Solar Activity 
root mean square 
Regional Warning Center 
sudden storm commencement 
Solar-Geophysical Data 
sudden impulse 

sudden ionospheric disturbance 
South Pole 
special publication 
Solar Proton Monitoring Experiment 
Solar Terrestrial Physics World Data Center 
United Nations Economic, Scientific and Cultural 
Organization 

International Union of Radio Science 

United States Air Force 

United States Coast and Geodetic Survey . 

universal time 

very high frequency 

very low frequency 

versus 

World Data Center A 
World Magnetic Survey . 

World Warning Agency 


formerly 


198 



INDEX 


A 


Page 


Page 


a, B coordinate system 170, 174 

AA 125 

Abbreviated Calendar Record 163, 190 

Absorption 115-118, 125-128, 185 

Activity indices 76 

ADALERT 175 

AE 84, 85 

Aeronomical rocket 185, 186 

Aeronomy 177 

Aerospace Corporation 68 

Airglow 160, 186 

Air showers 7, 9 

Ak 80 

ak 80 

Alaska, U. of 84 

All-sky camera 161-163, 184 

Alosyn 113 

Alouette 1 111-113, 130 

Alouette 2 Ill, 112, 130 

Alpha particle 66 

American Relative Sunspot Number ... 45, 186 

Anemometer 138 

Ap 81, 82, 152, 164, 

165, 190, 192 

ap 80-82, 192 

Applications Technology Satellites .... 144 

Arc 162 

Ascaplot 162 

A-scan 104, 116 
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